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UNCLASSIFIED ABSTRACT 


\ i 4 
(u) “this technical report describes analyses and techniques used I 
dn the design and evaluation of advanced decoy concepts. The work _ 
described addresses both the design of specific penetration aid | a 
elements and the formulation of techniques for their evaluation.) \ ; | 
The three major technical areas covered in this report are: Te ee 


= —— 


lL. Investigation of a penetration aid technique that degrades 
the measurement capability of the radar sensor. 


2. The design of a computer program to solve the decoy design 
problem with flexibility in the selection of optimization 
eriteria and constraints. 


3, Studies of the use of certain discrimination techniques 
for a hard point defense system. 


This appendix to this report contains detailed description of the 
optimum decoy design program. 
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4, Numeric al Procedure 


First the coefficient counter M and the accumulator ANS are 
initialized to zero. ThenITPis tested. If ITP is greater than one, 
statement 10 is executed next, setting NAZU equal to NAX and control 
then goes to statement 30. In this case, a polynomial in the three 
unknowns X, Y and Z will be evaluated. If ITP is less than or equal to 
1, statemert 20 is executed next where NAZU is set equal to 1 and 
statement 30 is reached, In this case, the polynomial will be a tunction 


of only two unknowns X and Y. 


At staternent 30, the DOloop ending at statement 130is entered to 
perform the calculations needed to evaluate the polynomial. The value 
of this polynmial is returned to MISC in the variable ANS. When the 
evaluation is complete, the RETURN statement following 130 is executed, 


returning control to MISC, 


5. Other Information 


A. P@LCAL is called by SUBROUTINE MISC only. 


B. P@LCAL does not call or reference any other subprogram. 
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Comparison of the Decoy with Reentry Vehicle 


The comparison of the decoy with the reentry vehicle consists of the 


evaluations of certain differences and the determination of integrals 
of special functions. SUBROUTINE F123 performs the difference 


calculations using SAVEDV to store the values, before calling INTGRL. 


SUBROUTINE INTGRL evaluates the integrals using in the process 
FUNCTION ADD, and subroutines INFCOF, INTERP, LINFIT, and 


i WRITEM, F123 also calls in subroutines HEADER AND TIMERS. 


These integrals include the corridor functions user to constrain the 


decoy within specified corridors and the effectiveness integrals which 


are used in the probability of discrimination calculations in SUBROU- 


TINE EFFECT. 
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SUBROUTINE F123(XXX) 


h | H 1, Purpose 
| F123 can optionally 


ti ht a. Read the data for a reentry vehicle 'cajectory from input or call 
| SUBROUTINE: VIXEN to calculate this trajectory data, 


Call VIXEN to calculate the data for a decoy trajectory. 


Compare the data for a reentry vehicle trajectory with the data 


for a decoy trajectory and save tables of differences, 


Call INTGRL to calculate and plot a corridor integral. 


Call INTGRL to calculate the sigma square integral to be used by 


SUBR®UTINE EFFECT, 
CALL INTGRL to calculate a table of influence coefficients and 


produce influence coefficient plots at specified altitudes. 
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2, Input 


Common Source of 
Wame Location Input Description 


*YNOTE* The numbers in the Common Location column refer to positions in the 
@CCUR array. 


| ALST 122 READIT Stopping value of & envelope, in degrees 
CASE 128 READIT Case number 
DATE 127 READIT Date 
DELY 219 READIT Linear thrust offset along Y axis in inches 
DELZ 220 READIT Linear thrust offset along Z axis in inches 
be Sp EM 129 READIT Memo number 
GAMF0 105 READIT Initial flight path angle in degrees 
. ca oe ie af Ratio of specific heats 
LAI 138 READIT Axial length of initial configuration, in 
inches . 
LA2 144 READIT Axial length of vehicle after shape change 
at ZTURN, in inches 
| LP 4000 VIXEN Error code 
as sa pate se Molecular weight of air, gram/mole 
PHIO 112 READIT Initial Euler angle > , in degrees 
PSIZET 223 READIT Thrust misalignment angle, in degrees 
PSIO 114 READIT Initial Euler angle , in degrees 
RB1 136 READIT Base radius of inital configuration, in 
inches 
ry 
; 1-459 
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2, Input (cont'd) 


ommon Source of 


Name Location Input Description 
eae aren loc tee 
RB2 142 READIT Base radiua of vehicle after shape change 
at ZTURN, in inches 
RN1 135 READIT Nose radius of initial configuration, in 
inches 
RN2 14] READIT Nose radius of vehicle after shape change, 
at ZTURN, in inches 
aS ue a oy Collision cross section in angstroms 
THEALO 113 READIT Initial Euler angle ©, , in degrees 
THEZET 224 READIT Thrust misalignment angle in degrees 
*HHNOTE He Unless otherwise indicated, the entries in the Common Location 
Column refer to positions in the PCCUR array, 
BCB, 40 6181-6220 ZREADX Array of lower corridor values for 
ballistic coefficient 
BCD, 40 6101-6140 ZREADX Array of lower corridor values for 
é deceleration in g's 
BCV, 40 6021-6060 ZREADX Array of lower corridor values for 
velocity in ft/sec, 
BCWLI, 40 6501-6540 ZREADX Array of lower corridor valuee for wake 
i length in meters for first radar frequenzy 
BCWL2, 40 6381-6620 ZREADX Array of lower corridor values for wake 
length in meters for second radar frequency 
BCWL3,40 6661-6700 ZREADX Array of lower corridor values for wake 
length in meters for third radar frequency 
BCWRI, 40 6261-6300 ZREADX Array of lower corridor values for wake 
radar crosa section at first radar frequency 
BCWR2, 40 6341-6380 ZREADX Array of lower corridor values for wake 
radar cross section at second radar 
1-460 frequency 
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fos 
q Input (cont'd) 
q 
* Mia ee eS a ee oe eee eee ee 
Common Source of 
Name Location Input Description 
1 BCWR3, 40 6421-6460 ZREADX Array of lower corridor values 
| for wake radar cross section at 
third radar frequency 
BETAPL, 160 641 - 800 ZREADX or Ballistic coefficients for the input 
1 VIXEN or calculated trajectory points 
DVAR, 50 11707-11756 SAVEDV Array in which design characteristics 
| are saved 
| DVH, 50 11607-11656 ZREADX Input values of design variables for 
che second perturbation of com- 
parison decoy 
DVL, 50 11657-11706 ZREADX Input values of design variables for 
. the first perturbation of comparison 
decoy 
H ,40 5841-5880 ZREADX Input altitudes corresponding to 
corridor points,in feet 
* 
SB, 40 6821-6860 ZREADX Standard deviation of radar measure- 
ment errors for ballistic coefficient 
SD, 40 6781-6820 ZREADX Standard deviation of radar measure- 
ment errors for deceleration in g'a 
SV, 40 6741-6780 ZREADX Standard deviation of radar meaaure- 
ment errors for velocity in ft/sec, 
SWL1, 40 6981-7020 ZREADX Standard deviation of radar measure- 
ment errors for wake length in 
meters at first radar frequency 
SWL2, 40 7021-7060 ZREADX Standard deviation of radar measure- 
ment errors for wake length in meters 
at second radar frequency 
SWL3, 40 7061-7100 ZREADX Standard deviation of radar measure - 
ment errors for wake length in 
meters at third radar frequency, 
7 
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i 
| 
| i P : 
y Input (cont'd) |! 
fl | Common Source of ‘ 
i Name Location Input Description 
TPLOT, 160 1-160 ZREADX or Times for the input or calculated 
VIXEN trajectory points in seconds 
i VOGPLT, 160 481-640 ZREADX or Deceleration in g's for input or 
| VIXEN calculated traje ctory points 
| VPLOT, 160 321-480 ZREADX or Velocities in ft/sec, for the input 
VIXEN or calculated trajectar y points 
| WLIP, 160 1281-1440 ZREADX or —_ Wake length for the first radar i 
VIXEN frequency input or calculated i 
trajectory points in meters 
WL2P, 160 1441-1600 ZREADX or Wake length for the second radar 
VIXEN frequency input or calculated 
trajectory points in meters 
WL3P, 160 1601-1760 ZREADX or Wake length for the third radar 
VIXEN frequency input or calculated ' 
trajectory points in meters | 
WRIP, 160 801-960 ZREADX or Wake cross section for the first | 
; VIXEN radar frequency for input or 
calculated trajectory points, see 
IDBL for units 
WR2P, 160 961-1120 ZREADX or Wake cross section for the second i 
VIXEN radar frequency for input or cal- if 
culated trajectory points ii 
WR3P, 160 1121-1280 ZREADX or Wake cross section for the third : 
VIXEN radar frequency for input or cal- 
culated trajectory points 
ww anenen LINFIT Coefficient of the independent 
variable in a linear fit 
xx woneae LINFIT Constant in a linear fit 
ZPLOT, 160 161-320 ZREADX or _— Altitudes for the input or calculated 
! VIXEN trajectory points, feet 
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i 
i 
* i 
= 
Input (cont'd) i 
Common Source of \ 
: Name Location Input Description i 
SWRI, 40 6861-6900 ZREADX Standard deviation of radar measure- _ 
ment errors for wake cross section 
at first radar frequency 1 
SWR2, 40 6901-6940 ZREADX Standard deviation of radar measure- by 
ment errors for wake cross section [ | 
at second radar frequency r 4 
: ; 
SWR2, 40 6941-6986 ZREADX Standard deviation of radar measure- _ 
ment errors for wake cross section 
at third radar frequency hf | 
TCB, 40 6221-6260 ZREADX Array of upper corridor values for | | 
ballistic coefficient : 
TCD, 40 6141-6180 ZREADX Array of upper corridor values for } 4 
‘ deceleration in g's | 
TCV, 40 6061-6100 ZREADX Array of upper corridor values for 
velocity in ft/sec, 
,) a 
; TCWLI, 40 6541-6580 ZREADX Array of upper corridor values for 4 
wake length in meters for first 4 
radar frequency Pei 
TCWL2, 40 6621-6660 ZREADX Array of upper corridor values for mt 
wake length in meters for second Pl 
radar frequency a 
TCWL3, 40 6701-6740 ZREADX Array of upper corridor values for a 
wake length in meters for third Bi 
radar frequency ‘i. : 
fiw 
TCWRI, 40 6301-6340 ZREADX Array of upper corridor values for || | 
wake radar cross séction at first Ni 
radar frequency hd 
TCWR2, 40 6381- 6420 ZREADX Array of upper corridor values for 
wake radar cross section at second 
radar frequency 
TCWR3, 40 6461-6500 ZREADX Array of upper corridor values for 
. wake radar cross section at third 
radar frequency 
1-463 
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Input (cont'd) 
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Name 


#eo NOTE The following quantities are integers: 


IREF 
IMP, 90 


LOPT 


LPLOT 


MHEAT 
MODE 


MOPT 


NCDMDY, 50 


NCP 
NDECOY 


NDVCH 


NOSEDP 


NTHRST 


Coinmon 
Location 


IDCCUR(301) 


IDCCUR(1-90) 


NOCCUR(07) 


IPCCUR(302) 


NOCCUR(10) 


IDCCUR(303) 


NOCCUR(03) 


JOCCUR(91-140) 


IDCCUR(304) 
IDCCUR(305) 


- IMCCUR(306) 


_ NOCCUR(05) 


NOCCUR(22) 


Source of 
Input 
ZREADX 
ZREADX 
READIT 


ZREADX or 
VIXEN 


READIT 


ZREADX 
READIT 


ZREADX 


ZREADX 


ZREADX 


ZREADX 


READIT 


READIT 


T~464 


Description 


Trajectory processing option code, 
see ZREADX 


Input integer codes, see users 
manual 


Trajectory option code, see READIT 


Number of points in the input or 
calculated trajectory 


Option code for mass loss calculation 


Fundamental option code, see 
ZREADX 


Option code for aerodynamic 
heating calculation 


Identification code numbers (indices 
in the OCCUR array) of design 
variables to be perturbed under 
MODE =2 influence coefficent 
calculations 


Number of corridor points 


Vehicle type identification code, see 
ZREADX 


Number of entries in NCDMDV 
table, i.e,, number of design 
variables 


Option code for noseblunting cal- 
culation, sce READIT 


Option code for thrusting, see 
READIT 
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CDB 
CDFINF 
CDI 
CDP 


CDP90 


CM 


CN 
CODLAM 


CPSZET 


CTHZET 
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Common | 
Name Location Description i} 
Ms 

*k*XNOTE*** Numbers in Common Location Column refer to positions in i 
@CCUR array for the following quantities: i 
| i 
ALMAX, 200 1044-1243 Maximum in e@at the Vixen point altitudes in } 
radiug i 
ALMIN, 200 1244-1443 Minimum in o< at the Vixen point altitudes in i 
radius | 

ALST 122 The stepping control on envelope @ , preset 1 | 
to 0,2 degrees | | | 

it 

CASE 128 Case number th 


99 
2793-2800 
100 


98 


215 


| Ie eM oneal lo gett 
ID Laer wna 
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Base drag coefficient 

Skin friction drag coefficient 
Induced drag coefficient 

Pressure drag coefficient including 


Pressure drag coefficient for zero angle of 
attack 


Moment coefficient about c.g., Ce 


Normal force coefficient, C, 


17,6489 x 10° io, used in PRELIM to 
define CODRAG which is then used in DRAGCO 


Cosine of the thrust offset angle ¥ 
PSIZET $ 


Cosine of the thrust offset angle g 
THEZET $ 
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Common 

Name Location Description 

DECDFP 235 The pressure induced skin friction drag 
coefficient 

DECFTC 236 The transverse curvature induced skin 
friction drag coefficient 

DELCDP 234 The induced pressure drag coefficient 

DELY 219 The linear thrust offset component along 
Y axis, in feet 

DELZ 220 Linear thrust offset along Z axis, in feet 

FACTRY 197 Numerical factor used in DRAGCO 

FMAX, 200 2244-2443 Frequencies at the altitudes corresponding 
to ALMAX values, cycles/sec 

FMIN, 200 2444-2643 Frequencies at the altitudes corresponding 

wei to ALMIN values, cycles/sec 

GAMFOo 105 Initial flight path angle in radians 

LAI 138 See input, in feet 

LA2 144 See input, in feet 

LP 4000 Error code 

MDOT, 32 2708-2739 Distribution along the body of ablative mass 
loss rates, in lbm/ft -sec 

MX 210 thrusting moment about X axis 

MY 211 thrusting moment about Y axis 

MZ ' 232 thrusting moment about Z axis 

OCCUR, 4000 1+4000 Common array 
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Output (cont'd) 


Name 


PEPSB, 8 


PHIO 
PSI0 


SINGO 


SMRO 


SPSZET 
STHZET 


TAMAK, 200 


TAMIN, 200 


un spgengrnerenerer eter 


Common 
Location 


Tonle seni Tuer neces nse aemenenaaemneeneenenll 


2801-2808 


112 
114 


2676-2707 


110 


60 


136 
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Description 
Distribution along body of the ratio of 
edge pressure to stagnation pressure 
dimensionless 
Euler angle d , in radians 


Euler angle p , in radians 


Distribution along the body of cold wall 
aerodynamic heating rates 


Initial value of the angular rate Q 
in radians/sec. 


Time rate of change of base radius due to 
ablation, ft/sec 


See input, in feet 
See input, in feet 


Time rate of change of nose radius due to 
ablation, ft/sec, 


See input, in feet 
See input, in feet 
Sine of initial flight path angle 


Initial value of the angular rate R in 
radians /sec 


Sine of input PSIZET 
Sine of input THEZET 


Times corresponding to maximums in 
o{ , in seconds 


Times corresponding to minimums in 
a in seconds 
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tt Output (cont'd) 
1 f = . 
aa} Common 
i | 1 Name Location Description 
i | TBMAT, 3 3643-3645 Components of thrust in body axis system, 
lb 
| } THEALO 113 Initial Euler angle 9. ,in radians 
Ki WDOT 86 Rate of change of vehicle weight due to 
| ablation, lbm/sec 
ZMAK, 200 1844-2043 Altitudes corresponding to maximums in 
o , in feet 
ZMIN, 200 2044-2243 Altitudes correspondint to minimums in 
o , in feet 
***NOTE*** Unless otherwise designated, numbers in Common Location column 
refer to positions in the PCCUR array. ‘ 
HCB, 40 6181-6220 Array of lower corridor values for 
ballistic coefficients 
BCD, 40 6101-6140 | Array of lower corridor values for ‘ 
deceleration in g's 
BCV,40 6021-6060 Array of lower corridor values for 
velocity in ft/sec 
BCWLI1, 40 6501-6540 . Array of lower corridor values fo wake 
length in meters for first radar frequency 
BCWL2, 40 6581-6620 Array of lower corridor values for wake 
length in meters for second radar frequency 
BCWL3, 40 6661-6700 _ ‘Array of lower corridor values for wae 
length in meters for third radar frequency 
BCWRI, 40 6261 -6300 Array of lower corridor values for wake 
cross séction for first radar frequency, 
see IDBL for unites 
BCWRz2, 40 6341-6380 Array of lower corridor values for wake 
; cross eection for second radar frequency, 
see IDBL for units 4 
1-468 
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Output (cont'd) 

" Common 
Name Location 
BCWR3, 40 6421-6460 
DB, 160 2081-2240 
DD, 160 1921-2080 
DV, 160 1761-1920 
DVAR, 50 11707-11756 
DWL1, 160 2721-2880 

” 
DWL2, 160 2881-3040 

7” 
DWL3, 160 3041-3200 
DWRI1, 160 2241-2400 
DWR2,160 2401-2560 
DWR3, 160 2561-2720 
HOLD, 150 7101-7250 

' : SAVE, 4320 7251-11570 

T469 
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Description 


Array of lower corridor values for wake 
cross section for third radar frequency 


The array of difference between the R/V 
and decoy ballistic coefficient 


The array of difference betwern the R/V 
end decoy deceleration in g's 


The array of difference between the R/V 
and decoy velocity in ft/sec 


Stored values of certain elements of the 
@CCUR array 


The array of difference between the R/V 
and decoy wake length in meters for first 
cadar frequency 


The array of difference between the R/V 
and decoy wake length in meters for second 
radar frequency 


The array of difference between the R/V 
and decoy wake length in meters for third 
radar frequency 


The array of differeuce between the R/V 
and decoy wake cross section for first 
radar frequency, see IDBL for units 


The array of difference between the R/V 
and decoy wake cross section for second 
radar frequency, see IDBL for units 


The array of difference between the R/V 
and decoy wake cross section for third 
radar frequency, see IDBL for units 


A 3 by 50 array of stored design 
parameters 


A9 by 160 by 3array of stored differences 
in design parameters for determining | 
influence coefficienta. | 
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Common * 
, Name Location Description 
SB, 40 6821-6860 See input 
SD, 40 6781-6820 See input re 
SV, 40 6741-6780 See input 
SWLI1, 40 6981-7020 See input 
SWL2, 40 : 7021-7060 See input 
SWL3, 40 7061-7100 See input 
SWRI, 40 6861-6900 See input 
SWR2, 40 6901-6940 See input 
SWR3, 40 6941-6980 See input « 
TCB, 40 6221-6260 See input 
TCD, 40 6141-6180 See input , 
TCV, 40 6061-6100 See input 
TCWLI, 40 6541-6580 See input 
TCWL2, 40 6621-6740 See input 
TCWL3, 40 6701-6740 See input 
TCWR1, 40 6301-6340 See input 
TCWRz2, 40 6381-6420 See input 
TCWR3, 40 6461-6500 See input 
TT» 200 5561 = 5760 Combined altitude array 
TTM, 200 TTMCOM Combined time array 
TS, 160 3201-3360 Array of saved R/V trajectory times _ 
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Description 
Array of R/V or decoy trajectory times 
Array of stored wake length values in 


meters at first radar frequency for 


Array of stored wake length values in 
meters at second radar frequency for 


Array of stored wake length values in 
meters at third radar frequency for 


Array of stored values for wake cross 
section in units specified by IDBL at 

first radar frequency for reentry vehicle 
Array of stored values for wake croas 
section in units specified by IDBL at 
second radar frequency for reentry vehicle 
Array of stored values for wake cross 
section in units specified by IDBL at 

third radar frequency for reentry vehicle 


Array of altitudes for R/V or decoy 


Array of saved altitudes for R/V trajectory 


Second subscript of the HOLD array 


Second subscript of SAVE array 


Number of altitudes in the combined 
altitude array, TT 


Trajectory option code 


Location in the OCCUR array of design 


1-47] parameter being varied 


1000} Ltn 4 ANNO HAA ee Og 


W 
Output (cont'd) 
Name Location 
TPLOT, 160 1-160 
WLIS, 160 4481-4640 
reentry vehicle 
WL2S, 160 4641-4800 
reentry vehicle 
WL3S, 160 4801-4960 
reentry vehicle 
WRIS, 160 4001-4160 
” 
WR2S, 160 4161-4320 
‘ 
WR3S, 160 4321-4480 
ZPLOT, 160 161-320 
trajectory 
ZS, 160 3361-3520 
*eeKNDOTE** The following quantities are integers: 
II IDCCUR( 308) 
INALPH NOCCUR(30) Option to Input o< 
JJ IDCCUR(310) 
K IDCCUR(31 3) 
é LOPT NOCCUR(07) 
NCON IDCCUR(312) 
» 
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NDDCOY 


NGEOM 


NLOW 


NTHRST 
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Output (cont'd) 
en NaN CAEL CELL CAAA 


Common 
Location 


a eS re a ae ea el RE a ri 
+HXNOTE** The following quantities are integers (cont'd) 


INFTST 
N@CCUR(L5) 


NOCCUR(22) 


ORL LLCO CODE ECE CO, ee vee too werent et TE Ae AP ee 


Description 


Control code used in INTGRL 
Input geometry option code 
A control code used in calling HEADER 


Input thrusting option code 
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* 1 
' 4, Numerical Methods ia 
The foliowing abbreviations will be used: te 
rvt means reentry vehicle trajectory. | 
\ 
bdt means basic decoy trajectory. 
cdt means comparison decoy trajectory. ' 
F123 saves three sets of tables for use by other subroutines, In the lorad { 
set of tables the values of up to nine characteristics of rvt are saved for each ; 
altitude at which SUBR@UTINE VIXEN outputs data, To obtain the second set | 
of tables, the values of the corresponding chavacteristics for a bdt or a cdt is : | 
subtracted from that of the rvt for each characteristic at each altitude. The i 
Z third table saves up to three sets of these differences in the SAVE array, 
In the table below, the first column contains the name of the array in which the 
Ld 
values of the characteristic for the rvt are saved, The second column contains 
the name of the array in which the differences are saved, The third column 
contains the name of the characteristic whose differences are being saved, 
rvt Values Differences Characteristic Name 
Vs DV velocity 
DS DD deceleration 
BS DB ballistic coefficient 
WLIS DWL1 wake length at radar frequency one 
WL2sS DWL2 wake length at radar frequency two 
WL3S DWL3 wake length at radar frequency three 
WRI1S DWR1 cross sectionat radar frequency one 
WR2S DWR2 cross sectionat radarfrequency two 
WR3S DWR3 cross sectionat radar frequency three 
td 
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\ The two statements:- 100 CALL TIMERS(1, TIMTAB) 
CALL TIMERS(-1, TIMTAB) 


is | provide for the possibility of printing out the elapsed time between execution 


of the two statements, In the present deck this subroutine is a dummy 


producing no output. 


After initializing NDC and NDVC, F123 calls SUBRQUTINE SAVEDV to save 

the values of SCCUR(133-145), BCCUR(205-209) and PCCUR(222) in DVAR(1-19) 
respectively, These values will be restored to the CCUR array later, Next 
IREF is tested, If IREF is not equal to two go to 497, The transfer to 497 


means that the rvt is to be input or calculated when IREF is three or one 


respectively, If DREF equals two, NDECOY is tested, If NDECOY equals 


one go to 497. In this case, either a bdt or the rvt is to be processed, If 


NDECOY equals two, one cdt will be calculated for each design parameter being 


perturbed, If NDECOY equals three, two cdt's will be calculated for each design 


parameter being perturbed, If NDECPY is either two or three, the D® loop» 


nd ing at 111 is entered to save the values of the design parameters being per- 


turbed, These are saved inthe HOLD array, HOLD(1,i) saves the value for 


b, the basic decoy. H@OLD(2,i) saves the value for the first comparison decoy 


and HOLDX3, ij saves the value for the second comparison decoy, NDVC is 


the number of design parameters to be perturbed in one case, The input arrays 


DVL and DVH contain the perturbed values for the first and second comparison 


decoys respectively, 


At 497, the DO loop ending at 777 is entered and IREF is tested, If IREF # 2 


contro: passes €0 498. Control passes to 401, 402 or 405 when NDECOY * 
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is one, two and three respectively. At each of these three statements the values 


of NL@W and NHIGH are set and statement 404 is executed next. At 404, 
NCQ@N is defined. This is the location in the OCCUR array of the Ith parameter 


being perturbed. Next, the call to SAVEDV restores the values saved in DVAR 


to the OCCUR array. Then going to 499, NL@W and NHIGH are set for the rvt 


and 499 follows. 


‘ 


At 499, the D@ loop ending at 444 is entered to obtain and process trajectory 
data. If KK equals two, OCCUR(NCQN) is set equal to DVL(II) for the first 
comparison decoy. If KK equals three, OCCUR(NC@N) is set equal to DVH(LI) 
for the second comparison decoy. Next, if M@DE is not equal to one, the case 
number is increased by .001 to aid in identifying printed output and plots. 

Next, the case, date and memo number are printed and IREF is tested again, 
and if = 3 passes to 28, The sequence following this test and ending at 56 
tests the input options and may reset some of them. The next three set values 
for numerical constants. The teticwing three set input options equal to storage 
quantities in order to carry them over to the following input case unless they 
are reset, The next eight statements change units from inches (input units) 

to feet (units Sige internally), The next ten statements set drag coefficient 
coefficient quantities to zero. Then L@PT is tested. If LOPT equals zero go 
to 201, If L@PT is not equal to zero, the next nine statements set rotational 
quantities equal to zero, If L@PT = 0, the quantities related to maximums and 
minimums in ef are zeroed. The next four find the trigonometric functions of 
the thrust offset angles. The next six change angles from degrees to radiana, 


In the next thirteen, the masse loss and heating quantities are set equal to zero, 
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i The call to TIMERS next prepares to measure the time taken to calculate 


trajectory, Then VIXEN is called to calculate the trajectory and TIMERS is * |i 


| i called twice, the first time to print the time elapsing for the trajectory 


calculation and the second for the time elapsing between the execution of i! 


td statement 100 and this call to TIMERS, The following six statements set the 


thrust quantities equal to zero, The next eight reconvert geometric quantities 


from feet to inches, Next the value of ALST is changed from degrees to | i 


radians and the next three statements set options equal to stored input 


quantities to prepare for the next case, 


Now IREF ie tested, If IREF equals two, go to 30, If IREF equals one, the | 


data for the rvt is saved in the DD loop ending at 22. Then LP is checked, 4 | 


If LP < 6 control passes to407, If LP equals six, the rvt calculations have 


failed, A méssage to this effect ia printed, SAVEDV is called to restore |" 


values from DVAR to the OCCUR array, and control returns to the main program. 


At 30, HEADER is called to indicate the type of decoy isajectory being cal- 


culated and print selected quantities associated with it, Then, in the DO 


loop ending at 33, the differences described in the table above are saved, Next 


LP is tested as above; if equal to 6, a message is written, SAVEDV is called, 


and the program goes next to 777, If LP is less than 6, control passes to 507, 


In the sequence of statements beginning at 507 and ending at 50, the altitudes 


in the input array H and the altitudes at which VIXEN produces output are 


combined into one array, monotonically decreasing, If the same altitude is 


found in both arrays, it is counted only once, The altitude values are put in : 
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the TT array and the corresponding times are saved in the TTM array, The 
number of elements in each of these arrays is K, Next in the DO loop ending 
at 333, the differences needed by the influence coefficients SUBROUTINE 
INFCOF are saved. The nine statements between 333 and 444 provide for the 
calling of SUBROUTINE INTGRL which may calculate corridor integrals and 
sigma square integrals as well as call the plotting routine, AVPLT, and the 
influence coefficients SUBROUTINE INFCOF, After 777 SAVEDV is called 


to restore values to the @CCUR array and control returns to the main program, 


5, Other Information 
A, F123 is called by the main program. 
B, F123 calls 
(1) SUBR@UTINE HEADER 
(2) SUBROUTINE INTGRL 
(3) SUBROUTINE SAVEDV 
(4) SUBROUTINE VIXEN 


(5) The dummied SUBROUTINE TIMERS 


C. F123 calls the IBM supplied finctiona SIN and COS, 


fr 
i 
t 
1-477 
* 
alti Pre ree) en 11) ea HAE YwPeAN 1 P anit "eg nate oar er otawe Ey HS mB 
Aue d eee as ha \) IL Ah tn ee 


i i {| il 
TUM ATT biti ifs 
il | | wai te 


of MMVI Tee 11 crvEReAPHA eto fnertHAe OEOGON HEYA Hipo9t yemings oe ay saete ves avn ‘secturbe} tote YT ET TED MEYER ABARAT, AT Tad tocar cme omeno ame TAP OT ag 06 
mi 
h 
mi 
et 
Wl 
i 
lh 
1 I 
: th 
boa 
owl 
fit 
| qh 
| lft 
il 
Ih 
iH 
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SUBROUTINE INTGRL (KA, KB 


l, Purpos e 


INTGRL is intended to perform one of four main tasks depending 


upon certain input options. 


A. To evaluate a corridor integral and print a tabular array of 


data at each of a number of altitudes. 


B. To evaluate a sigma square integral. 


C. Call the plot routine to obtain a corridor plot. 
: D, ‘Call the plot routine to obtain an influence coefficient plot. 
ih 
* 
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| 2. Input ! 
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Name Source Common Block Descript 
| ee ee eee 
AB? SR2490 or ZREADX PCCUR (11571-11597) See text 
B, 40 LINFIT PCCUR (4961-5000) See text 
C, 40 LINFIT PCCUR (5001-5040) See text 
D, 160 LINFIT PCCUR (5041-5200) See text 
E, 160 LINFIT PCCUR (5201-5360) See text 
F, 40 LINFIT PCCUR (5361-5400) See text 
G, 40 LINFIT PCCUR (5401-5440) See text 
H, 40 SR2490 or ZREADX PCCUR (5841-5880) Altituded 
IOP, 90 ZREADX I@CCUR (1-90) Option ¢ 
KA F123 Used to 
array ni 
KB F123 Used to 
array nj 
KD F123 See text 
KE F123 See text 
KX F123 IDCCUR (313) The nun 
altitude 
LPLOT VIXEN I@CCUR (302) The nur 
ae 
LX F123 Flag sp 
array al 
integral 
NCP SR2490 or ZREADX IOCCUR (304) Number! 
NDECOY F123 INFTST Test pat 
@CCUR, 4000 Unlabeled Common OCCUR (1-4000) 
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Description 


See text 


See text 

See text 

See text 

See text 

See text 

See text 

Altitudes at which corridor limits are input 
Option code 


Used to determine which elements of the IOP 
array need to be tested 


Used to determine which elements of the IOP 
array need to be tested 


See text 
See text 


The m.nber of elements in the T array of 
altitudes 


The number of elements in the Z array of 
altitudes 


Flag specifying which three elements of the A 
array are to be used in evaluating the corridor 
integral 


Number of elements in the H array 


Test parameter for influence coefficient plots 
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Preset 
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A(T) = 1 
If T=1 (MOD3) 
otherwise A(I)=0 
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2, Input (Concl'd) 
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Name Source Common Block 

Q, 40 SR2490 or ZREADX 

R, 160 $R2490 or F123 

S, 40 SR2490 or ZREADX 

SIG, 40 SR 2490 

T, 200 SR2490 or ZREADX PCCUR (5561 - 576U) 

TPLOT, 160 VIXEN PCCUR (1-160) 

Z, 160 VIXEN PCCUR (161-320) orate 
Array g 
data | 

we ae aie afc 

NOTE 


DC is a special input of alphanumeric data 
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Description Preset 
aon ce en ener etanennaneetainasnaentatmt anen en 

Array of lower corridor limits 4) 
See Text 0 
Array of upper corridor limits a 
Standard deviation in error of measurement 

pol - 5760) See Text 0 

r 160) Array of times corresponding to the Z array 

b 1-320) of altitudes 
Array of altitudes for which VIXEN stores trajectory 
data 
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Hl 
| i 3, Output 
| 4 Name Common Block Description Preset 
t | a A, 27 See Input See ZREADX for AA See input 
B, 40 See Input See LINFIT description 
Cc, 40 See Input See LINFIT description 
D, 160 See Input See LINFIT description 
E, 160 See Input See LINFIT description 
F,40 See Input See LINFIT description 
G, 40 See Input See LINFIT description 
HDIF PCCUR (11758) H(1) - H(NCP) i) 
OCCUR, 4000 = unlabelled Some values are changed 0 
common for use by the optimizer 
Ql Lower corridor limit at 
initial altitude 
Q2 Lower corridor limit at 
any other altitude. 
Rl Value of R at initial 
altitude 
R2 Value of R at any other 
. altitude 
\ §1 Upper corridor limit at 
initial altitude 
52 Upper corridor limit at 
’ any other altitude 
. T, 200 See Input See text 0 
x Tl Used by ADD, see Text 
T2 Used by ADD, see Text 
XX, 200 See Text 
+ ¥, 200 See Text 
ZZ, 200 See Text 
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4. Numerical Procedures 


Before the logical flow in INTGRL is discussed, it will be useful 


to describe in some detail what we call the corridor integral. The input 
array H contains a set of altitudes at each of which an upper corridor limit 
and a lower corridor limit has been input for one or more of nine different 
arrays, The appropriate set of upper and lower corridor limits are sent 


to INTGRL through the list in the arrays S and R respectively. The 


number of elements in each of these arrays is NCP, also an input quantity. 


The array Z contains the altitudes at which VIXEN outputs trajectory data, 


The array R contains the difference between the value of some specified 


characteristics of the reentry vehicle and the value of the same character- 


istic of a decoy at each of these altitudes. The number of elements in the 


R and Z arrays is LPL@T, which is defined in VIXEN. 


* SUBROUTINE LINFIT is used to fit straight line segments between 
the pairs of points (H(I), Q(I) and H(I+1), Q(UI+1)) 121,2,..., NCP -1l, 
The coefficients of the fit are saved in the arrays B and C, Similarly, 


the points (H, S) are fitted and the coefficients saved in the arrays F and G, 


and the points (Z, R) are fitted and the coefficients saved in the arraya D 


and E, We can then express the values of Q, R, and S at any altitude 
ALT by 

Q = ALT*B(I) + C(I) 

R = ALT*D(I) + E(I) 

S = ALT*F(I) + G(1) 


The area between the segmented lines formed by Q and R is called 


the corridor, So long as the value of R is greater than or equal to the 


value of Q and less than or equal to the value of S, no integration is done, 
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be omitted. LINFIT is then called three times to save the fitting 


test parameters NX and ITEST, the accumulator for the corridor integral 


‘ptatetnent at 10, 30, or 20 will be executed to write the first line of the 


If R is less than Q on an interval between ALT! and ALT2 the function | 
F(Z) = A(LX)#2" + A(LX + 1)*Z + A(LX + 2) is integrated over the area 
defined by lines R, Q, ALT! and ALT2 with respect to Z and the absolute 


value of the result is added to DA, the running sum of the integral. If 
R is greater than S on some interval between the altitudes ALT and 
ALT2, the same function is integrated over the region defined by the 
lines S, R, ALT1 and ALT2, these two being specific values of Z, The 


function subprogram ADD is used to perform this integration, | 


The subroutine flow proceeds as follows: Immediately upon 
entering INTGRL, I®P(KA + 2*KB) is tested. If IOP(KA+ 2*KB) equals 
0, no corridor integral calculations are desired and statement 1120 is 
executed next, If it equals 1, KQ is defined and tested in the next two * 
Statements. IF KQ= 1, the heading for the tabular output associated 


with the corridor integral is printed. The next two statements define 


IG@, IF JG@ is 0, the slope print out for the corridor integral will » | | 
coefficients described above. The next four statements initialize the 


DA, and the altitude at which R first leaves the corridor DB. Ql, Rl, il 
and Sl, calculated next, are the respective values of Q, R, and S at the ij 
initial altitude Z(1), These values are then saved as XX(1), YY(1), and | 
22(!) for the corridor plots. Next the value of NX is determined. NX 
will be 0,1, or 2 depending on whether R1 is less than Q1, between Ql 
and $1, or above Sl], If Ql is already outside the corridor, the next two 
statements will change ITEST to 1 and set DB = Z(1), Then if KQ = 0, 
statement 40 will be executed next. If KQ= 1, NX is tested and the WRITE 


corridor integral table as NX is 0, 1, or 2 respectively, Then statement 40 
ia executed, 
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t Following the statement 40, NS and NR are defined and the 
} D@ loop ending at 1111 is entered to calculate the corridor integral. 
| 


se 


The array T, used here, was generated in F123 by combining 


the Hand Zarrays. If an element is common to both arrays, it is 


included in the T array cnly once. The T, H, and Z arrays are 


ee 


monotonically decreasing. It will usually happen that for some J, T(M) 


will lie between H(J) and H(J+1). If so, the value of J is found and 


gaved in the D@ loop ending at 44. If not, J is set equal to NR. 


Similarly, it will usually happen that T(M) lies between Z2({1) and Z(1+ 1) 


for some I. If so, the value of I is found and saved in the D@ loop ending 


at 55. Otherwise Lis set = NS. 


The six statements starting at 60 set Q2, R2, and S2 equal to 
the values of Q, R, and S respectively at the altitude T(M) and then save 
« these values in XX(M), YY(M) and ZZ(M) respectively. The next three 
statements find the value of NY which will be 1, 4, or 7 accordingly as 
R2 is less than Q2, between Q2 and S2, or greater than SZ, NX and NY 


now describe the behavior of the R line between the altitudes T(M-1) and 
T(M), IG@ is set = NX + NY andthe computed GO T@ is executed, sending 
control to the statement having the statement number = 100*IG@. This 
transfer will result in the execution of some or all of the following six 


types of statements. 


A. T2, if not = 0, is defined as the altitude at which R crosses 
either Sor. It will then be either the higher or the lower of the two 


altitudes used to define the region of integration. 


B, Tl = T(M-1). 
Tl, when used is the higher altitude defining a region of integration, 
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Cc, DA=DA+ ADD : 


The function subprogram ADD performs the integration and adds 1 | 


the result to DA, 


DBD, CALL WRITEM 
Wf KQ = 1, a line of output is printed in the corridor invegral table 
for the altitude T(M). 


E. IF ITEST =9, G@T®950 
If ITEST = 0, statement 950 is executed next. If not, the following 


atatement is executed. 


F. GOTO 999 


Statement 999 is executed next. 


A description of the region of integration at each of the statements 
100, 200, .,. 900 follows, The following notation is used in the description 
The general form is. at XXX LI, L2, L3, L4, L5, where XXX ia the 


statement number, L1 is the lower line of the region of integration which may 


be either a triangle or a trapezoid; it is always either Ror S, 12 is the 

| “pper line for the region; it is always Qor R. 13 is the higher of the two 
altitudes defining the region and L4 is the lower altitude. L3 and L4 are 
parallel. 5 may be any of the four RQB, RSA, RQI and RSI meaning 
respectively R crosses Q from below lower corridor, R crosses S from above 
upper corridor, R crosses Q from inside corri‘or, and R crosses S from 


inside corridor. If there is no crossing, L5 is omitted. 


At 100 the lines R, Q, Tl, T2 define the region of integration. 
At 200, R, Q, T2, T(M), RQI define the region of integration. 


At 300, two regions are needed since R goes from above S to below 


Q. These regions of integration are described respectively by the lines 
S, R, Ti, T2, RSA and by the lines R, Q, T2, T(M), RQL * 
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At 400, R, Q, Tl, T2, RQB describe the region of integration, 


At 500, no integration is required. R is between Q and S for 


the entire interval. 


At 600, S, R, Tl, T2, RSA describe the region of integration. 


At 700, two regions are required since R goes from below Q to 
above S, These regions of integration are described respectively by the 


lines R, Q, Tl, T2, RQB and by the lines S, R, T2, T(M), RSL 
At 800, S, R, T2, T(M), RSI define the region of integration, 
At 900, S, R, Tl, T(M), RSI define the region of integration. 


At statement 950, ITEST is changed to 1, then DB is set = T2 and 


statement 999 is executed. 


At 999, NX is redefined to specify the relation of the beginning of 


the next R line segment to the corridor. 


After statement 1111 has been executed KX-1 times, DA and DB 
are put into the OCCUR array for possible use by the optimizer. Next, 
if KQ = 0, the values of DA and DB are printed along with the name of the 


variable in the R array. 


Next, IMP(KA + KB) is tested, 


executed, 


If it equals 0, statement 2300 is 

We use 
the trapezoidal rule to integrate the function (R/s1G)" from H(1) to H(I+ 1) 
I= 1,2,...,NCP-1 with respect to H, 


If it equals 1, we compute the sigma square integral, 
The running sum is kept in the 
variable SUM. The values of SIG are input in the SRS array. As a preliminary 


RL®@ is set = cRiv/sic(iy)? and SUMis sets 0, The values of R at 
the altitudes H(2), H(3),...H(NCP) are found by linear interpolation in 
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the D@ loop ending at 2222. SUBROUTINE INTERP finds these values, 
RHI is defined as (rin /sicu)?. The value of the integral from H(I-1) 
to H(1) is computed and added to SUM. RL®@ is then set equal to RHI 
and the loop executed again. SUM is then put into the OCCUR array 
for possible use by the optiniizer. Next I®P(KA) is tested. If it equals 
0, statement 4000 is executed. If it equals 1, AVPLT is called to 


produce a corridor plot. 


At 4000, if NDECOY is less than 3, control returns to F123, 
If it equale 3,1@P (KA+5*KB) is tested. If it equals 1, AVPLT is called 
to produce an influence coefficient plot and control returnsto F123, If 


it is equals 0 control returns to F123 immediately. 


5, Other Information 


A. INTGRL is called by F123 only 


B. INTGRL calls 
1, SUBROUTINE LINFIT 
2, SUBROUTINE WRITEM © 
3, SUBROUTINE INTERP 
4, SUBROUTINE AVPLT 
5, FUNCTION ADD 
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FUNCTION ADD (L,R, S, A, B, C, D, E) 


l, Purpose 


$ 
Add evaluates of F(t)dt over a trapezoidal region in the t, y 
R 


plane defined as the area enclosed by the four lines: 


I t=R 


. t= 5 
y= A*t +B 
. y= c*t+D 


” 2 
{t) = + 
Fit) Et + EL 41 t E+ 


ee Ww & 


Name Source Desc ription 


INTGRL curve fit coefficient 


A 

B INTGRL curve fit coefficient 
Cc INTGRL curve fit coefficient 
D 
E 


INTGRL curve fit coefficient 


INTGRL AA, weighting coefficients for corridor 
integrals, 11571-11597 


L INTGRL integer indicating which AA's are to be used 
R INTGRL altitudes in feet 
i s INTGRL altitudes in feet 


| 
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i] The E array is input under the name AA, It reached SUBROUTINE 
INTGRL through the PCCUR common array and INTGRL then sends it 
to ADD through the argument list. 


are preset to 1,0 and the rest of the E array is preset to zero in subroutine 
SR2490, Rand S are altitudes in feet. 


Name Description 


ADD corridor integral increment 


4, Numerical Procedure 


Since the height of the trapezoid is(C-A)*t + (D-B), we may write » 


s 
ADD = “df F(t)*((C-A}*t + (D-B))*dt, 
ry ae 


a | r ’ meat tk ‘This can be integrated explicitly as a function of t to give the integral I, 


2 (C-A) (E, 


t 
4 
3 t 
(D-B) (E, + +E St +E t) 


The value of I at t= R is subtracted from its value at t = S and the result 
is placed in location ADD. 

_--s PUNC'TION ADD is called by SUBROUTINE INTEGRAL only, 
FUNCTION ADD does not call or reference any subprogram. 
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, SUBROUTINE INFCOF (DC, 12, ZQ) , 


1, Purpose 


SUBROUTINE INFCQ@F calculates a table of influence coefficients 


’ and possibly second derivatives, at each altitude for which VIXEN stores 


trajectory information 


2. Input 


*indicates integer quantity 
Name Common Block Source Description 
, einai ites 
DC, 2 - INTGRL name of delta 
H, 150 PCCUR(7101-7250) F123 holds design variable values 
Il IMC CUR(308)* F123 second subscript of H 
IZ, 19 - INTGRL locations, in OCCUR array, of 
19 special parameters 
JJ I@CCUR(310) * F123 first subscript of SAVE 
LPLOT I@CCUR(302) * VIXEN number of output altitudes 
NC@N I@CCUR(312) * F123 location in the @CCUR array 
ek cece parameter being 
NDEC OY I@CCUR(305) * Ret test parameter 
S, 4320 eet 5 , F123 array of saved deltas 
T, 160 PCCUR(1-160) VIXEN time array 
XQ, 40 - INTGRL array containing design 
parameter names 
Z, 160 PCCUR(161-320) VIXEN altitude array 
bal 
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3, Output 


SUBROUTINE INFCOF produces only printed output. 


4, Numerical Procedure 


The following definitions help to explain the meaning of the influence 
coefficients and the second derivative and the method of calculating them. 


The three deltas which may be printed out at each altitude are 


First delta equals S(JJ,1, ') 
Second delta equals S(JJ, I, 2) 
Third delta equals S(JJ, 1, 3) 


The array called S here is called SAVE in F123 which produces it. The 
first subscript, JJ, specifies the appropriate of the nine possible vehicle 
characteristics for which deltas have been saved, The second subscript, 
I, indicates the altitude at which the deltas were computed, The third 
subacript is 1, 2, or 3, when the decoy being compared with the reentry 
vehicle is respectively a basic decoy, a first comparison decoy and a 


second comparison decy. 
First influence coefficient equals 


S(JJ,1, 2) ~ S(JJ,1, 1) 
H(2, 0) H(1, I) 
Second influence coefficient equals 
S(JJ, 1, 3) - S(JJ, 1, 1) 
H(3, II) rr H( 4, I) 
The array called H here is called HOLD in 1123 where it is generated, 
The subscrip: I comes from F123 and specifies which design parameter is 


being perturbed, The first subscript of H has the same meaning as the 
third subscript of S above, 
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* Second derivative equals (second influence coefficient minus 


first influence coefficient) / (H(3, II) minus H(2, II)} 


The flow proceeds as follows: The D@ loop ending at 1? tests to 
determine if the design parameter being perturbed is one of the nineteen 
for which names have been stored in the XQ array. If so, JP specifies 


which of these names will be printed as part of the table identification. 


Otherwise, JP is set equal to zero and the value of the integer NCQN is 


printed instead. Next X21, the difference in H, is defined to be used in 


calculating the first influence coefficient. Then, if NDECQ@Y equals two, 
statement 100 is executed next. If NDECQ@Y equals 3, X22, the difference 


in H, is defined to be used for calculating the second influence coefficient and 
X31 is defined to be used in calculating the second derivative. 


- At 100 AMID is defined as the arithmetic average of H(1,II) and 


H(2, II) for print out. Some preliminary information is printed, then a 
table heading having the five columns: ALTITUDE, TIME, FIRST DELTA, 
; SECOND DELTA, INFLUENCE COEF is written, In the D® loop ending 


at 222 one influence coefficient is calculated at each altitude and a line of 


five numbers is printed out at each altitude. The control returns to INTGRL, 


At 500, the procedure is quite similar, The three arithmetic averages 


AMID, BMID and CMID are computed for preliminary print out. The table 
heading this time has eight columns: ALTITUDE, TIME, FIRST DELTA, 
SECOND DELTA, THIRD DELTA, FIRST INFLU COEF, SECOND INFLU 
COEF, SECOND DERIV. 


In the D@ loop ending at 555, two influence coefficients and a second 


derivative are computed at each altitude and the eight specified values are 


prin.ed under the appropriate headings. Then control returns to INTGRL, 
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5. Other Information i q 


A, INFCQOF is called by INITGRL only. 


B. INFC@F does not call or reference any other subprogram. 


F ' 
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SUBROUTINE INTERP (M,N, W, X,Y, Z) | 4 


1, Purpose | 4 


Given the two tables X and Y(X) and, in addition, a particular value 


of X, called W, SUBROUTINE INTERP finds Z equal to Y(W) by linear 


interpolation, 


2. Input 


Name Source Description 


M WAKE dimension of X and Y arrays in WAKE 
N WAKE size of X and Y tables 

WwW WAKE particular value of X 

xX WAKE table of independent variables 

Y WAKE table of dependent variables 


Name Description 


Zz value of Y(W) found by linear interpolation 


4, Numerical Procedure: 


First N is tested to determine if it is greater than zero and less 
than or equal to M. If itis, statement 20 is executed, If it is not, an 
error message is printed and SUBROUTINE EXIT is called to terminate 


the computer run. 


a 
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At statement 20 X(1) is compared with X(2) to determine whether 


the X array is increasing or decreasing. The X array must be either 


monotonically increasing or monotonically decreasing. If it is not, 


results will be unreliable. If X(1) is less than X(2) control is transferred 
to statement 50 and the value of I, such that W is between X(I-1) and 
X(I), is found. If W is outside the range of the X table, I will be either 


2 or N, accordingly as W is less than X(1) or greater than X(N) 
respectively, The value of Z is then computed at statement 40, If 
X(1) is greater than X(2), the value of I is found in the DO loop ending 


at 22 and Z is again calculated at statement 40, 


5, Other Information 


A, SUBROUTINE INTERP is called by SUBROUTINE WAKE, 


B, SUBROUTINE INTERP calls the IBM SUBROUTINE EXIT. 
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SUBROUTINE LINFIT (ND, NENTRY, X, Y, A, B) 
- 
1, Purpose 
Given the NENTRY points (X, Y), LINFIT finds the equation of the 
line passing through the points (X, Y), (X,, Y), in the form y = A, x + B., 
fori=1,2,..., (NENTRY-3). These coefficients can then be used 
by the calling program for rnaking a linear fit when the i:terval required 
is known in advance. 
2. Input 
Name Source Description 
ND F123 and/or INTGRL adjustable dimension 
NENTRY F123 and/or INTGRL number of elements in X, Y, A and B 
¥ 
X, ND F123 and/or INTGRL array of independent variables 
Y, ND F123 and/or INTGRL array of dependent variables 
w 
3, _Output 
Name Symbol Description 
A, ND A. coefficient of X in linear fit 
B, ND B. constant in linear fit 
4, Numerical Procedure 
A,FY,4, 7 yl 4, 7%) and B.*y,7 A, %, are the equations 
evaluated by SUBROUTINE LINFIT. 
5. Other Information 
A. LINFIT is called by SUBROUTINE INTGRL and by SUBROUTINE 
F123, 
‘ B. It does not call or reference any other subprogram. 
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SUBROUTINE WRITEM(X, Z, T, Q,R, 5S, DA, DB, SLOPE, N) 


4 1, Purpose 


Each time SUBROUTINE WRITEM is called, a single line may be 
printed in the corridor integral table. 


2, Input 


| The source of all input quantities is SUBROUTINE INTGRL 
| “indicates integer quantity 


Name Description 


DA cumulative value of the integral 


DB altitude at which delta enters or leaves corridor 


[x* test parameter for the write statement to be used 


N* control code; if N equals zero, print out is deleted 


Q lower corridor value at altitude Z 
R value of delta at altitude Z 
& upper corridor value at altitude Z 


SLOPE slope of delta on latest line segment 
T ‘ time at altitude Z 
z altitude 


3. Qutput 


Only printed output is produced as described in Section 4. 
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4, Numerical Methods i 


First N is tested. If N equals zero, no printing is done and 
control returns to INTGRL. Otherwise, IX is divided by ten so that it 


can be used as the index of a computed go to statement. 


Control is then transferred to 10, 20, 30, 40, 50, or 60 when 
IX is 1, 2, 3, 4, 5, or 6 respectively. At each of the six write 


statements at 10, 20, 30, 40, 50, and 60, the first six columns of the 


integral table are filled. These contain the values of altitude, time, 
lower corridor, delta, upper corridor and integral value. At 20, 


the leave corridor column is also filled, At 30, the enter corridor 


m column is filled. At 40, the slope column is filled. At 50, both the 
leave corridor and slope columns are filled. At 60, both the enter h 
corridor and slope columns are filled. After the line has been printed, | ( j 


. control returns to INTGRL. 


5, Other Information 


A. SUBROUTINE WRITEM is called by SUBROUTINE INTGRL only. | , 


B. SUBROUTINE WRITE calls no other subprograms, 
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4 SUBROUTINE SAVEDV(K, DVAR, OCCUR) 


1, Purpose 


SAVEDV stores the values of @CCUR(133-145), ®CCUR( 205-209) 
and @CCUR(222) in DVAR(1-19) respectively and later restores these 
values to the OCCUR array. 


2, Input 
*indicates integer quantity 
Name Source of Input Descrip<ion 
’ 
DVAR, 50 F123 array to save selected OCCUR values 
K * F123 test parameter 
OCCUR, 4000 F123 - common array 
* 
3. Output 
Name Description 
DVAR, 50 See Input 
@CCUR, 4000 ; See Input 
4, Numerical Procedures 
If the input test parameter K is one, statement 100 is executed. If 
K is two, statement 200 is executed. 
A 100 the values of @CCUR(133-145), @CCUR(205-209) and @CCUR(222) 
are saved as DVAR(1~-19) respectively and control returns to F123, 
LS 
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At statement 200 the values stored in the DVAR array at 


statement 100 are restored to the appropriate locations in the OCCUR 


array and control returns to F123, | ! 


5, Other Information 


A, SAVEDV is called by F123 only. 


B. SAVEDV does not call or reference any other subprograms, | 
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4,3 Miscellaneous Printout i 
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SUBROUTINE HEADER(K) 


l, Purpose 


HEADER prints the type of vehicle whose trajectory is being 


processed along with the values of nineteen design parameters, the 


value of NGE@M, an input option,and the value of LP an error test 


parameter. 
2, Input 


penta 


erent neers snetennnnerret rrarent-snennnernmerneca na in 
Name Source Common Block Description 
ee tepennnntennmetntnareneteten 


K F123 test parameter : 
N@CCUR, 30 READIT N@OCCUR input options 
OCCUR, 4000 @CCUR common array 
LL tN ne Retr ent eteratentant tetas ninteteettesatteenneeereeenaeteetieene mene 
3. Output 

SUBROUTINE HEADER produces printed output only, 
4, Numerical Procedure 

Immediately upon entering HEADER, K is tested. The statements 
10, 20 and 30 are executed next when K is 1, 2, and 3 respectively, 

At 10 the heading "BASIC DECOY TRAJECTORY" is printed and 
control is transferred to 40, $ 

At 20 the heading "COMPARISON DECOY CHARACTERISTICS" 
is printed and control is transferred + fo 

I~ 
* 
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At 30 the heading "REFERENCE REENTRY VEHICLE CHARACTER- 
ISTICS" is printed and statement 40 follows. 


At 40 the values of nineteen design parameters, the option NGE@M 


H | and the test parameter LP are printed with appropriate headings. Then 


control returns to F123, 


5. Other Information 


A, HEADER is called by F123 only, 


B, HEADER does not call or reference any subprograms, 


‘ | lt 
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SUBROUTINE TIMERS 


1, Purpose 


SUBROUTINE TIMERS is a duminy subroutine, 


None 


4, Numerical Procedure 


The original TIMERS has the facility to print the time elapsing 
between the execution of any two specified statements in a subroutine, 
It has been replaced by the dummy because this print out is not desired 


at present. 


5, Other Information 


A, SUBROUTINE TIMERS is called by F123 only. 


B, SUBROUTINE TIMES does not call or reference any other 


subprogram. 
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5,0 EFFECTIVENESS OPERATIONS 
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Effectiveness Operations J 


The implementation of the decoy effectiveness model 


technology is contained in SUBROUTINE EFFECT, 


a 
4 
* 
id A 
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SUBROUTINE EFFECT 


1, Purvnose 


SUBROUTINE EFFECT calculates quantity PD which is the 


probability that a decoy will be discriminated. SUBROUTINE EFFECT 
receives from input a quantity PFD (the probability of false dismissal 
of a reentry vehicle and calculates a quantity SIGMA. Next T is 
found such that 

2 

2 


t x 
4, l — 
PFD = om e dx 
-o 


‘ This is done by a Newton-Raphson iteration scheme. Finally PD is 


calculated as 


T+SIGMA | 2 
Se 


PD * in e dx 


- 
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2. Input 


a 


a ERSERRE EERE ERCROnNE 


Common 
Name Source Block Preset Description 


ccsenameer ee 


mel 


LMP, 90 SR2490, ZREADX IOCCUR (1-90) 1 option parameter 


HDIF INTGRL PCCUR (11758) 0 total length of corridor 
PFD §R2490, ZREADX PCCUR (11757) 0 probability of faise 
dismissal 
SRS, 9 SR2490, ZREADX PCCUR (11598 - 0 number of smooth radar 
11606) samples 
3.  Qutput 
| 
Name Common Block Description 


PD @CCUR (3962) probability of discrimination 


SIGMA @CCUR (3965) The sum of nine weighted squares 
divided by HDIF, SIGMA equals the 
square root of the quotient 
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4, Numerical Procedure 


The D@ loop ending at ll accumulates the sum of nine weighted 
squares, Then the two statements following 11 divide this sum by HDIF 
and set SIGMA equal to the square root of the quotient. TOLD is set 
equal to -2 as the initial guess for T. Y is the corresponding value of the 
integral at this value of T. Kis an iteration counter. The IBM function 


subprogram 


1 r 
2 -x 
DERF(T) “Sit | e dx 
0 


determines the error function used in the integration. If we substitute 


XA@ tor X then the desired integral can be expressed as 


Wada ss z 
INT = 5 + (1 + DERF (7m) ) 


The iteration starts at statement 10 where K is increased by lL. 
If K reaches 20 the program goes to99 where PD is set equal to one, an 


error message is written and control returns to FEV. 


If convergence is reached, the program goes to 99 where the value 


of the integral of (T + SIGMA) is calculated and control returns to FEV, 


If K is less than 20 and convergence has not been achieved, TOLD 


and Y are reset and the prograrn goes to 10 to increase K and try again, 


5, Other Information 
A. SUBROUTINE EFFECT is called by SUBROUTINE FEV. 
B. SUBROUTINE EFFECT calls the IBM routines DABS and DERF. 
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CLASSIC FUNCTIONS FOR TESTING OPTIMIZATION 


TECHNIQUES 


6.0 


t) Classic Functions 


The following section contains the description of SUB- 


ROUTINE CLASSC, which contains classic functions to 


be used in testing the optimization procedures. 
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1, Purpose. 


SUBROUTINE CLASSC contains classic functions which are designed 


to provide difficult tests of how well a minimization routine works, 


2. Input 


*indicates integer quantity 


Common Source of 
Name Location Input Description 
A, 514 @CCUR(3ul - 814) READIT curve fit coefficients which 
are read in ¢ 
I¢P NOCCUR(3)* READIT code which determines the 
tl function to be minimized; 
input as MOPT 
X, 4000 @CCUR (1-4000) FEV the OCCUR locations defined 
in FEV 
3. Output 
LP @CCUR(4000) * error control code 
th X, 4000 @CCUR(1 - 4000) the output @CCUR locations 
| correspond to the function 
} value 
i “| rs 
Hii 1-514 
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4, Numerical Procedure 


salad 


SUBROUTINE CLASSC begins by testing the control code IOP. If 
I@P equals 3, 4, or 5 control passes to statements 300,401, or 501 respectively. 


If IP has none of these values it evaluates the equation 


2* 2 
F = 100, (X,-% ) + (lL ~ X)) +A 


1 ] 


then sets X00 equal to F. If the absolute value of IMP is 1, control then 


passes to statement 999, If not, the following equations are evaluated before 


control passes to 999, 
B= A), +(X,-A,){ Ay + (X, - Ad) [ 4, + (X)-A,) Aglt 


ial has “Ai 2) Aye t (XA IL Ayg + -A,)) eh 


Xi. 7 Az(X, - B) 


Koz * Agl*, - ©) 


Statement 300 begins the calculations for another function with the 
zeroing of the quantities FA, FB, FC, and FD, The quantities FA, FA2, 
FA3, and FA4 are then evaluated from the equations below and their values 


assigned respectively to X x X02 and X 


200’ ~°201’ 203° 


FA=x?+x? +2,x? 


2 
1 2 3 1%, - 5X; - 5% 


- 21,.X,+7.% +A 


3 4 \7 


2 2 2 


2 
FA2 = 8, + X,°X, +X, - Xi: x, - x, 
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. | Control then passes to statement 999, 


Statement 401 evaluates using the equation below the fourth function 


whose value is then assigned to X00" before control passes to 999, 


Fe4.X) +X, + Ay, (1./X) + 1/X,) a Ax 


The fifth function is evaluated by statement 501 and its value assigned 
to X00 before statement 999 is reached, 
Fz A te #aS X_ +A Ro HK X +A XU +A,. 
eee | ah ook ee. 3 2 41 aay! 6 
Statement 999 sets the dummy variable F321 to zero, then LP is A 
set equal to 1 before the return to FEV is executed. 


5. Qther Information | 


A. SUBROUTINE CLASSC is called by SUBROUTINE FEV. 


B. SUBROUTINE CLASSC calls on the IBM supplied routines 
IABS and FDXPI (exponentiation), 
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1.9 PLOTTER INTERFACE SUBROUTINES 


1~317 
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Plotting 


The subroutines utilizing in plotting, AVPLT, MAXMIN, and 


PLT are described in the followi. g sections. 


1-518 | 


_neeresnqepeneerprers = fe eee eee ce erteree wee + ----= tees * 
es nae cme 


N 


 mnrernenrnrrnremennerenemnenee 


a al ecle: \ AAMT 


SUBROUTINE AVPLT (IZ, XQ, X, Y, Z, KY, IC) 


1, Purpose 

AVPLT can produce any of four kinds of plots, depending on the 
value of KY. if KY = 1, influence coefficient plots will be prepared at 
a specified number of altitudes, NPA, NPA may be a maximum of 160. 
Three points are plotted on each graph. These points are not connected 
by straight lines. If KY = 2, one corridor plot is produced. Three 
curves are plotted on the same graph. The points of each curve are 
connected by straight lines. If KY=3, one plot of Y vs. X is produced. 
The points are connected by straight lines. If KY = 4, one plot of log, 0 
Y¥ vs. X will be produced. The points are connected by straight lines. 
If KY = 5, the plot file will be closed and the computer run will be 


terminated, 


T-519 


i 
| 


ia Input 


“indicates integer quantity 
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Name 


HOLD, 150 
Ic * 


I * 
IMPLOT * 
I@P, 90 * 
IZ, 19 x 
JI * 

KX * 
KY * 
LPLOT * 


NCON * 


NDECOY * 
NPA* 


NPV, 160 * 


SAVE, 4320 
T, 200 


XQ, 40 


ZPLOT, 160 


at ls Toh RP weft i lh 


Source of Input 


F123 
Calling Program 


F123 

SR2490 

SR2490 or ZREADX 
Calling Program 
F123 

F123 

Calling Program 
VIXEN 


F123 


§R2490 or ZREADXK 
SR2490 or ZREADX 


SP.2490 or ZREADX 


F123 
F123 
Calling Program 


Calling Program 


Calling Program 


Calling Program 
VIXEN 


Salsas holisslauauiibiistcjatsaianesielar aannanAbaliniaiHh 


Common Block 


PCCUR (7101-7250) 


IOCCUR (308) 
IOCCUR (309) 
IOCCUR (1-90) 


IOCCUR (310) 
IOCCUR (313) 


IOCCUR (302) 
IOCCUR (312) 


IOCCUR (305) 
IOCCUR (307) 


IOCCUR (141-300) 


PCCUR (7251-11570) 
PCCUR (5561-5760) 


ve 


PCCUR (161-320) 


hpi HAA A itt ep el RR 


ock Preset Value Description 
17250) i) design variable being perturbed 
- If KY = 3 or 4, test parameter for log plots. 
| If KY = 1 or 2, contains the vertical uxis title. 
) - Second subscript in HOLD matrix 
) 0 Test parameter for opening plot file 
0) 1 Option parameter 
- Array containing plot title information 
) - First subscript of the variable SAVE 
) - Number of points for each curve on corridor plot 
- Test parameter for kind of plot 
) - Number of altitudes at which information is output 
from VIXEN, 
| - Location in PCCUR array of design parameter 
being perturbed. 
] Test parameter for corridor plot title 
1 The number of altitudes at which influence 
| coefficient plots are produced, 
|-300) 1 Array containiny indices of altitudes at which 
influence coeffic.ents flows will be produced. 
~1 1570) 0 Dependent variable for influence coefficient plots. 
|-5760) 0 Expanded array of altitudes 
| - Dependent variable for all plots except influence 
coefficient plots 
- Array containing horizontal axis 
- Either second dependent variable for corridor plots or 
dependent variable for third or fourth type of plots 
- Third dependent variable for corridor plots 
320) 0 Array of altitudes at which VIXEN outputs information 


vith {eye FR I PHBE AEA YAH a HET eB PAN i HH i i ele 
: Lidia di pm | 


i iW 
' 


| ‘i " | 
iM i 4 Mee io eet pe ed 
1] 
qi Wi 
é : : 


an “” oe Z bh 
Ah abe bs a Geet atten iaci' ran nite 
wa 


yf -sponlnmmeteeacatorne ane acing FE MOEA EAE REO. 


4 


BLANK PAGE 


Oe 


7 


sil i Oeil MG ui : 


1 


ult Ce : ie 


it 


3, Output 
Name Common Block Preset Value . Description 
IMPLOT* IOCCUR (309) 0 Test parameter for opening 
plot file 
4, Numerical Methods 
Immediately upon entering AVPLT, the parameter IMPLOT is 
tested. If IMPLOT= 0, IDFRMV is called to open the plot file. IMPLOT 


is then changed to 1 and statement 10 is executed. If IMPLOT = l, 


statement 10 is executed next. 
At statement 10, KY is tested. If KY =1 or 2, statement 29 is 
executed, If KY = 3 or 4, control passes to statement 110. If KY = 5, 
* control passes to statement 230 where the plot file will be closed and 
the computer run terminated. After statement 20, the vertical axis 


plot title is defined. 


Next, in the D® loop ending at statement 33, the program starts 
to compare each element of the 1Z array with NCON. ‘If any IZ = NCON, 


statement 40 is executed next, the required subscript of IZ having been 
NCON, JP is set to zero, 


saved as JP., If none of the 12's and state- 


ment 45 is executed, 

At statement 40, JS and Jr are defined. These specify which 
two elements of the XQ array will be used for the horizontal axis title 
of an influence coefficient plot or of a corridor plot if a design variable 


is being perturbed, Statement 45 is executed next, 


1 aera 
Hey. i; ah 


T521 


HT caalsva/aldvaniwi alii fnisevsandan evil ais aaa gaal aan easinanneaennn iva 


ih 


Di inl a WAN GMRAUSGALUNTIILATINVUTwvanLIEiTUbso tas UUlcassbesy aul Ki UI VIVAL TRL LHVDADEDBL ULI eGLLUITVUVVHGGBIGd Lad and 


At statement 45, KY is tested. If KY = 2, statement 70 is 
executed next. if KY = 1, the three statements following statement 
45 define the W array of independent variables wht 1 then remain constant 


for each of the NPA influence coefficient plots to be set up in the D® 


loop ending at statement 66. The W's are the three values of the design 
variable being perturbed. W(1)is the smaller of the two values for the 
| comparison decoys. W(3) is the largest of these two values. W(2) is 
| the value for the basic decoy. The H®@LD array values are set in F123 
| in the D® ioop ending at statement 111, The values in the SAVE array 


were computed in F!23 in the D@ locp ending at statement 33 and then 


saved in F123 in the D@ loop ending at statement 333. The J defined at 
the beginning of the D® loop ending at statement 66 is the index of the 
altitude at which the values put into the R array were saved. JJ specifies 


which of the 9 possible arrays of differences is to be used. The statements 


r\ J following the definition of R(3) define the plot titles. Then IOP(76) is tested. 


if Consequently, if 1OP(76)=1 , the program jumps to statement 60 


where the AVCO plot routine EZPLOT is called to set up a plot and then 


statement 66 is executed. IF IMP(76}= 0, XL, XU, YL, and YU which 


| give necessary scaling information for the Aerospace plotter are calculated. i 


PLT is then called three times to produce a plot and statement 66 is executed 


next. The statement after 66 returns control to INTGRL, The sequence 

of statements beginning at 70 and ending just before 110 produce one corridor 
plot, Three curves are plotted on thesame.graph, Xve. T, Y vs. T, and 
2Zve. T, Xis the expanded array of lower corridor values. Z is the 
expanded array of upper corridor values. Y is the expanded array of deltas. 
| These deltas are calculated in F123 in the DO loop ending at statement 33 

io} and expanded in INTGRL. TT is the expanded array of altitudes calculated 

i in F123 in the D@ loop ending at statement 50, XL and XU are respectively 


the smallest and largest values in the T array. YL is the smallest element 
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in the X, Y and Z arrays and YU is the largest element in these three 
arrays, These values are needed by the plotter for scaling purposes. 
The plot titles are defined next, then at statement 95, IMP(76) is tested. 
lf itis zero, XL, XU, YL, and YU are changed as required for the 
Aerospace plotter. The 3 calls to PLT are then executed and control 
is returned to INTGRL. If IOP(76) = 1, statement 100 ig executed i 


next, then control return to INTGRL. 


The sequence of statements beginning at 110 and ending just 
before 230 can produce any of 28 possible plots depending upon the | 
meaning of X, Y, KY and IC(1). | 


If N= 1, a plot of Y vs. X will be produced. If N = 2, a plot of 
log to the base 10 of Y vs. X will be produced. At statement 110, N 
is set = 1, thenif IC(1)= 0 and KY = 4, Nis changed to 2. The DD 
loop ending at 111 then tests to determine if all of the Y's are ecual, 
If all the Y's are found to be equal, no plot is produced and control 
returns to VIXEN. If they are not all equal, statement 120 is executed | 


next. The values of Y are then put into the W array in case it proves 


necessary to take logs. This will be necessary if the Aerospace plotter 
is called to produce a log-linear plot. The log if each W is taken in the 
D@ loop ending at 122 if necessary. Statement 125 and the subsequent 
write statement put the plot title into the B! array, Then if IOP (76) = 1, 
statement 130 is executed calling FZPLOT. Control then returns to 
VIXEN. If LOP(76)=0, XL, XU, YL, and YU are calculated, and the 
Aerospace plotter PLT is called to produce a plot. Control then returns 


to VIXEN, 


The statements beginning at 230 close the appropriate plot file and 


terminate the computer run, 
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5. Other Information 
A, SUBROUTINE AVPLOT is called by either SUBROU .iNE INTGRL 
or SUBROUTINE VIXEN. 
B. SUBROUTINE AVPLT calls the SUBROUTINE MAXMIN. 
C. SUBROUTINE AVPLT calls the Aerospace routine PLT and 
te following AVCO library subroutines: 
1, SUBROUTINE IDFRMV 
2, SUBROUTINE EZPLOT 
3. SUBROUTINE FRAMEV ‘i 
4, SUBROUTINE ENDJ@B 
5. SUBROUTINE PLTND 
» 
D, SUBROUTINE AVPLT calls the following IBM functions: 
1, DMIN1 , 
2. DMAX1 
3, DL@G10 
E, Statements of the form WRITE(9, XXXX), where XXXX is a 
format number, place information concerning plot titles on tape 9, This 
information is read back by statements of the form READ(9, YYYY), where 
YYYY is a format number, 
« 
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SUBROUTINE MAXMIN(NTABLE, NDIMEN, X, XMIN, XMAX) 


1, Purpose 


MAXMIN finds the largest value and the smallest value in the X 
array containing NITABLE numbers. The smallest value is placed in 


location XMIN and the largest is placed in location XMAX, 


2, Input 

Name Source Description 

NDIMEN AVPLT adjustable dimension for X array 
NTABLE AVPLT number of elements in X array 
X, XDIMEN AVPLT array of numbers to be searched 
3. Output 

Name Description 

XMIN smallest number in the X array 

XMAX largest number in the X array 


4, Numerical Procedure 


Both XMIN and XMAX are set equal to X Then X., X 


ir 2 Xy--: *NrABLE 
are compared with XMIN. If any of these numbers is smaller than 


XMIN, then XMIN is set equal to that X. Parallel to this Xo X, ye XN TABLE 


are being successively compared with XMAX. If any X is greater than 
XMAX, then XMAX is set equal to this X. 
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5, Other Information 


A. SUBROUTINE MAXMIN is called by AVPLT 


B, SUBROUTINE MAXMIN calls no other subroutines or functions, 
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JBROUTINE PLT 
1, Purpose 
SUBROUTINE PLT is a dummy subroutine included so that the calls 
to the AEROSPACE plotter can be simulated. At AEROSPACE, this should 
be removed before the deck is compiled or it will replace the AEROSPACE 
plotter, 
2. Input 
* None 
3. Quiput 
. None 
4, Numerical Procedure 
When PLT is entered, control is returned immediately to AVPLT, 
5. Other Information 
A, SUBROUTINE PLT is called by AVPLT only. 
B. SUBROUTINE PLT does not call or reference any othersubprogram, 
i 
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8.0 LIBRARY SUBROUTINES 
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8. 1 Arc Cosine Function - AC@S (AC@OSR, ACGSD) 
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Are Cosine Function 


a oneotoiandl 


xX = AC@SR(W) xX returned in radians 


Kk = ACOSCW) X returned in radians 
X * AC@SD(W) X returned in degrees 
| 
' 


In all the above cases, W must be between -1,0 and +1.0, Alsa, W must 
be given in REAL*8 and X is returned in REAL*8, 


EXPLANATION: (Method ) 


J 1-R° or 
sing 
L-W* or 


sin@ 


I 


-W or Wor 


“COBY cos 


The arc cosine is computed from the arc tangent i.e. 


Cre arceos (W) = arctan a 


Since fiw? is always positive (because the Fortran square 
root function always returns positive numbers) and W can be positive or negative, 
the angle returned will be between O and TT {as shown in the graph). 


NOTE: This routine was compiled under @PT = 2, 
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My DeBolt IBM 360/75 __ 26 June 1967 _ | ASIN(W) ,ASINR(W) ,ASIND(W) 
my 
| 
LDENTIFICATION: Arc Sine Function 
USAGE: 
| 
, X = ASIN(W) X returned in radians 
X = ASINR(W) X returned in radians 


X = ASIND(W) X returned in degrees 


Im all the above cases W must be between -1,0 and +1.0. Also W must be given as 
REAL*8 and X is returned as REAL*8, 


EXPLANATION: (Methdd) 


Wor sin’ 
-Wor + sin” | 
1 


W 
= ARCSIN(W) = ARCTAN i aed ) 


Since fi >We is always positive (because the Fortran square root function always 
returns positive numbers) and W can be positive or negative, the angle returned will 
be either in the first or fourth quadrant (as shown in the graph), The angles in 
the fourth quadrant are returned in their positive form because if the angles are 
negative, they are added to 277 . ald example for the angle -TT/6 the answer is 
returned as 


NOTE: This routine was compiled under @PT = 2. 
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LDWENTIFICATION: Arc Tangent Function 
A ® ATANGCY,X) A is returned in radians 
A ™ ATANGRCY, XK) A is returned in radians 
A ™ ATANQDCY,X) A ie returned in degrees 
X and ¥ must both be REAL*8 and the answer A is also returned as a NEAL Ae 
lf X and ¥ are both equal to zero, the program will verminate and an error will oe 
given (1HC2651). | 
| 
EXPLANATION : | 
The arc tangent is computed by using the Fortran supplied routine DATAN?, J.a. 't | 
finds the angle whose tangent is Y/X (arguments as above), The oniy<¢ ffi rence : 
between these routimes (DATAN2 and ATANQ) is that ATANQ always relurvns Lhe angle 
in its positive form from 0 to 211. It chetks to see if an angle 1.8 negative aud | ‘ 
if it T#j it is added to ZT. 
{ 
NOTE: Thies routine was compiled under @PT = 2. . 
| 
' 
| 
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8.4 FIXIBCGM (MIBCOM)} has been dummied out of the program. 
Originally this routine was a modification of the IBM system 
routine IBCQM. 
* 
a 
I-535 
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it was related to the error traceback of MIBCOM, 
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8,5 FERR@R has been dummied out of the program, Originally 
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Li 4 IDENTIFICATION: 
i | i Title: Variable Field Input Subroutines | 


| | PURPOSE: } 
| — 


To perform the same input operations on 03/360 as BCDC@N/SYMBLS* does 
on the 7094. ‘These routines are also more flexible in the inputting 
i of subscripted variables. 


| USAGE: | 


rentitedeamrrat 


Calling Sequence: 


CALL WHERE (TABLE) ; 
CALL SETUP (8HNAME , L, XY, Nl, N2,....) 


CALL READIN (INC@I1) ‘ 
CALL ESDING (NTAP2) I 
Complete description of procedures for proper use of these subroutines: \ 
L \ 


location of a table which will be generated by SETUP and used 
by READIN. The size of the table must be set by the programer 


| 
| 
A. WHERE: This routine ig called to inform the routines of the |! 


such that the number of words in the table must be greater than 
or ecual to the number of input variables times seven. The beat | 
way to do this is to use the explicit type statement: 


REAL* 4 TABLE (ND), 
where ND 2 7 times the number of input variables. WHERE must be 


called before any of the other routines sre, and should never 
he called more than once. 


B. SETUP: This routine is similar to BCDCQM in that it generates a 
table of 6 character EBCDIC names versus associated core addresses. 
Variable length information must be included in the argument 
list and subscript information 18 included. as arguments rather 
then in a hollerith field, as is the case with BCDCON. It is 


* Reference: Writeup F2-21 of the Programmer's Handbook. 
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called once for each variable as follova: 


The first argument ie the extermal name of the variable and must 
have at least eight ‘characters; however, only the first eight 
characters will be used as the table entry, so any symbole greater 

than: eight characters having their first eight the same will 

| generate the same symbol in the table and only the first one 

| encountered will be effective. If the name has less than eight 

characters, it must be left adjusted and the eight character field 

; filled with blanks. ‘ 


The second argument, L, is the length of the variable, in bytes; doa. 
1 for LAGICAL*1, 


2 for INTEGER*2, 


4 for LOGICAL*4, INTEGER'4, REAL, 


so -- 8 for REALS, - + “vee 


The third argument is the internal veriable name (i.e., the core 
location with which the external name is to be associated). This 
can be subscripted to associate the symbol with core locations 
other than the beginning of an arrny. 


l The remaining arguments are the subscripts. (These are omitted 
1f the variable is not dimensioned), Except for the right most 
: mubscript, these must match.the diméneiods of thé variable 
exactly. The value for the right most subscript need not be 

l correct, but a value must be given in the argument list. ‘The 
maximum number of subscripts 1s seven; if more are given, only 
the first seven will be used. ‘The array length is also some- 


what restricted: e.g., if a variable has n subscripts then the 
4 product of the left~most (n-1) of them must be less than or 
equal to 32767. Yor @, & singly-subscripted variable my 
y 

! 


have dimension 50000°(say) but « doubly-subscripted variable 
may not have ite first subscript greater than 32767, although 
its second subscript may be as large as desireds « 


This routine, unlike BCDO@N and BALEAM e distinguish between 
single-valued variables and one -dimensi arrays. This can be 
eeen ty the following: 

Given: CALL SETUP (Gur , 8, X, 10) 


CALL SETUP (Sux , 8, Y) 


— 
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and the data card: X 5.0 3.0 Y 8.0 6.0, the equivalent Fortran 
statements are: 


X(1) = 5.0 


X(2) = 3.0 


Y = 6.0 


Therefore the &.0 will be overwritten, instead of the 6.0 going into 

the next double-word location above ¥. Incidentally, if Y were sub- 
seripted on a data card, READIN would terminate execution with an 

error message since Y is not supposed to be subscripted. | 


SETUP may be called anywhere in a program (even in other sub- 

routines) as long as it 1a called after WHERE has been called 

and before READIN is called. However SETUP must not be called ® 
more than once for any one variable since the calling instructions 

are not overwritten as BCDC@N does. 


The only error ie if WHERE is not called before SETUP, in which i 
ease SETUP 1s never executed, and when READIN is called, P 
execution will be terminated with an error message. | 


C. READIN: ‘This routine reads in the input data (from Fortran Data 
Set Reference No. 5, so the "data" card is K «FTOSFOOl DD *) ’ 
ligts the cards on the output tape (FIO6FOO1), obtains the location 
where the data is to be stored by finding the symbol punched on | 
the card in the SETUP-generated table, converts the data as pre- | 


scribed, and then stores it into the core location obtained above. 
This process is continued until® 1-9 punch in column one is 
encountered, at which point this number is stored in INC@LL, tke 
rest of the card is processed, and control is returned to the 
cailing program. 


NOTE: A zero punch in column one will set INC#L1 equal to zero 

but will not terminate a case. If the routine encounters an E 

(part of END-9F-J9B) in colum one, execution is terminated | 

inmediately by calling EXIT. 
| 
| 


Either set of characters (i.e., either the 029 or 026 keypunch) 
or a casbination of both may be used on data cards. However in 


albanien. fields, if EBCDIC characters (such as parentheses, 
equals, apostrophe) are desired, the 029 keypunch should be 
used. 


“soem een ery” 20) 0/0 hominy mises lai MA A | 
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The breaking character used to separate items (fields) on the cards | 
is a blank. Therefore the way in which data is presented to the | 
subroutine is completely varieble as it appears in colwms 2-72, | | 
with the exception that a symbel or a piece. of data may not be | 
continued from one card to another. The items on a card may be i) 
divided into two categories, symbols, and data fields. 


1. Symbols. These define core locations and are of two types, 
name symbol. and subscript field. 


a) The name symbol is a set of characters, the first eight of | 
which are used to look up the entry in the table generated | 
by SETUP. The name must begin with an alphabetic char- 
acter (or $) and is terminated by a blank or left paren- 
thesis. From the table entry the base core location, the 
variable length, and the dimension information is obtained. 
Execution is terminated with an error message if a synbol 
is not found in the table. 


b) The subscript field, enclosed in parentheses, modifies the } 
core location defined by the last name symbol and sete the 
addressing increment for successive storing of data. 
Maltiple subscripts in the field are separated by commas, 

z and as with 7094/SYMBLS, blanks may eppear anywhere within | 
the parentheses. Data may be stored sequentially (i.e., | 
varying the first, or left most, subscript) or by varying i 
any one of the other subscripts. If the subscripts are i 
not given, the data will be stored se La. note this 
ie different from 7094/SYMBLS). If any r of sub- 
scripts (less than or equal to the number of subscripts 
as dimensioned) is given, the data will be stored by __ 
varying the right-most subscript (the same as 7094/SYMBLS). 1. 
However, storing on any one of the subscripts may be | 


iH 


forced by putting the subscript number tely after 
the right hand parenthesis cf the subser’ * | 


Examples: DIMENSIPN X(4,6,3) i 
1) X 520 6.0 8.0 
x(1) 5.0 6.0 8.0 are all equivalent to | | 
X(1, 1, 1) 5.0 6.0 80 
? X(1, 1, 2) = 5.0 | 


tee X(2, 1, 1) = 6.0 
ry 4 
ee Tie or, aig og 
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i X(3,.1, 2) 8.0 | 
| i 
| 44) XGQ,1) 50 6.0 60 7 
| | r are equivalent to | 
| x(1,1,1)2 5.0 6.0 8&0 | 
x(1, 1, 1) = $.0 1 
| x{l, 2, 1) bad 6.0 
X{1, 3, 1) «8.0 | 
444) X(1, 1,1) 5.0 6.0 &0 
X(1,1,1)3 5.0 6.0 6.0 are equivalent to 
| or even X(1)3 5-0 6.0 8.0 
| || 
| | x(Q, 1, 1) = 5.0 ‘ i} 
| X(1, 1, 2) # 6.0 } 
} | } x(1, i; 3) = 8.0 | 
| | This xeally comes in handy vhen one wants to store by varying the * i 
oad second subscript (say), starting at some point in the middle of the | 
a array, Bxemple: to start at X{3, 2, 2), In 7094/SYMBLS one would hi 
7 have to write X(27,2) . since 27" 4x6 x1+ 3 1 
ie | and mist be c d from the dimensions. Now one only needs to | 
| write X(3, 2, 2)2 . i 
/ | As in 7094/SYMBLS there is a provision for using subscripts to define i 
Lah! /erens to be set to a value defined by a piece of data which follows, | 
/ Example: | 
i! 
! X(3, 2-4, 2) 5.0 is equivalent to | 
| X(3, 2, 2) # 5.0 | 
i X(3,, 3, 2) * 540 
| x(3, 4, 2) = 5.0 { 
i Any of the subscripts may be of the form Kl ~ K2 (where Kl < K2). | 
int Also, if another data field appears next, it will be atored at the 
AN location corresponding to K2 +1, with the other subscripts the same. 
il| " : 
| vhs i Saelld C ceal cailanl : 
Heed A 
i = lina 1-542 ~ , 
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Example: If "6.0" follows "$.0" in the previous example, it will be 
stored in X(3, 5, 2). 


If the number of subscciptse on a date card exceeds the number given 
in the argument list of SETUP for that variable, execution is 
terminated with an error message, 


2, Data Plelds. These can be divided into eix items: 


a) Integers are defined by lack of a decimal point, B, D, or an 
imbedded plus or minus sign. They can be positive or negative, 
and a plus ia assumed if there is no sign in front. For a 
length 4 variable the value must not exceed 2147483647 (231 ~ 1) 
and if the value exceeds 32767 for an INTEGER*2, the value 
stored will be modulo 32768. Also one must be careful when 
reading an integer into a length 6 variable since the second 
word of the double word will not be zeroed. 


—————S = = ~~ 


r >) Floating point numbers are defined by a decimal point in the 
field. However the decimall point may be omitted if an E or 
D, or imbedded plos or minus or both aypears in the field, as 
XXEYY, XX+YY¥, XXE~-YY. 


in In this case the decimal point is assumed to be just to the 
left of the exponent field. 


The maximum field length to the left of the decimal point is 
nine digits. Of the total number of digits in the field, only 
the first sixteen are converted. Any digits after these are 
ignored until an exponent field (Z, D,+, or -) or a field 
terminator (blank or parenthesis) is encountered. An E, D,+, 
or - defines the exponent field, and of course the exponent 
must be integral. Note that E and D are interchangeable since 
all conversion is done in double precision (only the store 
depends on the length of the variable). In gener, the per- 
missible forms are 


' tf X. OOOOCOOOKOOOKEYY 


+ X.so00NOOOOOOOE fry 
£ X. XOOONCOOOAXAAETY 


£X.XKXKXKXNXXKEKK, where there can be one or more X's, there 
may be either ons or two Y's, the deciwal point may appear any~- 
where to the left of the tenth X (measured from the left), 
and if no sign appears, it ia assumed to be positive. 
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¢ } Hexadecimal numbers axe input in a manner analogous to the way 
octal numbers are input by way of 7O94/SYMBLS. A@ (or X) In | 
column ore of a card signifies to READIN that the card contains 
only hexadecimal data and nothing else. ‘The data may be read 
into length 4 or length 8 variables only. A maximum of 8 digits 
may be read into length 4 variables, and a maximum of 16 digits | 
may be read into length 8 variables. If more than this maxi- H 
mun is given (before a blank), the first 6 (or 16) will be stored 
correctly, and the digits left over will be right adjusted an. 
stored in the next core location as defined by the last core | 
defining symbol. If less than the maximum is given, the field if 
‘da right adjusted and stored, and if other than 0 - 9 or A -F } 
appears in the field, execution 1s terminated with an error ii] 
message » 4 


4) Alphanumeric data may be defined in three different ways: 


1) Hollerith card for HEDING. This type of card has an H in I 
colimm one. Colwms 2 - 60 of this card (not 2 - 72 as ii 
with 709%/SYMBLS) will be read into o specially set as‘de - 
block of storage as Hollerith characters. Nothing else can 
be done with the card. 


ii) Alphanumeric data intended for arrays may be read in using an | 

A-card, ‘There are two cards associated with this option. 7 

The first card must have an A in column one, a core location i 
defining symbol and a length of field count anywhere in ; 
colume 2 + 72. This count is in units of the length of the dl 
variable. For instance, to read 24 characters the length | 
would be 6 for a Real*s variable and 3 for a Real*6 variable. il 
fhe second card, which must follow immediately, will contain 
the alphanumeric data to be stored. The characters on this 

| card, Jateg noe one will be stored until the length 

| ; of fi ry ted. The rest of the card is then ignored. 

| If the number of characters as cbtained from the A card 

exceeds 72, execution is terminated with an error message. | 

1. Example of A-cards: i 


on the 7094, A TITLE 12 would cause colums 1 ~- 72 of the 
next card to be read in. On the 360, 


A TITLE 9 (TITLE Real*8) : | 
3 or | 
, A TITLE 18 (TITLE Real*+) ! 


would aleo ceuse colums 1 - 72 of the next card to be read in. 


i 


| 
el 


i 


| 


it 
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141) Alphanumeric data intended for single storage. This type 
of data has a alash at the beginning of the field, and 
the first blank encountered terminates the processing of 
the field. The maximm mmber of characters is 8 for a 
length 8 variable and 4 for a length 4 variable. If more 
than the maximum are given on the card, execution is 
terminated with an error message. If less than the max 
are given, the characters defined are left adjusted und 
the remainder of the word (or douie word if Real*&) is 
filled with blanks. Also, to fill X with blanks, one 
can use 


e 


x / 


e) Logical data: To input logical data, all that is neces- 
sary 1a to use a decimal point followed by T or F; all 
characters after this are ignored until a blank isa 
encountered. For instance, the following date cards all 

Fr will produce the same result 


i | 


One must be careful when reading a logical constant into a le- 
ngth 8 variable, since the second word of the double 
word will NOT be zeroed. 


f) Table aig ss ale data. This feature is exactly the same 
as it is in 7O9K/SYMBLS. I. e-, data may be generated 
with a minimum of keypunching when numbers are desired 
ranging from a Lower limit to an upper limit in steps of 
certain given deltas. ‘The format is as follcws: 


This will generate from X, to X; in steps of Xp , 


Xy to Xe in steps of X, , etc. 


where the X, may be floating point or integers, but not 
mixed, There may be as many deltas as desired and they 
may be positive or negative. There ld be no 

' imbedded blanks in the expression, bu y are permitted 


" 1-545 


\ i wn nea eA rial al iI HINARI PHP “ee I GP H Hh shbiaC Hinson At dopa ririnennterrinentsntrh=eemmmnsntinaall | | 
VAN I i, ii 
| il Mi | y (hi 


bid 


(nm 


ni ttn it 


gs /360 PoPtran 
IV and BAL " 


=e Tangange: 


before a number and after a parenthesis 


io®., x, ( Xo) X, ia permitted. 


3. Tranefer card. Thies card will cause return to the calling pre- 
arem after the remainder of the card has been processed. It 1s 
signified by a 1 - 9 punch in colum 1. Also, this value 1s 
stored in INC@LL. 


End of job card. This card signifies to READIN thet there is no 
more data tu be processed, and that the caller wishes to terminate 
hie run. READIN will then transfer control to EXIT. 


Intexynal and External Controle. 


A) A * or $ in colum 1 will cause the entire card to be ignored 
(it 1s printed on the output tape, but NOT on line). 


B) There is no ®PRINT or *NPPNI option as there is in 7094/SYMBLS. 
All data cards are automatically printed. There 1s also no if 


CONMAN option. 
C) There are no external controls. 

ERRORS: Execution 18 terminated immediately with an error message, if . 
1) SETUP has not been called before READIN. 


14) A name symbol on a data card is not found in the table generated 
by SETUP. 


411) A field contains an illegal character (such as a decimal point 
in *. subscript field). 


iv) Aeariable is given more subscripts than allowed sccording to 
the table. Also, if an error is found, processing is term- 
inated Lmmed ia . The reat of the data cards up to the 
transfer | are not scanned for errors, as ie the case with 
7094 /SYNBLS. 


D._ HEDING 


This routine will write (on NTAPE) the hollerith card read by READIN 
and eject a page. If an H card was not read in by READIN, the printer 
will merely restore = page. 
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' IV and BAL 
: The only error is in the case that the Fortren data set reference 
; number (TAPE) is not defined on a DD card, in which case the 
| result is the same as if one tried to execute a Fortran WRITE with 
the same number. 
| 
Summary 
1. No complex variables allowed. 
2. Either keypunch may be used. 
3. Core storage: 292A, ¢ bytes or 26351 words. 
4, External subroutines used: 
a 
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TBCAM# 
| 
! FOXPL# 
. | EXIT 
\ 
I 
i" | 
F : 
hs 
{: 
ls 
3 
Is 
i? 
J 
| 
ith 
yoo" es ee etic, Sica 
- peaaa i 
a 1-547 
H| i shy atest \ H TDM VaBtal PD Pe TR PP AU A Tetons Ii | 
| alli \ 


ania bpuatietnaietnetaainemmameeteenammbanonadammenssemmeaameemttai ee ete ee 


" . th 9 4s , : | i 
Mees | reezarch and cdvanced development division | 4-005:1 iat 


| Pe ee ene oe PROGRAMER’S HANDBOOK _ Pe) eee 
| |PROGRAMER MACHINE TATE | BE HEY IERE SEPQL 
7 Gd —Aukingon, R430 1 OS 360 Te cite ayer 198 Ve a REA DENMEDING 
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| | 
| | IDENTIFICATION 
| | 
| Tiie’ Varsaole Field input Subrou:ines (Addenda: 
| USAGE | 
This wisie-upis an attempt to clarify some po:nts in fhe prev. cus 
write-up (F4-005) and also a discussion of twe modifications recenily 
mace to READIN 
A CLARIFICATIONS Referring to page one of the previous wr. ¢- | : 
' up. inthe discussyon of WHERE ND must be greater than or | | 
equai 10 seven limes the number of inpul variable names’ » e | | 
. greater than cr equal tc seventimes the num of calls to | | 
_ SETUP | 
I! » 
| The first sentence in the fourth paragraph on page :wo should vcad ; 
"The remaining arguments are the dimensions". ; 
in 7094/SYMBLS a dollar sign anywhere on a data card meant ‘ha: at 


the rest of the card was acomment This is NOT true for READIN 
A dcliar signis treaied as a letter, except when 14 appears in «clumn 


one in which case the entire card is treated as a commen: and 
ignored 


module which has been compiled under OPT=2 The compiler wil] 


{ 
Care should be exercised when calling these routines in a source | 
l 
invariably compile the module incorrectly unless 
! AJ) mput variables are in COMMON 
4 


2 The module does no arithmetic calculations but instead « Ils 
other subroutines whi 1 perform the calculations 


Therefore unless one (or both) of the above two conditions 18 sat- 
isfied the programmer should not use OPT=2 for the module which 
calls SETUP and READIN 
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READIN has been modified so that it will treat a comma and 
an equal sign (from either ah al pik as a field terminatcr be 
variable names (symbols) subscript fields, and numer:ca 
_ However, only a blank (or end of maar vermina'es 
alphanumer:c and hexadecimal data, This was done to enable 
READIN to process NAMELIST data cards, The only changes 
that need be made to the data cards are: 


1} The asterisk (for repeated ‘ields) mus. be removed sine 
READIN wall treat the field as illegal, In leu of the asterisk 
the dash may be used inthe subscript eld as follows 


= 200 0 becomes A(1-20}= 0 0, | 
2) The quctes must be removed from alphanumeric data ard a 


slash (or A-card) used instead, as NAME: ‘ALPHA’ te- 
comes NAME=/ALPHA 


3) Logical data must have a decimal point in front of che T or F, 


An option of returning control to the calling program afrer 
encountering the END-@F-J@B card has been added This has 
been done so that the programmer can regain control in order "¢ 
close any data sets (such as the plots, by calling PLIND: which 
woulc otherwise be closed abnormally by the system However, 
the programmer must be extremely careful not to try to read 
“tape! 5 afterwards After calling PLTND for whatever} he 
merely needs to CALL EXIT. 


This option can be exercised if the programmer adds a secund 
argument inthe CALL to READIN. whichis a non-standard re: 
turn and sets INCOL! (the first argument} to the value 8888 upon 
first entry to READIN, This is a flag to inform READIN sha there 
18 a non-standard return An example follows: 

DATA INCOL1/8888/ 

CALL SETUP [etc } 
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DATE SCO NAME WHE ei ETUP, 
| DIN, me 


1 GALL READIN (INCOL1, $99) fora99 if EBCDIG! 


(Calculations) 
GOTO 1 
99° CALL PLTND 
CALL EXIT 
RETURN 
END 


It ehould also be mentioned that no changes are necessary to existing | 
sgvayeme which do not use this option as long as they do not set 
INC iL to 8988 initially. 
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IDENTIFICATION: 
Title: Variable Field Input Routines (Addenda) 


Kl \ 
fh : 

READIN has been completely rewritten in BAL in order to increase efficiency 4 | 

and to save core storage, In the process several modifications have been made. " 


(1) Special characters are now permitted in the names of variables; only mt 
the first character need be alphabetic or a dollar sign. Any special characters, i | 
except an equals sign or left parenthesis (either keypunch), may be used, | 


== 


(2) A blank is now permitted between the E or D and the exponent in 
floating point data fields. READIN treats the blank as a plus sign. Thus ' 
READIN is now able to read cards written by a FORTRAN program, q 


(3) In logical data fields, . Y and .N (for true and false} are now permitted i HI 
in addition to. T and. F, Therefore the user can have logical print, iaphi plot i! 
flags in his program to control the output, aa follows; | 


PRINT , YES PUNCH .N@ PLOT .YES a | 


As with.T and.F, all characters after the Y or N are ignored until a i, 
blank or comma (or the end of the card} is encountered, ' ih 


(4) The suppression of printing of the data cards is now possible by punching 4 

N@PRINT in columns 1 through 7, To restore the printing mode, the user | 

| punches PRINT in columns 1 through 5, Anything else punched on these cards 
| is ignored, The user does not need to start with a PRINT card since READIN i 
is normally in the printing mode. | 


(5) READIN normally ejects a page for each case. This may be changed 
to a line skip for each case by punching N@FLIP in columns 1 through 6. To j i 
restore the page ejection mode, the user punches FLIP in columns 1 through 4, 
The remainder of these cards is ignored by READIN. 


(6) The core storage requirement of SETUP/READIN haz been reduced to | 
LIBB ¢ bytes or 1135 single precision words, a 


|More H ‘ rp ISPAe DFP HTH ep ayteetenrcmerrteeirominsete eaten AN, 


tl 


a 


The dummy subrouiines EZPLOT, ENDJOB, FRAMEV, I. 
and PLTIND replace the AVCO plotting routines of the same 


names and must be removed when the program is used at 


AVCO. 
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CHAIN - MAIN 


¢ AEN : 
ot casas tic ech seh, woh aweedsosanioac ospenocra a jetted sstsufar GiobaioaeS it Gite deo fisGbaak sr Sie naan de jamais neon MORO 
¢ HAIN 
Se Seis acinus wee eeneten Reagan Aonuma teas EMA 
¢ ‘PURPOSE. MALM 
a Ee SRE Ree ee eT Ee en Oe 
C MAIM OLRECTS THE CALCULATION OF THE OPTIMUM VALUES FOR THE MAIN 
0 for omen SIGN PARAMETERS. OF A DECOY. —. THE TECHNIQUE USED 15 10 a 
C MINDMIZE THE DIFFERENCES GETWEEN CEATAIN CHAAACTERISTICS UF A WAIN Lh 
+ Cnn KLEEN TRY VEHECLE AND SHE « SAME CHARACTERISTICS FIR A DECOVS MAIN Law 
t THE FOLLOWING NINE CHAHACTERISTICS MAY BE USED = VELOCITY, MAM hd 
vind Gos minora DE GCELEMA TION, G4LK09TIG GJEFFICTERTs WAKE LENGTH AT EACH OF MAL bo oe 
C THREE RADA FREQUENCIES AND WAKE CROSS SECTION AT EAC’ OF THESEMALN i> 
oe Poe wow FREQUGAC TES « FOUR METHOOS FUR FINDING THESE MINIMA ARE MAIN he 
"7 (1) A FIBOMACCH UNE VARIABLE SEARCH (2h A FIBONACCI Two MAIN LT 
ws Coven we toe WARK DALE SEARCH F391. THE DAVIDUN METHOD (4) THE ROSENHAULK ALN dd wn 
C ROTATING COORQINATE ME THOU) A QUANTITY PO, THE PXOBADIEL@~ Malm 29 
Came Y THAT A DECUY WILL BE DISCRIMINATE AS HeLNG A DECUY,s CAN «Main . 20. 
C ALSO BE DETERMINED, MAIN 2h 
ee : : MAIN Bc ow 
c MAIN #4 
ie (oa oat Benet es ee ia - MAIN 24 
C MALN 25 
ee er Si Eas a es wee «Ape gn UN MEE Wescorotececs 
COMMON /XKSAVE/ MISAVE 140}, X2SAVE(4O) » KOUNT MAIN 2T ‘ 
COMMONS TOPT/ LPKOC oe LNy NCONS pI PNT, LEMS LIMIT TAAND wee | MAIN 2d 
COMMON/FOPT/ERR «PAANG, FAC, DELTA MAIN ¢Y f 
mown COMMONS MTNA AL OW (20) UP 120) AMULTO 20) CALOWEZO1CTPC20) —... -. MAIN 30. 
COMMON UCCUH, NOCTUR MAIN 91 
1. Me «COMMONS CCKN/O em . nwa 2601 A Me eae ta cables wl beeen MAIN. od Peers 
COMMON /TONOSS IDLL50}, 102150) MAIN 33 
+ ono COMMONS OWL AZOVERC TI 20} sWRFYLAED MAIN 34... 
COMMON/OLKI/ HHA 1 B00} y X440) 9GL 40) 9 $E40) pXP (409 GPCO) TCS 99 GBIGOI4ALY 35 . 
weoweronnm he Ge GSP y GTP GSSs GIT, GSHeFeFPoFUeFOrEs PeTOyASsShe LeQeAeELy MAIN 36... 
RUELZZ FAK eC 140,10) MAIN 37 
common mies CHIMP SENS EZ ES ENG g EG EME cee mets cena mime me cate ema TAIN BD ets wane 
CUMMON /NALTFG/FGS4yNALT MAIN 39 
miviioommmme COMMON FELCCURS LUC CURLER DD. wo. - ee wee . » BAIN 40 
COMMON SEND ITCAM MAIN 4) 
vom aeovcmetmemee GCMMUNA MENG EA Ke ALB C20) UBC ZO) KR CLON eB MINH cone sect 
COMMON/OO0PT/DEL KI 20) WAIN 43 
Eb ck halle Bs, SPR aerate os case te ase ae asia eS MADN 94 coe 
DIMENSION TOP(90) MAIN #3 
wim EMENS LON GUPT ZO) monn mnn a satte Sah eee cea ae etna Se SNE AAG) Gece 
sig ek OCCUR (40001 5 NOCCUR{ 301 MAIN 4&7 
Tee ‘ inte ail ai ieaieen ae HOEK Coe, a 
EQUIVALENCE {NOCCUR (14) ¢NPRINT? MAIN 49 
ntti EE QULW AL: ENCE 6 OC GURL OUI. Dg LPP nt ernment etme LA Nn On nem 
_ EQUIVALENCE FTOCCUR(QO0OL), OP} MAIN SL 
ior pereemsorenb UTA ENCE 11 OCCUR UDO RODTs. WPL OT banana cnn ee MAEM Og Hw 
EQUEVALENCE (IOCCURIOO301), IREF } MAIN 43 
womamooos BQUEVALENCE (LOCCUNIDO303 be MODE © Jon. ee et ee | MAE Om 
¢ MAIN 55 
ARGGIN tlh eon: > . Hite senha -it e Hetetmetsnacinninesvaitenanien aati ants neve tenn MAE BG mene 
C ITERM is DEFINED IN FEV i MAIN 57 
C TF ITERM @ whe FUNCTION = ZERO es ae meas ise MAIN 58 
CG TF ITERM = O, FUNCTION 1% NUT ZERO : MALN 59 
a 
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C OIPOITERM = Le FUNCTION 4S$--UNOER INGO 9———————— nen ee MAIN BQ noe 


C MAIN G2 
“ » CALL WHERE TAB 9.9) not soe ie i i nttcenecnnnane enn amet romeenins avicny vim UML Al ob eecence reese 
OCCUR 42) 3414159265 3569793 MAIN 63 
CALL $R2490 se ee creme tnt cies at mm vem ee ene tee siesta nes ete MATIN | G8 cremnsemnes 
GALL ZPRM MAIN 65 
MODE mem eee meen anemones med meieerens sence mote sve can oc TN Go ccnmes an 
Poe 0,00 : MAIN 67 
pe Nee eacecs ASO SO 0e x Ss Stacia ee rant nee ee enemies ratetamamtnee MiB | Nw 6B nneeern ne 
DEL TA#1.000 MAIN 69 
FC BT SOOO. sri iomt® ssn poner CORAM macaw OE ede Suwa tue a, = MAIN 70 
K=O MAIN FL 
Te NL oe Meets cmercses vine ec ele is ats sate ree eee Nhe a Sn, Se mum MAIN TF2.. 
TMPLOT#O MAIN T3 
nie enanosi ell RE OVD ooo ~ aceanteennnnrenennecnemenmit enters nines ewer en ect wan ante wen ae MALIN To, cesmeccon 
0 707 Tw |, on MAIN 75 
Ngo i einen Seer. eee Osu a3 fe MAIN 76 
707 to2eh) = O 4AIN TT 
c : anwaate Srnec sate arcu al iam Mla oy oe F ae ice MAIN To. 
C MALN 79 
wanes ov inne Li A EA OL TILL LD LMP LOL. area nent ee MAIN. BO. WW. 
ILLDLia2 MAIN b1L 
IFEMJOE2} LOOe1LO,120.. 0-000 20°. : aes MAIN b2 
6 MAIN 63 
100 =IKEFS = IREF a eee wwagens me Sees MAliy 84 
[REF = ] MAIN 65 
= NOPE DE) tO er weenie cere reenpereeeneenit eens eens eet een erage mews wees eee AU, BY cecceecene 
CALL FL23 (VAL) MAIN 87 ‘ 
RRC ee UGG ee cilia tate a wrctaincmetinens oe ern ieee MAIN 88 
GU TO 2 — MAIN 89 
1100 CALL FI230VAL) os See Te Seer ree bee MAME OBE is 
Go TO 2 MAIN GI 
Pic! ; ; arg Renn anatase dae ncn, ee Ge Be Se ere il Tee Gee menienn 
€ MaIN 93 
520  LEUTREP cE qe 2). GO TOR... ne eee sree ee ee eee 
CALL FL23(VAL) ' MAIN 95 
GO Tog + mecinetats ni penance ne tt en ntesdeae te beens e eesen: meee teste ices 9 MALIN POs oe cerns 
C MAIN 9T 
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CTPSAV = CTRL) eae MAIN 100 
I¥ = ID201) MAIN 104 
-VSAV = OCCURCIYD..0 ww ow i ee moat oh mcienice an cc) dn EE, + 
tr MAIN 103 
saeene penetrometer nna nero e an MAIN 10%. 
C SETUP INITIAL STATE , MAIN 105 
seem DOLL ID 9 TIN imei co set memes eco inaree meme at oe ie tenn teen ey eeeemenenetetom anes anv MAIN 100... .00- 
12s IDL(T) MAIN 107 
ornare oe MECURE [Z)eOVECTILS as ait hay) maim ream iee RD aN woes tiene woah Sim re seemed MTOM RNG Dee 
X(T )sOVECTU) MAIN 109 
eee Da KK ODL ROVECT OD nner MALN, LEO 
C MAIN 111 
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IFLITTERM.ZEQ.1) GO TO 200 . noo. MAL LLB 
A{L1} = XMIN : MAIN 119 
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C Fuh} & 
C FUN Y T 
€ FUNCTION FUNT COMPUTES THE PULSE SHAPE RADAR RETURN AS A 7 BUNT TAO 
G FUNCTION OF OTSTANCF SFHIND THE LEADING FOGFE OF THE RADAR FUN 9 
c PULSE, FUN? 1f 
C FUND 2 
(cs FUNL 22 
SIG1= 0.9nN FUNL 13 
“= S3G2= O,6n0 (cD tee : Sn EASE awe S Ore eae fre aes a FUN 14: 
Stc3s 9,00¢ “PUNY 15 
SiGé= 9,9N0 FUNY 16 
MXTT= XN = XTT ‘FUNT 17 
C # FUNY 1A 
Cc Yl~= REGION AFFEC TEN BY IVE RDEN SE SCATTERING FUNT Lo 
TOL YL = OMINTCWSI,xODPY - EST Sree os) are ages CAEN Me ee et nes 
¥)o = OMINIC¥ IT, XXTT) TUNT PY 
WS? = ATT & XNMP GC WSL = XK FUN'T 22 
VP = OMINDTC¥ 1,052) FUN) > 
FEUMMTT OIA, FLA, TAT FUND = 24 
TOV IONXTT~ “HSA KOOP } 104 9 Tht, TOD FUN Y 9B 
c Ne Sse RMe RES? ott as cl tial eas oS FUND 96° 
C SUG?) = PARTITION OF PLL SF SHAPE NUE TO NVERNENSE RETURKA Fil = > 
704 SIGL = VL * OWST ® SpH? FUNL oa 
Whe WS? = PON & XLMS FUNT = 20 
c FUN T 434 
C YZ = REGION AFSEC TEN BY FURST PART OF MULTT PLE SCATTERING FUNT = 3Y 
YO = DMNINYOWST WS?) “FUNL 3? 
WS2 = XXTT = xrine FUNT 33 
¥2 = OMIN1tW52,¥2) FUNT 34 
WS? = XTT & XMSP £ WST =KX FUN = 95 
¥2 = NMINI¥2?,WS52) FUNE 36 
TEQXXTT = XMOPD ITS: Wb, 706 FUN E ees 
TOE TECKKTT = WSL = XMSPV TORS TNR TNS nen ee ra “FUN T "3R™ 
€ S1G? = PAE TION MF PULSE SHAPE QUE To CONSTANT PCRTTGA OF FUN t 30 
ig FUNT  4n 
C MULTIPLE SCATTERING REGION FUN. 41 
TOT? STG2 = N,500 * Y? * NWST & SPH? FUNT 4? 
705 WS3 = 1,500 * WS] FUNL 43 
eT TEC S DL SEX OK TT SMM SPD POR TOD TVG rn oe nt ee ee FUNT “4%” 
c PUNT 4&5 
C Y¥21 = UPPFR LIMIT OF REGEON AFFECTED ®Y SECOND” PART Che Ss FUNT 4&6 
ee 6 MULTIPLE SCATIFRING FUNL 4? 
TOR VIL © XXi= WSL = XTT FUNT 4 
69 TO P14 FUNT 4° 
TN4 ¥Y?1 «XM Sp * Pee a ae! Se nc Cane ne ar emp, es Pe ters neem cine Sd mss jee aon tie m=" LIN 1 “eM, 
The TPORK=ONTT & KSSP)I7116 7115 12, TUN. $1 
c Po aeees ; ia ~" FUNT | §2 | 
Cc ¥22 = LOWER LIMIT OF REGION AFFECTED RY SECOND PART CF FUN? 43 
C MULTIPLE SCATIFP ING FUN TL Se 
71h ¥3? © XXTT TUN. 55 
co tT F123 ss “a tle Merete c=, Bitint ase ars iia PUN” 58 °° 
Tle ¥37 = xSsp FUNT ST 
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‘ yl = GPP sb bap ber een pom Wet Creo WY SINGLY SCATLERING run] 
flh YS. * XSSP fib oJ 
cn Te Fe ae 

PLA W412 XNTTH~Ws? FUN 

c FUN 
G Yé2 = LOWER LUMIT NF REGION AFFECTED BY SENGLE SCATTERING FUN | 
716 -¥4? = XXIT FUN L 
UWixxtTy ~ XNSO) 720, 215 7?3 FUm 1 
T?71 IPCXXTT ~ CWS] 6 X88P)3797, NT) 729 FUN 
TAT EXY @{NExXR(-92e ¥2) mor xpt-nye YR? pp ® URAL. $00 FUN 3 
EX 2 =O AEXP (4.905% SLY VAL) =~ NEXP ($4,903 091 * ¥32)) FUN 
C FUN YT 
C Teo FOLLOWING TRST TS USED TA PREVENT UANESFL rw FUND 
LFEDARSLEX 2) ~ 2.CD= A) G07, N07 918 FUN 
OONT EXe = Aone FUN 
QuuP EX2 2 0,280 # JORG & EX? # 2CR4 FUN 
TFUXR 25) 72% P19, TLE FUN 
c = FN 1 
c S1G7 - SAORTION NF PLLUSE SHAPE DQUF TO VARTABLE PCRYLTON OF FUN] 
€ : “SULTEPLE SCATTER ENG REGION FUNY 
TAS §$1G2 = NWASTE ® CFXT € N2) & (37 = VAL} *& SPHS 4 89 @ Cow FUN 
Gi Th 729 PUNT 
T1A@ STH? = ONPSTE & LF XY «& EM?) & 2,990 & WS4 *LOE KAI -9,.5C0 * AL * FUN] 
1 MA QCRY AP Ye SENPPRO,See SL KIZ9eVI279® BY * SPH24942,009— FUNT 
LT X92 PAC LARN AT) aOWST FUNT 
T20  TEOCXXTT = XSSP) 730) 722, 722 FUN I 
T22 EXD = MOXPLHK4.% le FIM Val d= DEXPL-4,N9078 RIM ¥42) FUN] 
rKa = MEXMPL-AL® Y4aL) - DF XPL-N LE YH?) FUN LE 
c FUN 1 
Cc $164 = PIRTION CF RITURN OLE TA SENGLE SCATTERYPAG REGLCN FUN 
$154 3 Me? SNTeOWASTI® ZOR4SE ML e2C2% A ONWAAL €1S PHA 16,0 LI FUNT 
{ eMUST/IOW E RNY wICR% & FN? . FUN 
c FUN 1 
t SIG - PULS® SHAPF = RADAR PFE TURN AS A FUNCTION CF OOLST ANCE FUN 
C ALHING THE LEAPING TOGE NF THE RADAR PULSE FYNY 
TM: STG = S161 & SIG2 & S154 & SIGG FUN 1 
TFC ING 2= RRP OCS LLC RS VO S FUN L 

LAWS Wipe besaenrs) Vp SIG V7 STG? ¥3) Fury 
Ma & FORMATE LI» TH VL ay LPS Loe Feb) SIS L= pl RELS oT POH Y2= pIRELS oTyLOHFUNI 
t SIGPs, TPF TS. ty Ali Y31=,4 VPFLs.7} FUN] 
Lk Were te, 20908) Y22, VSL, V4? 4 IG} 4 3164 FUN? 

PaaS FIPMATO LON, TH V2%25 TPT] Se 7s (4 YL sgl PF 1S 67 6H ¥42" pLRELS pTy LOHFUN I 
! SIGT=» IPE LS ATO SIGS pL PF1S. 7) FUN 1 
Lethe sw@eire bey2cors ) XX+ SIG FUN T 
pave POMATE 2X, 2X =, IPL 15, T,6H S1f= pLPE 19,71 FUN 
yTS GO TO CAPO, LONp NM FUN 
FIIN 1 

ne ge ee FUN 
FUND 
RETUSN FUN T 
C ‘s TS: FLUIN 1 
WI FUNL = tsicuns - - “sigiera FUN LT 
C oy FUN 1 
RETURN FUN 1 
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satus coc abe ue onto Ud ner svn He A rho tyiteny! el oo THEA ETL ET FE il : ] 
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uties Pi a I 
SUNROUTING PL23CXxXX) 7 ase a dt oe Sef ee ee ere Gt’. Ae ee 
IMPLICIT 8 FAL®S (heHO-Z) ish net) 7 at 25 Ae Veale 
‘ ¢ Fl2a 4 
¢ FIPe I 
t SUMBOUTINE FL23 CAN UPTICHALLY = (1) READ DATA FOR A FI[h & 
c REENTRY VEMICLE TRAJECTORY FROM PNOUT. CR CALL SUBROUTINE VIXEN F122 7 i, 
C TH CALCULATE THIS DATA. (2) CALL VEXOR TC CALCULATE THE CATA Fist ke 
¢ BOR A OFCOY TRAJECTORY. €3) COMPARE THE DATA Efe A KeENT RY ryai og 
¢ VEHICLE TOAIFCTORY WET THE DATA FOR A OECEY TRAJECTORY Ab Pipa its 
C SAVE THE DTFRFRENCE S$. €4) CALL INTG2L TC CALCULATE awe PLUL FI?™ Lt 
¢ A CURR TDOR INTEGHAL, (9) CALL TNTORL TO CALCULATE Thr SIGMA Fiz 2 1? ; 
¢ SOUATE INTEGRAL TU KE USED WY SUBROUTENE EFFECT. (eo) CALL ° FI22 13 | 
ic INTGRE TO CALCULATE A TATLF Ck OEAPLUENCE COEFFICIENTS ANU TU f1??~ 14 
c PRANUCL INFLUINCE COEFFECICNT PLOTS AT SPECIFIE® ALTIFULDES. ~~ F12% 157 i 
C F124 YA \ 
C F123 17 i 
COMMON ACCURSNOCCLR ‘FUP 1a I 
COMMON ACICCUR/ TNCCUREI2C) FI24 49 i 
COMMUN ACPECUR/S PCCURATLI70) oo 1 
COMIN STIM CUMS TIM ee Ge Be a eee) Ma ieee aed, 
CoeMn ZIMETSTY cocAY F124 099 
“OREAL ® 8 LAL, LAZ: MOOT? May MX&e MYy MZ PO ML Yao oa oe j| 
DIENS LON ALMAX COUCH PALHINE 200) 98C4 (40) BCD (40) ACV ISO} (KEWL IL aE? FL?4 24 
CIMENSIUN BCWL 2040) ¢3C WLIL4 7) pBCHRL (40) ACW? (47) gp ACWR3 (40) “PLIR 98 
CIMENSTON = TUMATOU, Y10260) FL23 94 
~) CIMERSTON BETAPE CLEC) eB SELEO) COPING (7 6242) DALLAND DOG LoU) TT ORDO i} 
CIMENSTION USC LEG) DYE L460) ef VAR (50) sOVH(ED) sOVL (59) sDNLE (Loy) F123 7A | 
, C IME NS TUN Dil PCL EC) s EWE 3 (LOO) gl) WRT 1153) eVWR2 (161) DW 834 Loy) F122 Qa" 
CIMeNS ION FMAXOP°O) -EMING 2°99) oH t4 4 eHOL OCF 5) pI CR (90) s MOLT C4—4) Fl?3 tin 
CIM ESS TUN NCOMOV CEG) pNSCCUR CAC) NPVILEC) pFCCURI4IC Ty PEPS ELE) ‘T1235 310 
CTMCH SION Q90TC4s be SAVE ED el Oly Vd eSB C49) o SO44} SVI49) pSaL1 (40) TRV2A See? 
CLMENSION SWL 204 CV 9 SW 4049) p SHOAL (49) p SHR? (iid SWEF C4) pT AMAA (20) F123 430 
CIMENSHON TAMING 2 CU) ¢ TER C40) 4 TCO (4,) pTCV ES) p TOWLE (40) p ECWL2 044) F123 34 { 
2 CTE ES TON TOWL CGC) p TOWN LCOS) 9 TOWRA? Co) pTEWPS (69) pT HLL ou) FLYPA 3K | 
CEMENS TON TPLUTELEC) »TSCLGG) TIM C2 ru) VOLOT L603 VS EAU) FI24 94 ‘| 
LCIMENS ION VUGPLTCLEC) phLIPCLGO) pwLl SCL53) pWL2PCLA9) eWE2SC Low) Fl?23° 37 7 ia 
CEMEMS TON WE 2° CLEC) eWLISCL GG) eWREPCLS) gWRIS (LOD) pHRAPE Low) F124 38 t 
" CIMENS LGW’ WR 2SCLEC)» WR IPLLOAD phUBSLIAM) 92 MAXC29 DS PHING 2007" ES Raa ao 
LIMENSIPN PPLOTELEC) 280160) £123 0 4h 
“ EQUIVALENCE (OCCLECON22S5), ALALT } 3 aia ee RFLPY 41 
HQUIVALENCE (OCCUR (O1CG4), ALMAX ) ‘ FLIPS 4? 
“EQUIVALENCE COLCUP(CL244)5 ALBIN ) ais aia FI230 430° 
COUIVALFENCE (OCCLREICCEO2), ALPHA ) FI22 44 
EQUIVALENCE “LICCUR (LOCO 9) ALPRIMp OTE eter iment et eg ag 
CQUIVALENCE (NCCUR(LGL22), ALST } FL23 46 
" EQUIVALENCE (COCCURCCO1L28), CASE yp fi es “FL23 CAP” 
TFQUIVALENCE (MCCLOCACCSS),» COB >} F122 49 
is FOUITVALENCE (OCCLR(Q2753)—9 COFINFR © 00 0 Ummm eg pg 
TOUTVALENCE (COPCURLEC1LOS), COL ) F127 40 
TOV AL ERE E LOCC UP USEC SAN GNP: he 278 ee te no aoa eteane 
EQUIVALENCE (UCCUPECEICLY, COPA} peat FI23. 52 
EQUIVALENCE (NCCURCCC2r2), CM ) - ee EES SRV 599 
FAUEVALFUCE (CCCUTECOZCI3D, CN ) FI?3 4 
ERULVALERCE (OCCLRI606S4), CCODLAM) : “err ns “elaa 6a” 
HRUTVALENCE (OCCURICCZLG), COSZETD ncn Be ae 
PUWEVALENCE (OCCUR (40215) CTH2ZCT) : : FIP% 57 
FCUTVALENCE (NCCURIQCL27), DATE } FI2a oR 
TOUTVALENCE (OCCUR CEC235), DECHEP) 7 |” , : Fi2? 59 
ane : 
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j Mn Te VA AUOU PMU DIN NEUE ape Mea omega 4 A 1 ugg STOEL TU ee ee pve papoay (0g 1.00 aL tf ACMA 
iy == va ih 
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| ith tf 
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| | iI / | 
| :' FQUIVALFNCE (OCCUR (CO2tée), Ofer tC) Fi2? #0 | 
a COUVALENC® (OCCUR (CUP 4)» OF LOOPS biPa Al 
| LOULVAL FAGE EMCCLECCO2ESS. OTLY ene . FL29 A? | 
a EMTIVALENCE (OCCUR ICT2 74), PELZ EL?4 63 
ae BUIUEVAL ENCE CORCUPONCE2GY—@ FAN 1 FLD3 AG " 
nl! FUUTVAL ENCE COP CURECOIG7I« FACTAD) Se 
ml FUUETVAL NCE (OCCUR EOF P4S) « FMAX i 
En COUEVAL NOT (OCCUR CE2&4al) FAEN D f124% 47 
an LUULVALPNCE COCCURECEI CS). GARE) 7 EVI Y Ah 
Li it COU IVALTHCE (OCCUR EGIC2 A), GAMMA D F173 AN 
an , PIMEVALITNER EMCCURE COTS, LAT 2 e177 th 
ae SUUIVALENCE CUCCULCUtL aa LAR) Fi22 71 
ee VOUTVAL ENGR LOCCURLEPTF OD, YET 2 rl2. om? 
i | EQUIVALENCE fOCCUT COLT) MS } FI27 73 
EOUYVAL ENCE (MCCLURE NG2DTD, MX ) FIPT 74 
EQUIVALENCE (OCCUR COC7ILe MY ) Fi?% 75 
EQUIVALENCE CNCCLROCOP12Z$, 42 ) FIO% 7% 
LOU AVAL FNCEUUCC UK C40lt be LO) F173 «77 
CQUIVAL INCE COCCUREL2EC1), PEPSA ) TL? 78 
EQUIVALEMNCR COCCUREOCIL2), PHT) F122 7a { 
_GQUTVALFNCE COCCUPLSF22°3b» PSIZETY oo. oe ee FL27 0 ar q 
FOUTVALENCE (OCCURCCLE14), PSIG) F)23 AL 7 
TYUIVALENCE (UCCUREC 2676), ANAT 9 FI7a A? q 
EQUIVALENCE (UCCUPULCYIN, AF) F124 93 q 
FOULTYALENMCE (OCCUR(CECEG), RADUT FIP% fH if 
EQUIVALENCE (QCCUR{EC1I6), RAL) Fl22 95 ,. 
JROUAVAGONGE (CCLET OTS O Ts PR Ne oe 6 8 iat ie eo oe Lee 4 
FGUEVALEMCE (QCCURUIOCSS), @NDOT I FI?% a7 1 
CQUIVALENCE CQCCLP(GL135), RNI FI23 9A - 
CQUIVAL EMER COCCUNCCO141 9, ANZ) FL?4 99 » 
COUIVALENCE (MCCURCUCLIE), SIG  ) FL?2 90 i 
ECUIVALENCE COCCUPCUL221), SENGO } f12% 91 
SEQ UVALENES COGEUPEOGURINGCENBN A se oa liaitlaees ) aoEe Oe ae 
EQUIVALENCE (1CCUSCLC2TAle SUSZET) F123 93 iy 
EQUIVALENCE (UCCUR(GPA17)— STNZET) F123. 94 , 
FUULTVALENCE (UCCURIL1444)5 TAMAX ) F173 95 «| 
EOUTVALFNCE (OCCUR(CLE44), VANTIN} F122 96 q 
COUIVALENCE €UCCUP (03642) _¢ TRMAT } F123 97 
— BOUIWALPNCE (NCCUR( CU 201)» TH Tees fo eae des Ole Ee eae 44 
EWUIVALENCE COCCLR(UCLLID, THFALO} F124 99 iq 
EQUIEVALPNCE (OCCURCCC224)_ THEZET) aaa FI?3 100 a 
EQUIVALENCE COCCLPIGC266), IZTF ST) F139 1 . 
EQUIVALENCE (OCCUR(COCRED, BOOT ? F123 102 
EQUIVALENCE (MCCLUR(91644), ZMAX ) FI2" 193 ia 
AUT VAL FNCE (UCC UF et Eee Gl De i a i ence nen een see FIZ3 1°4 i 
EQUIVALENCE (PCCLR(CAIATT, BCR} F123 196 | 
FOUIVALENCE (PCCUR(KAICITs "COD F123 196 | 
EQUIVALENCE (PCCURIOGC2L), NCH) FLZ2 1A? | 
PAU IVALIERICE (PCCURT CEST Ds BCWLL Do rs — FL24 19K | 
BQUIWAL ENCE (PCCUT(CEEEN)» CHL? I Fi2Zz3 100 a 
EQUEVALENCE CPCCURTGGK4LDs FCRLS DO a i, 
EQUIVALENCE TPCCURECEZELD, 3CHRL F173 111 i 
FAUEVALENCE CPCCURCOGT4IY, BOW DO by Set echle de hile iH, 
EQUIVALENCE G€PCCURCEA471), BCRRI) F12% 113 1 
CYVIVALENCE CPCCURICCES)» BTTAPLA FB DD ; 
CQUEVAL ENCE (PCCUPCOIF414, BS ’ F123 115 i 
BCUTVAL ENCE (PCCUR(GZCALRs OR ce cen ee PERS LLG { 
EOUTVAL FACE (PCCUR(CLS21)» DO } ; F123 117 i 
FYUIVALENCE (POCUPCUTERT + OS } FI? 118 i" 
EMIVALENCE (POCURLLIL7ET)» OV ) . F123 119 i 
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TUMTVAL ENCE (PCCLPOLTTOT s fNAR OD}, FLI3 129 il 
LQUIVALERCE (POCURELIOO 7)» UV) FIPS 171 i 
“UEVALENCE (POCURCITERT « OWL} a eee il 
- LOUIVALENCE (PCCLPLO2TALds DRLE Fiza 173 i 
CAUEVAL ENCE CPCCURCOPPR LD, OLD . FIPA 124 : 
PEUUMAL ENCE (PCCUPLUAC41), MALT} F124 125 f 
CUUIVAL ENCE EPCRUROCOA4L)» ORFL FLA 196 ; 
COUEVAL ENCE (POCUR ECA4CL) » MWB?) FI2Z% 127 ; 
EQUIVAL ONCE (POCUMEG2EO1Y, NWR?) he Sc Pee \ 
“CQUTVALERCE (PCCLRCUSEAL I » I ) o i F174 129 
POUTVAL TCE €PrCurtciulds Het J oat we eet FLea 144 h 
FUMEEVALCHCE (PECURL C7251), SAVE) / r1?4 3} 
FOUIVALENCE APCCUR (CAEL), 62 ~ = F123 122 , 
LOUTVALENCE (PC CUR COATED)» 89 ) F122? 133 ! 
FAUTVAL INCE (POCURMCE T4129 SV ) | FL?3 134 
EQUIVAL INCE (PCCURECESID), Shtt oo} : 7 ~——™~<CS«~—~S Ps; l 
CYULVALENCE (PCCURCE7EALD» Set?) QQ FI2a 174 
FCUTVALEN OS (OCCUR LCTLEE) » ShlL3 bamenen F173 137 | 
FOULVALENCE (PCCUR(DERGLY, ShhL  } Cc} F123°«138 { 
EQUTVAL ENCE (PCCUPCEESCL) > She? } FI23 139 | 
LOULVAL FACE (PECUCCSEC4D DP, Shae? _) FI73 14" 
PG TVALENGE (PECURC 14 2214% Wik oF : ji FI2% 141 
CUEVAL ENCE (POEM CCEL4L Ds, TOR) Kn) FLO% 142 
FAN IVALEICE (OCCUR ECKL), WV) . : F123 143 ' 
COVEVALENCE (PCCURECHESD)» TCHLL D ee ae FI7% 144 ‘ 
CQUIVALENCE €PCCURCCEE2L), TOhL2 173 145 ; 
FIUIVALENCE (PECUREQETCLD» TOhLI ) a FI72 146 _ ' 
" FOUTVALENCE (PECURC SEIT, WHRL VY ia ae j Fi23 147- 
FOUTVAL ENCE (PCCUKCOB3ALY, TCHR? } Lam, F123 14A 
: beUIVALTNCE (POOLE EC O4ELDs IC4R3 a ; FI2% 149 
CNUIVAL INCE (PCCURCURSETDD, TT } ’ F173 159 
FQUIVALENCE (POCO CL 2201}, TS ) Os f12% 151 
TOIL ENGR €PCCUR ES IRITD, VS ) Ee igs ri?a 15? 
EUUEVALINCE (PECLRCCCGLL), TRLOT PO FL?3 153 
EQUIVALEHCE (POCCUTEOCSSLD, WeCPLTD Sn F124 154 
+ FOUIVALENCE (PCCUR(CC321)9 VPLOT } F123 195 | 
FQUIVAL ENCE (PCCURCGT 2929, WL19 2 , — FI2% 154 
LQUIVALENCE (PCCUR( C1441), bL22 ?) wa A F123 157 
POUAEVALENCE (PCCUR(KLECT)'» WLIP om eee PPA 
EQUIVALENCE (POCUPR(44481)_ KLIS)” = an “FL23 1507 
COULVALENCE (PCCUR(C4ES1L), WL2S J es we FI23 160 
“ FQUIVALENCE (PCCUP(C4tu1), WhL3S ] F123 161 
EQUIVALENCE (PECURTIGROL)s vPIP ? Fi2% 162 
EQUIVALENCE (PCCURLECSEL)» WRh2P 3 _ ~ "  "FY23 163 
EQUIVALENCE (PCCLA(MLI7L)» WRAP _ > end FI22 144 
SUTVAL ENCE (PCCURECS4CCL Ds Wels )) ee eee ELF 168 
SOUTVALENCE (PCCUPCEGLGL)» wh?2S  } F123 166 
EQUIVALENCE (PCCUR( 94221), Wh3S } - ct Van age F123 167 
: tPIVALFENCE (OCCLE(OCLEV) > ZPLOT ) ; F123 148 
COULVALENCE CPCCULECI2EL) > 2S ’ a F1I24 149 
FOUTVAL PACEETOCCURECTC3ICPLs, TY ) FYQZR f Te 
PLUIVALFUCE CIOCCURCCCIE1L), PREF OU Se ere toe Fo aS 
FOUTVALENCE CICCCURC ROR Ly, TOP) FL23 172 
kell TVAL CICEC INC CLR (ULC beds Sd ) Diy mean = F124 173 
E.IVALCNCE (INCCURL OC IULAb, K ) F127 174 
EVUYYALENCE (TOCCLPCIC2Z)« LPLOT ) r a FI24 175 
EOUTYALFNCE CIOCCURCGLAC3P, HONE} a ria 176 
POOTVAL ENCE CINCCLE CHES] NCOMDV) a dice "he ee £2 eg ca 
POUTVAL ECP ELOC OLE EMCO U 2) CUNY FL? 178 
POUVALENCE CPOCCLEL Way, HEP ) : F173 179 
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EQUIVALENCE (TOCCUR(CC365), NOECOY) 
CWUIVALENCE(TOCCURL CT ILG), NEVE) 
PIMPIVALENCE (NOCCCURE AC)» TPNALPH) 
PNGTWAL ENCE ONAICCUR( C7, LOPT ) 
FOUEVAL RICE (ROCTURE TE Ds MHEATD 
EOUTVALGRCE (NACCUREO IAD, AOPT 
PGEVALENCE (NECCURC ECS} ) NOSEOPD 
EUULVALENCE (NOCCURC LEDs NGEOM) 
FOULTVAL ENCE (NOCCURCA 2), NTHR ST 
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¥XX 3 6,00 

ry yt be 3 
€ CALLING THE “dimer SLEROUTENG TO RECORD THE STAPT CE. 
L CALCULATIONS ).BYi THE PRIGRAM. 

Loe CALL TIES ky UPR TES 
C Tad 


4s dn vnc * NUECRY a anak 
VeVC = MDWCH n I 


"oh yt 
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CALL SAVEDV( LeDVARYOCE LP ) 
IFCIREP ANE 62) GO INees? 
LE(NFECOY «FQ 6 Lt GUT 497 , 
Co LVL 1 ® de ND VOM 
NCGI- wIRCOMOVGCTIME ae a ne. eyes, ee Oe dye 
TRE VALUE OF TRE Tit clSopn VASTABLE TO Af PERTURBFD IS SAVEL IN 
pani 
FULT ly LT). ‘THIS VALUE eer) LSED TO CALCULATE BASIC CECCY TRAJELTLRYs 
PULDC 1) 1) = ACCUR CNC 


m9 


C CALCULATING THE TRAJECTORY, oF THE FIRST CUMPARISCN DECCY. 
ROL OC 2s Td = OVLET) fom 


€ FULEC3, 1b SAVES TH pen ThA peo VALLE OF TITS VARTARBLE TC BF USEC FUR 
C TRE TRAJECTORY OF THF SEGA COMPARISON ITCOY. 
LLL PGE OL 2) 1) = DVHC TD goer 


C 
C FOR UNE PERTURBATION ONLY s THIS VALLE MUST NE SAVED TK HECUlSs 1) 
‘ig Ih T Tr} PERTURBATIONS APE MANE, THE SMALLER VALUF MUST BE SAVEC IN 
HOC OC 2st) AND THE LARGER THR HOLO(3 eT). 
4st By PUP TE e Lehtive 
| TPAC IR EE RE 2) CO TO 498 
“ PRENDECOM= 294019 402944 


a UU 
WKIGE # } 
Gu TH 404 nie, 


t 
waite > 


 COUNE COMPARISON DECOY TRAJECTORY WILL BE CALCULATED, 
#2 wt = 2 
NeUGF # 2 | we er 
Gu TO 4uA4 : . 


Cott cave av tsi necny TRAVFC TORIES WELL OF CALCULATED. 


4% MLUH #9 
NETGr a 2 
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eee Oe es Ee teh eee - Sat perth 
" C SAVE HicCURC 192-1481, (OCCUR (Zc5=206) AND GéCUR {2221 IN DVAR( 1-19) 


: 
_C ROLEC2e 1) SAVES THE neevithen VALLE ME THIS VARIAALE TC BE USEW FOR F 
C 


§ he wer oa 
C A HASTE CECcoY TRAJECTORY KILL BE CALCLLATFD. \ 


7 Hm AHL mel I! Helter nee er erentesietiroremytep ttn 810 a unt pe pmeeesrmergent Spying] gegen 


FI23 1 


F124 


FI? 


EI? 
Fi? 
F.% 
FI?3 
F173 


ene 


FLo% 
F193 
FY?4 
FI?4 
FY\23 
Fl?3 
FY?3 
F124 
F122 
bk ie 
fir 
F1?3 
FY?3 
F174 
F1?% 
F124 
F123 


F122 


F123 
F123 


FY24 


Fi?3 
F173 
124 
F173 
FI? 
F123 


FL23 2 


FL?3 
F123 


“EL? 


F123 
F124 


FL a 2 


F123 
F123 
F123 
F124 
F123 
F123 
F122 
Fred 
ri23 
F123 
ri12a 
F123 
F124 
F123 
FL?4 
F124 
F124 
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ile ‘3 : 
Aaya sored es er aR not co MAE ABC vec pomp UME ELTA A Ml a eonlahutid eV AAMT eee TE TT i) 
’ 
c FIP4 aan 
1 Te GUROVCUND DS THE LOCATION UF THE OTS TSN VAQTAQLE BLING PERTUKLED. [193% 241 
4 6 NCON * NCCUMAYVEED) i F734 2a 
f . ee eee 
ri EONLETURE EVARCL§ 19) TO UCCUR CEI %=145) 9 DECUNE2 85-299) AND MLGIRE 222) F193 244 
CALL SAVENVE 24D VAR QCC UP} FL24 245 
CH Th 495 FI24 24h 
C ; FLO? 247 
CoIf IRTE TYUALS ONE A PENTKY VEMICIE TRAJECTORY WELL SF CALCULATEC. = FIT PAA 
Coat TkeF EQUALS THREE RERFRENCE VENTCLE TRAJSOTRRY DATA MUST ste INPUT. £122 249 
AQ NLGW & | se File3 225% 
“hIGR #4 F172 25} 
ge fC 444 KK @ NLOW MHIGH | F194 282 
IFOTR TENE?) Gu TO Sl F123 253 
c FIAA 254 
C PUT LOWER VALUE OF TDTH ME SIGN VARLASLE TO BE PFRTURRED IN LELUR ARKAY F122 255°” 
TFEKK fe 2) OCCUR INCOND = DVLOIED ; FI73 254 
C , P1273 257 
C PUT UPPFR VALUF OF TETH NE SIGN VARTASLE TO RFE PERTURREN TN CCCUR ARRAY'F123 25A 
Th( KK of0.3) OCCURINCON) = OVACLT) F122 759 
Sob CONTINUE F172 740 
IFCMUOC NF el) CASE = CASF & .CORDO Ce a a EM make 1 OM ip aaae co alone 
WRITECGs LLCOICASE DATE pE MO F122 262 
TFETSFR «TO .3) GO Th QF i : OE F123 263 
(e LPR 2H 
¢ TESTING THE INPLT QUANTITIES. FI?3 245 
IFILJPT.LE.49GG TO 53 F123 26 
WUITEC Es LOE FLOPT a Oh ice Sa aera a Sia ah Be to) ad 
LOPT a} F192 PAA 
: KA WANTESST ALES ZIGH TU 54 F124 240 
WE ITTEC Gy LOESINTHOST F129 27F 
\THEST #5 ; ; i F123 27 
b4 CINTINUL ERS Dad tee nara a net ee a ee ee 
IFINGFINM4 aL Te TIWRTTEC Ge LIAZPNGEQ¥ es eee F123 273, 
IFINCCOM eG. IIWRITE( Ge LIC 2) NGEOM FI2% 274 
A a TROINGECM OL Te LINGCOMAL 2 ee die Sse EP a aaa “| Inge 2) 
IE(NCEDM GT * FINGEOM =] FL72 276 
“TP CINALPRSEQ.OIGO TO $ ; ie Wek 8 ee Ree aor, 
IF CQ LUPTHCLOP T~2)) EQe Ob RI TF (691077) LOPT FI7% 272 
BR cut ue 7 DAW a gd Mag AN ghhe eat Et eel aN (0 OAR Osa 
IFCSOPT INE «OIG TN 5A FL24 pre. 
TFOMFEAT NEC SOIR ITE CG, LOSS) MFA YT Sal ds A HROEY ogee apes OS lee F124 7A” 
IFEROSSOP NE SCIWRIIC (6, 1C56)NCSENP ‘ F123 299 
55 IFIMREAT NE.OIGH IQ 5é acai as rea ais aaa 20 Ueto) ae 
LFINGSPOP NE SCIWRITE (CG, LCSTINGSCOP F173 284 : 
erat A Real ott ON ce meen tnos macnn mmm men meet ET gg ggg 
C F123 284 
C DEFINING NUMERICAL FACTORS. ee ES SS SR ea 
COML AM=7.64850= 9%MW/ 0 SIG* SIG) F123 2AR ; 
FACTE OSS RTE GAMMA } ag aie 5 te ee AG pe . 
AL AL far coe F124 206 H 
‘ ee eemennteoe ee near eee Bg Oa : 
C SETTING INPUT OPTIONS FOLAL TO A STORAGE QUAATITY IN URECER F172 292 
C TO CAQKY THEN OVER 10 THE FOLLCWING ENPUT CASE» UNLESS THEY ARLFI22 293 
C REDEFINED. : FI27 204 
LUPT aL opt co ee ee i ar + “F123 295 
MTR Tete Tue §T F123 296 
LOM Poe THALP HS | > I? ig: Yee aoe cae ie FI9a gay” 
\ F193 994 
t. CRANGING THE UNITS FROM PNCHES (fePuT UNITS) TC FEED (UNITS 317% 799 
| f,.... KNOT. REPRODUCIBL 
. 
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Peer n tern ht tee ee he be mem ng tae te erm wise . r = mere ey oe 1970-21 9 05000 9069 UR eR ep : 


‘ ‘iis MOAT A Use Wh cece 1 AAT A os 000 OA AA Me 
wiht tL ALU F TE CE 


| 
| 
t |! Ml 
i | 
1 t WELT INTERNALLY}. Fiz% 209 
i WN TERN TAL Ze te Fiza ™) 
RHLERBLALP GG en ree, eee ee a ashes et 
LAL AT L2.uud FY23 364 ‘ 
PA SEPA TAL? Cnt rI9s W4 
MUSEO 2712 Ue FI? 35 
LATHE ACS L2 2096 ; ; FI2% 3°46 
Ce L2etkl s/s i2etOe ‘ FY?4 347 
a DELYeOELYZLeew0O a Se ee ear: ae PIZD 7B 
¢ FY? 349 
¢ SETTING THE NKAG COFFFICIENT QUANTITIES TC ZERC. ; / FLFR ALY 
CORINE( dels PPO Py FY] 
CREINF( Zee LP FOLDE d FY?% 312 
CCEINF( ce le 2b tue DC FY?3 313 
COG IMEC pep NeueNG eo: he oe, ee Be ' vp dene. eee Tae 8 lee 
COMse of ' FIL>3 2315 
Cr] Copal mh Nebhey £123 36 
a GOP ate ef) F123 717 
pal Cle. sah eg eget tt ele ee ‘FYP2 31A 
CHL CLP sr NG FI24 219 
aa) CECORD ECR OD te CG eee Bee as ab ehede goyises RW ee okt eek nent oo EN GG a 
pees * CRCEIC#A COS F124 32] 
TC) TUTEST ae Oe 1st. oe FL>2 32? “ 
id FYF4 223 
a c SITING THE ROTATIONAL QLANTITIES TO 7ERO. FY Ae 
LFILGPT.SQ CGN TO 2C1 F123 226 
wea) Chey r0C SSPE CY Pe eR URNS Rr Kio 
Day Oe ea ea a F12% 327 
-) PS [hae 06 FY?3 399 
teas DO F123 9 ’ 
alee SAPS as eI FI?% 3270 
A PHP 8c eb F123 331 
PEEL OPT EO eLYTHEALC SO. COC Rew a i ee FIFI 242 
[x] ao. CONTINUE FY? 3 244 
LFRUL OPT ME IGO TO 2cs : a FI?A 344 
falar PO Au Jade 205 FL?4 335 od 
TAMAX( J P20 19 Se a) . , FLAY 336 
TMAM( I PSC AG F123 237 
ne POR ck ne vets co fo wae ae a is Ne el ee F123 338 
ALOAM ES eG ee. F123 339 
TAPINEJ 230.08 tet wee ; : pant ae a ee  FL24 34¢ 
Zy THINGS )80 D6 F123 341 
food FRINGES. - F123 342 
Bn Bue ALM EWP 8900 FL2% 344 
7 26% CUNT NUT i Se Oe ee ee ee 2 aU 
Gora wi FLP% 346 
¢ FINDING THC TRIGONMME TRIC FUNC TIONS CF THE THRUST CFFSET FLIP 46 
¢' ANGLES « FLA3 A347 
CTR EPCS 60 0G 74832 CH WE ZF) F129 248 
y eae of) F487290 Ce WI ZF 1) FI2% 7490 
PEL CT AGG GL TAS AAG CHP SI ZE T) te ut, PAM mR esas. bd ek ee ee F124 360 | 
SPELELOS NEG: CL T4E37 ORNS STTE 1) F124 491 
C F173 752 
C ChAGING THE UNITS OF ANGLES FROM DFGREES TQ RADIANS. | F124 353 
GAN But. ab LTAS 225 CHGAME 9 ; ‘ : FIP3 354 
SIG L405 WIGAMF A) ‘ F124 795 
PHY Sel LTS IPSN CHPHTC Pel evo aear., * i #2) i ag DFAS 9856S 
ONL a GL PSE A260 CHP SET . f FIP 357 
THEAL sete 174892 S108 THEALY ; T1292 4849 
WSTe PU LISE2ASOCAAL ST, Fip32 3590 
De me i ty ie 1 1 rig ” bath 
ATU OHTA TOA, ai 
Nhs a Nga Mie ee i 1] to itt 
i Be ah tat ag ' tue | Ah | ane " ah Mh ‘ 
i ” 
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CALLING [8 “SUPROLTINE WEXEN WHICH PERE CRMS THE CALGJLAT IONS. &172 


CALL VIXEN F223 
a F174 
PERTCOM THE CALCULATIONS CALLED FIR BY VIXEN. Beg ere 
CALLING IN THC TIMF® SLARDLTINE TO MFASURE THE TIME Tu FID2 
CALCULATE THE TRAJFC TORY. F124 
“CALL TIMER SCS SP VIM TAL) , ae ELPA 
F123 

» F123 


CALLING TN THE TIMER SLAROLTINE TO MEASURE THE TOTAL TIME TO. 
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igo TT PMN Po 


F243 


REAL IN THE INPUT QUANTITIES AND PERFORM THE CALCULATIONS. F193 
CALL TIMER St= ly TIMTAR) — fen Yate ae wee Fie 

3 F123 
SETTING THF THPLSY QUANTITIES TO ZERC. ee ee. 

 THEG UR : aon : ay “FL23 
MX eC Ont ry 
“uy ateGin os Oe ee eer ens. Pre eee eas 
M7 Bu oO F123 
TEATOLS0.F0G CU al ch ee i as Re eae Ae 7) Ur Je 
TEMATC Pd8C FOO F123 
‘se TOMAT( 2)20,°79¢ SEGRE edad af eae WE OM EN Sin Pan ie Ba eae as aS £134 
; F123 

CRKANGING UNITS UN THE GEOMETRIC QUANTITIES FRCW FEET WACK TU F123 
INCHES. ; 6 inn en Beat diolise =e ee 
“eyiseylet?.one Te SPT Re he et eae Seria F122 
RAEN 981912 .600 ; F123 
LALSLAL#12.006 “0 or aac F123 
WY 2euM 2H 12.C0C F123 
LAZSLAZHY2.000 vikins cat ees “FL22 
RuPsl A712 UDG F123 
COLY SVELY# LZ.¢NG 9 UT Sa ts ali a Soo ee aoe RL oe 
POL 2 eDEL ZH T2000 F pen eee eae & Fire 
FL?4 

CKANGING UNITS FROM RADIANS TO DEGRFES TC FREPARE FUR NEW FL" 
cee lh ean aa Ns NU dt Patni tes Ba ag 

AL ST #AL ST/8 20174632500 F232 

: : ae aval nies Wk ea 
SETTING THE CPTIONS FQLAL TO THE STORAGE TNPUT VALUES TU ry72 
PREPARE FOR A NEW CAS. F124 


AD 
FAL 
AD 
7 
764 
FHA 
BAL 
aH? 
CEN) 
249) 
270 
271 
372 
a7 
27h 
ams 
ATA 
q77 
a74 
279 
aan 
aqy 
a9 
393 
ANG 
705 
aAG 
any 
St FY 
239 
age 
391° 
"42 


mt 


' 
sermon toneoey Ha MB be eo ee nh Wa ent we A a a A TE Ee nu vet nner tA ym NE aL 
: F123 
StTTING THE MASS LOSS AND HEATING QUANTITIFS TC ZURO. F123 
QUT CLs Lbee ne Mes. a On eee Woke 
WaT ae LCI , rya4 
Cl J Jet eee a 
ANOT Led 20. cde FIP? 
MUOTC2e Jd 1206 ONE E129 
MCOUEC IS P20. le F234 
Burne rhe arg cee ete 
"PFPSUL J) e2080 | a i a F124 
QruTlersIper ne F123 
Por ays Pel. IC F123 
FJ?% 
wT ay arth = FL?% 
PACU Lae ene FL?9 
RBEUT EE ONC | a = “ia eka ae ae 72 
F122 
CALLING IN THE TIMER SLBPOLTINE TO RECORD THE CALLING IN UF F122 
SUSPRUTINE VIXEN « ‘F124 
CALL TIMP S€2¢TIM IAB) F?4 
ELAS 


393°" 


30% 
3957 
396 
3497 7 
48 
399° 
tale) 
4ny 0° 
aye 
a 
Ane 
AD® 
496 
“yt 
4h 
aie 7 
41 
4117 
412 
413 
414 
415 
416 re 
417 
41R 
41" 
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vo EK HH PAA lt 


jRHE yen fate ony Oe Cee cc OCG ea eee tt Hoa U0 Led ond MA MMH gd | Yat 
MOTO CeO CE Cerca 


IAL Pus INALP Pr PIDF 429A 
QTPM STEN THR TC FLD% aD] 
Ad Ce ag TS Pes Peers, - FIP% 42? 
ILIAC oe?) GU TH 3F PEO. 47a 
vy CALL HEADLEY FLSA 424 . 
¢ P17 425 
C UAVE MENTARY VORIEE PF TRAJEC INRY DATA FOR POST TERE CCOMPARISGA WITH THE FL79 424 
t ro Ag ECT Ny OF A DECOY. E123 427 
hii Me te de LRLUT ; i 
eset) « 2PL0TCN) F122 420 
‘ TSe1) = IOTET FL?3 4%, 
PROMOTE eh) GU IO 21 FI?% 494 
IC 1OP( G4) 6FOe1) VEEL) © VPLOTED) FI2% 43? 
IFC IOPCOS) EQeL) OSETD & VEGPLTUDD FI2% 443 
PRCRUPC CObeE DED SCL) = BF TAPL ID} ee ore se 
1 LEC TOO 47.6920) GO TO 22 rL?% 435 
PEC TOP COT) CQ.) WLESETD & WLIPES) FL2% 436 
IPC MUP COR EVe DD WLASET) © WL2901) FL73 437 
We bwPrCES PoE Deld WLASCE) = WEIR ELD F124 444 
CC TOPO EQelL) WELSUL) = HPLOEL) F123 439 
6 VECIRR ETT shed lars TA SOROS ne Se ek oe hes Seamer ee 
LFCLCPY 72) oOo) WeASKT) = WRAP CT) F123 441 
22 CUNT INU ; . E123 44> 
C F127 443 
WRITEC Ry LOCCDEZSEDD, TSELD gS VSET) pOSELD RSET) gf=1 pLFLCT) FL2% 444 
ERCP TRUS aE tee ATUL CEN AL SCE TS UU Mes OS 445 
AL WEASEL DWL2SED De WLISEID® LabsLPLOTY A PIS C123 446 
Lear F124 AGT 
TF(LP oh 166969 TO 47 F124 448 
¢ : FLOAT 4% ’ 
C INTECPATION ROUTINE HAS FAILED. RESTURC DESIGN CHARACTERISTIC VALUES FI24 45° 
CTU CCOUR APRAY THEN TRY NEXT RV GR NECNY. F1P4 452 
wh ITE Beeb OUP ee = Nee eh F198 482 
47 Cry TINUE F232 454 
CALL SAVEDVE 2e"VAR OCCLRI a ELIA 4Ge 
eae & CT URN F124 455 r 
2. CALL RRAQERENLOW) F127 46 
“¢ T1?3 457 
¢ a: THE DTFPERFNCES BETKEFN THE RFENTRY VEHICLF AND THE CECOY FL?7 45R 
* HUANTITICS © Ve VOT OVER Gy ETA» WAKE LENGTH AND RACAR FIP3 469 
C CENSS SCC TION FOX THE TNPLT RADAR FREQUENCIES. "F124 460 
ty a2 f= Ve LPLOT FIPI 461 
TECTOPE SI) GeO) GO TO 43t £122 442 
TFUMP Coad el Qel) OVITD = VEIT) -VPLOTIN) F173 463 
L MPELORD EB) O01) OOETD 2 OSL eWOGPLTED Pe 466 
TEMPE ESD EO 6T) OBCL) @ ASCT) -RETAPLEDD F122 465) 
ad PPC TOPC 74) 6EQ40) GU FO 33 ME OAL Congtta ee nl ae 
LELTOP COT EQe ET OWLLEL) = WLISCUID wWLEP(T) “ FL23 467 
’ MPEIOPC EM eEOel) DWL2CLD = UL2SCE) mWL2P CT) it . F123 4648 
PEC IGP CG The O 61) OFL ACE) = bLISKL) -hL3P(T) . F124 46° 
Po PRCEOP ETD E01) GHATOT) OMNALS IT) —WPTOTE es, EPEAT 
bai is tap FEC EOP CTL EQ eL) OWR2(T) & WP2STL) ~WR2P CT) "PLR ATL 
lens FFL TOPE TZ) oEIe1) OWRTUT) * WRASET) =WRBP UT) MPD aa ae F123 472 
34 CUNT TUE F124 473 
a FI24 474 
We ATECEs TOCLNEPPLUTCL) » TPLATEL) sVPLOT CL) sVOGPLT ELD RET AFL(I)y 18 Ly F127 475 
4 dueueT) F124 476 
i COOPER TOP (74) ED LIWPETE Feel OGL) CZPLOT EDD eTPLOCTOLD pWREEOLD o@BZPCLD,  FL2% 477 — 
Lat 2PCEDgHLIPCT) WLP CT) @hL3PEL) of al eLPLUT) F124 478 
Leoop FIIA 47 
| 
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thi ll 
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‘in fi a ] mi emetntephAt FC ree rs pL We Pe oD rt S  uaiaiiaalacaisaiallaenan aaa aRaeial ' dd QAP Meee wees ws ent Ae Yt ibe, pee ett ere a eae a ea) L 
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ua 


wordt 
SOQQIIOMn rapa ea end co mae pus ouE ate Et 8a TIME; mc ae TEA ea ols ae MTL ATV * 


* 
IF(ILPALT.60G0 Th 507 Fiza 4°) 
( ; FL?3 44) 
LOINTECPAT WN POCUTIUL KAS FAILED. RESTORE DFSTGN CHARACTERISTIC VALUES F123 487 
: CTO OCCUR APE AY THEN TRY NEXT @¥ U8 DECOY. £124 6A 
we They ane) LP ryLea ans 
CALL SAVCDWE 2¢f VAP NCCLRY FI27 4As 
WEITCO A, 4666) COVETI OVC DDR CED sti @LPLOT) PIPR GPA 
S666 FUKIATOIN @X7HPELTA VPXTHDFL TA DAXTHDELTA B/ELH pPxFLy.7 I) F12% 4Aat 
mb Ti FTF ae ae ee ee 
BLD CONTINUE aa F123 4A9 
¢ ‘ FiLA% 445 
Lo TRE SCWENCE UF STATELENTS ENOLENG AT 5) COMSTACS THE H AKD ZPLLT ARRAY F123 49) 
CotRE 2PLOT APRAY HAS THE ALTE S AT Wilt TRAUFCT&PRY GATA Ts UUTPUT. F174 499 
CoTR! OTPLUT AWRAY FAS THE CO28° SPONDING TIMES, FL27 494 
CoTRE & AMPAY HAS THE ALTETLOES AT bHICH CORREO re VALUES ARE LvPuT. F123 494 
GoPe foe THE TT ARRAY STHON THF TIME CORRESPONDING TC EACH th 1S FLUNL bY F123 495 77 
CLEP GAR INTCePCLATION ANO COMUENEE. wEIH THE TPLOT APRAY TC GIVE The Fi24 49% 
CoTTr ARRAY TE TIMES COPE SPONDING 1:1 THe TT ALTITUDE ARRAY. F12% 497 
Jy a 8 Ey) ‘FE 17% 498 
Ko 4, F122 499 
my $6 1 = dy LPLOT ee bw] F123 50 
JA = JH ghee ale =) FI? 59] 
C1 44 2 = UA, NCP F123 59? 
k= kK E 1 p brug F123 8°37 
IF(ZPLUTC LDCS 9) 38, 76537 > | F123. 57% 
a4 Jaz J 6 3 => F174 695 
TTCK) = HOD) FI22 6NG 
“CALL LINE ITUS¢ ey ZPLOTUL) p FOLOTEL) eWheXXt “e). . F123 617>°° 
uae WHC SIEXX FL2% GAA 
" Cut) 44 . : te | ELTA EN 
36 R= Jy 61 F12% 510 
TTUK} = ZPLOTOT) an Qo F173 511 
TIN OK DETPLOTES) ; Fiza 512 
fee era Fen, en en en eae ieee ees gig 
17 TTEK) = 2PLOTED) Ry FI?P% S14 
* TIM(KISTPLOTID) art . te? an pee CR Oa eG F 
co TA 5 Re FII2 B16 
44 CONTINUE Sen Rig AN ge a RO Ske ais “7  " py24 817 
St CONTINUE j Fi?? 518 
i ee <i eee ne rice 
95 Cl 223 JJ = LyLPLAT . FLPa Son 
SAVE Leda eo KK) & CVE) a a ae ine a as se ey ae 
SAVEL FeSdgKK) = OD( 00) ' FL23 922 
SAVCEL AeJJeKKE = OPC) — eS RE OS Oa aoe 
CSAVEL 4 SdeKKY = OWLICID lenny F129 524 
SAVECSpddeKK) = Uhl esd) ae ee a 2s eh ORE 
SAVEC he SJ eRK) -® OWL ASU) F124 526 
SAV IC Tedd eKK I = OWRD CIS) il cae A a ee, ae yl F1l?3 537 — 
SAVE CS edd eKKD & WR ZC : F122 578 
2 SAVECS JJ eKK) = OhR ICIS) ‘es > a FI73 520 
333 GanTINye : F124 530 
pm Pe ee ae ee ae ate a ances ceineees iat eee fgg aay 
C as ak F122 632 
NCUGAY = 1 eens ws ata T1723 533 
ITE KK EG NEIGH SAND ANDECOYSEQe 2) NNOCOY = ? C123 934 
IT 1K eB aN RIGH SAND NMECUY E02) NDICOY = 3 F173 935 
Cc. ; FI2a 4% 
C OISFLLR GAs PROGUCE AND GUIPLT FIS ST AND SECOND MTFERREACFS EETWEEN F122 537 7 
C PASTE GECCY DATA AND COMPARISON DFCOY DATA SAVES IN DC 332 Luce F123 $29 
CFT SPECTR ERD ALTITUDES AWD PROOLCE SCATICR® PLOTS NE THE SAVEO F122 539 
* 
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| 
tet Pr 
« i “y Al 
iy" le 
i 
if 
it 
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| | m ee inert nee Te Te TiC ner Se CT RTT oa ee ee ee 
' iq 
fi 
il iy 
| I 
_ | i) | 
, i | if 4 : 
; | 
i) C CLUPSTINCES VS THE VALLE OF Taf PPSIGN VARLAMLE KETAG VARICE F129 540 ; 
[ok C Tee © WELL BE THREE POINTS ON EACH GRAPH, F124 S41 \t 
Pid IFEYUPETM FUL) GA TO 6Le , a PUPAE ly. it 
| INCL (OAD EQS) F124 543 
ae LCALL UN VCELC Dy Te Dp 391 09557 ACV IVs TCV ESV GBH VELCCTYD F124 544 » | 
FE TLIC HST QD) FL24 545 
| LOM EN TGPL (2675 Sp 19) 6p 3984 iC oODe TCISD RH OF CE LPN ELPA AK 
IE io eed Oe) EL2% Sat r 
LCALL INTSPEC 2s 30 Te B5E 793555 DCH eM CHEST MH BALLSTC) ow, FLA? SHR if 
bud IPCLUP ETD ede) GO TO 464 F123 549 
- VEC UWI (GTAP O.T ‘ EYP a Hr. 
CALL INTGELU225 69 1G) 3841 e 3656s BOWLE FOWL! eTCWLE SWLIQSH WAKE LL) = F123 S51 
TFC LIP (B30 Fel) FIDZ2 AG 
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LOO FORMATO CLS e146 TFL 522} : INFEN 
ENG INFCn 
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1: 
CINTERP INTER} 
. SUMROUTING INTER PUMeNe he Xe Voz) ; Bad a ies INTFR 2 
IPL ICTY REALM AC Aw, O-7) INTER OF 
CIMENSTON X(M}_ YOM) INTER 4 
Cc INTEe 5 
c INTER 6 
G GIVEN THE Twn TAPLES MAND Y(X) AND A PARTICULAR VALUE OF X INTFO OF 
C CALLED We SUSROUTINE INTEPP FINDS 7 EQUAL TO YW) bY LINEAR INTFP & 
C INTERPOLATION, INTER 9 
¢ INTER LC 
c M 75 THE NIMENSTON GIVEN X AND ¥Y TN THE CALLING PRUGRAM. INTEP 21 
G X 1S THE INDEPENNE NT VAPTABLE AP KAY. IT CONTAINS N ELEMENTS INTFR 12 
C Y [S THE DEPCMOFNT VWARTAMLE 44RAY. FT COKTATNS N ELEMENTS. INTER 14 
G w oTS A GEYER VALUC OF xX FOF allCH THE CORRES PCANTNG VALUE QF ~ INTFR 14 © 
¢ ¥Y 15 TO BF FOUNO, INTEP 35 
c 2 0S THE VALUE OF ¥Y FOUND TN CORRESPOND WITH WwW. INTFF 1A 
€ ‘INTER LP 
c THIS VALUE TS FOUND 9¥ LINEAR INTFREOLATICN. INTER 18 
C INTER 19 
nay comm" ve M TS LESS THAN 2 M2 GREATER THAN My AK ERROR MESSAGE IS “7 INTE® 24 7 
C WRITTEN ANQO THE RUN [TS TPREI NATED. INTE® FY 
G INTFO 2A 
C Tee ELEMENTS OF THE X ARRAY SHOUL 8f FITHER MCNOTONICALLY INTER 24 
C INCREASING NP MONT TONTCALLY CFOREASING. INTER 24 
G INTFR 25 
Ce HM CHPFCK TS “ASE TNR SUCH CONSTSTENLY., 9 9 7” : pm INTIP 76- 
C INTER ?7 
¥ € IF W LS QMUTSIDE THE PANGE OF x, THE VALUF CF 7 WILL BE FOUND INTEC OR 
¢ AY EXTRAPOLATION USING THE FIRST TWO OR THE LAST TWO POINTS INTE 29 
G AS APPROPRIATE. NTTe 3a 
C INTFR 31 
ee eT EM AGT is a Ateneo Cer h Olle EO) SGe ae ieee ie eT INTER 32° 
WRITE(C ATO) Me N INTER 33 
: 10 FORMATOCING LIHDIMPNSTON © 1129 LX P2NTABLE St7f oe 112) INTER 34 
CALL FXIT INTFE 35 
EtG INTER 34 
¢ INTER 37 
any 28 LOCK YL GOT ee ee ee ONT ED eae 
ro 22 1 = 2, N INTER 39 
IF(W.GE AXCTII GO Tr 4A INTER 4f 
2? CONT INUE INTEP 41 
C INTTR 42 
a) foe N INTFR 43 
Br Bw CLV TROT HTD OXCT ISK CLE TID ECW KET ART) EV OT STD OUTS INT ER 46 
C INTFR 45 
RETURN : eat. S ; ; INTER 4 
C INTER 47 
iam 6 i ae : poe i. F ; "INTER 43 
59 0 54 f © 2e INTER 4 
meme TEC OLE OXOTID OO CIM 4” ASE a rae a re a ee re I ora COs at 
55 CONT INUE INTER 5} P 
cn ta 3G ‘ rs ots : eee INTFR 52° i 
Cc : INTER §4 
rm) a, ‘ INTFR 54 
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CINTGRL INTG® 
} SUMROUTINT INTGRL EKA SKS LX eKM KE oQoite SeSI eMC) INTOR 2 
IMPLICIT REAL#ACA-H,0-2) INTOR 4 ‘ 
| REAL % GPC, XO INTGR 4 
COMMON OCCUR RACCLR INTCR 5 
’ COMMUN /CPCCUR YS PCCUR CLLTTO INTGP & 
1; | COMMON /CICCUR YS INCCUR (329) INTGR 7 
a COMMON ¢TIMCOMS TT ; Bh ae a = ; INTGR 9 
CoeMON “/IRETSTY NOP COY INTGR 9 
OLMENSTON CCOURC ACI), NONCUP (30) INTGR 15 
DIMENSION ACAI BC4G) pC L4G) OD C1 G9) SE CLOM) oF (40) eG T491)HI40) INTO? YY 
| CIMPASTON LOPCONDSLZOLSD eQE4 Pd yp VELGO @SU4 9) STG L402 eT 209) INTGP 12 
| CIMENSTON 220200) 4 DOU?) INTGP J 
| CIMENS ION TRLOTOLEUD p TTA(Z OGD pXOC60) @ XKI209) YY (7990270160) 9 > OINTOR 14 * 
LOW IVALENCE (PFCURCLLATID® A ! INTOR 15 
FOUT VAL FACE (PCCUECSP4LT, H ) INTGP 16 
FRUTVALENCE €POCUTCLI759}— HOTF HINTSP 17 
EQUIVALENCE (PCCUREOS EAL), 7 ) INTGP 1A 
EOTVALENGE (POCCURCOCET LT), TRPLOT > INTGP 19 
™* PNITVALENCE CPORUPENCI61), 2 ! sae INTGP 2 
FULIIVALOENCE CINCCUeENenal), 108) INTG® PY 
TCUIVAL ENCEC TOR URES T3819), Jud INTGP 22 
POULVAL ENCE CIMC CURT AC3LIAD, KK ) INTA? 23 i 
FOUTVALENCE (CTOCCURCECIS7), LPLOT ) ‘ INTGR 74 hh 
FOUTVALFNCE CINCCURCR Nand), NCP OD INTGP 75 
FOUIVALENEE (PCCUR(O6961)» ft pe) Be re eS SE SINTGRe 26 
EQUIVALENCE (PCCURCOSCOL), ¢ ) INTGP 27 E 
COMI IVAL PACE CPCCUR(C reread, 2 ) INTO? +?A 7 yi 
EQUIVALENCE (PCCUPRCCS201), F ) INT OP 79 
FOUIVAL ENCE {PCCURL AS36116 F ) INTGR AA | 
FOUTVAL ENCE LOCCUR CES 4Pr), 6 ) INT OP FAL 4 
DATA LIA LI% 135) 1959 L2G LAT oT BB ol 39 gl4O ol GL lO? slLodelGay Lady INTGR 4? 
# 2 Se POE INTs PUN, AL 9g 222/ INT OP 33 
MATA XO/LOCHWY THE TAL PNY RAL LAMDAL LAL w2 INTGR 36 ~ i 
1 THETA? ®N2 RA? LAMOA2 LA? ZTURN = ZON ZOFF INTGP 35 
2TH) TOM TOF F rsp ERROR / INTAP A } 
¢ INTGR 37 7 
nee : INTGR 3A Mi 
C SUBROUTINE INTGPL PERERMS ONT Of FOUR MATN TASKS CEPENCING INTGe 39 
a UPON CERTAIN INPUT OPTIONS. THESE ARF (1) TN FVALUATE A INTGR 4° 
C CORRIDOR INTEGRAL AND PRINT A TARULAR ARRAY CF DATA Af EACH INTGP 41 
C OF A NUMBER OF ALTITUDES, (2) TO FVALUATE A SIGMA SWUAKE INTGR 42 
fe INTEGRAL WHICH [S$ USTO IN THE CALCULATICAS MF SUBROUTINE INTGP 43 
Cc EREFCT. {23° TO CALL THE PLOT ROUTINE TC ORTAIN A CCRRIDOR ~ INTGP 46 
C PLOT. (4) TO CALL THE PLOT ROUTINE TO CRTALN AN INFLUENCE INTG® 45 
C COEFFICIENT PLOT. INTGR 46 
c INTGR 47 
ae INTGR 48 | 
jd = CL XE203 INTGP 49 | 
ITC TOP CKAG 8K DZFO. 2) ON TN 1120 - ae “INTGR 5° 
KO = IMP(KAE GEKA) INTAR 51. 
IT(KO.FI.1) WRITECE, TCC! OC INTGR 5? | 
JG) 2 0 INTGR $3 
TEC TMP CKAL KOI,EQ,2) ICO = 40 INTER 54 | 
CALL LINFITEGOENCE ote Qe 800} INTGR 55 
CALL LIME TEGO NCP ably Reh g GP wT mmm ote oe INTG? 54 
CALL CMP TTC POP LPL Te 2p Dold i INTO? 457 | 
moe INTAP 5R i 
test zo : INT GR ao | 
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fA = %,09 INTSP 4f 
rao = 9,0n ISTGR 62 | 
* Ob So RPL IESE, OPS z a Pik INTGR 6? i 
RY = Oerye7z7Od) ROE OV) INTGP® 64 i] 
St se (URGE WE GCE) INTG® 64 | 
xxOby = Oy INTGP 46 | 
yy({ i) 2M} INTG® 6A 
Ze) &# St INTAP BT 
IPtP 1eST.ST) NX ‘2 ? ae no : ae % i INTGR AP 
IP (RL TOV) NX FC INTG? 40 
IF(KG.FO.IGN TO 4f INTGR 7° 
TECMY NE SD) LTE ST = I INTOP PY 
TECNX NEST) OP = FOL) TIT GP 77 
IFONK@= 1) 1G, Ay OF INTO 73 
Ly CT MRITE hp LLCO) ZL, TPLOTE1) «O15 RL Sle DAs ZULD ean ee TINTON 74 1 
cP TN 46 INTG? 75 
ee WhITEC Re LIAO) 2010 — TPLOT(L) eO1 5 Ply Sle OA, ZOU) INTOR 7A 
co om 49 ‘wTao 77 
3) WwOLTh( 6, VION) FA1) 9 TOLONTOLDs Glee Rhee Sty OA INT GP 7A 
. INTGR 79 
oe ie snail 1) aR 2) aS a ie Rie a +, INTSA RO d 
NP = NCP =1 aN eal | 
C ; INTGR R? 
C + INTGR 2 
ro YALL M = 2, KX INTGR AG | 
C1. 44 K = 2, NEP INTO" 85 | 
a)e See K mJ eg hae a hae Be tae sie ea As d 7 ay : INTR? RA 4 
‘ LECT IMD .GE HICK) ) GO TO 56 INTO? 87 
ah POT INV | INTG? at 
Joe MR . - INTGR ga 
50 rm B86 K & 2, LT PLOT NTS? 99 
et ee, INTG? 9] 
IFUT(M) GF eZ(K)} GO TATE 7 ‘ ea ONT AR ao- 
* HD CONT IMATE INTO 92 
toe hs INTGR 94 
6a eo? = BISeTIM) € Cf) INTG? 95 
are OCTeTOMNY FOOT) INTAR 96 
"S? 2 F(SOS TCM) & GES) INTO? 97 
ae Al aie) eer (peas BEN eae wane ate” ke eA ae ea ee Wee ger ae < ‘ ae” . INT GP ga 
YY{M) = RP INTGR 99 
77(N) = S$? ; ; INT OP 164 
NY = 4 INTGOR 19] 
TE{C 2.57.52?) NY & 7 INT 1%? 
TF(R zal T6492) NY = ] INT GP 104 
160 am f/x f N¥ © RE ee eee) wee a aS ake oie Sa <3 a ee aa ee i “INT GP Ws 
Cc INTGR ICS 
C : Te irhate Lr. 2 INTC? 1n6 
CO TOL LOD: 20%, ACM, SOM, SAC ABLE 970998099997) 1G0 INTGP 147 
" 106 T2 2 ALON: ; . : ak SS eine "= INT GE LAB 
Laas INT CF 1N@ 
“PA e DA EG AMDIL XT TL TOM eB CIT OL) DOLE SE OV) QAP INT BR LTO © 
CALL WE ITEM( SIGUE Ty TOMY e TIM (I) 9 92 0% 29527 sD ATI AOE TL VeKQ) INTE tL 
CO TO 98s ; ~ OINTGR LL? 
700 T) = T(M-1) : INTG? 113 
T? = m(CHIIRFOLIIZONESIH-90T1} “"" INTGR L14 
CA e DA G ANOLE Xe TIe TOM) ob ED) oC Ci) oe DCT) oF IED AD INTGS TLS 
CALL WOTTE EC IGTR 2s TEMD TTR OM) 9 02 R29 52eDA GT? ADE To 2eKO) © " OINTRR 116 
IFCTTESTCET.%) GU TA SEF ; INTHE LL 
C1 oT S45 INTGP 1A 
age T2 & CP (P GID IZON OE HF Od) : INTHE 1190 
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| TL = TéMeL) INTO 124 
| CA 2 MA GB ALOU XK Te TA eV CLD oF CED AP CU) 0G CID AD INTGP LAY 
CALL WOLTER CSCO E 2g TEND g TIM(MD 9.92 922 952 oA 9T2 DL 1 ¥4KQ) INTO? 172 = 
>. J TEC ITEST FQ ed) GO Th 366 INTOR 123 
fi ITrSt «= t _ INT HP 124 
; i 4 nia 1? INTGE 125 
1 Wh TD ew MICO ECTYIAON EAN) INTO? 1246 
a | fA = DA & ADDCLX,T? TOAD BCU) oC Cd) oDETD oF CLD AD INTG? 127 
fy CALL WTITEME IGOR 3 TOM) yp TIMED 9 Q? 982 9S? MATA ONE LE VeKGR OUT INTGR 128 © 
a cn Tn 969 INT 2129 
a nd TL ® T(M=1) INT EP LAM 
| T2 & (CKD IMEC TIIZAIRC IV -O01 DD INTGF 141 
i CA = DA G ADMILK TL) T2pB CI) oC Cd) DET) VE CLD DAD “ NTA 132 
i} CALL WRITEMLIGNG 204 TEM) p TIM OM) 992 9R2 952 97AGT2 OL LE VeKO) INTGP 133 
IFC ITESTLEN.O) TN TO S50 RIE ne ei estes INT GP AAS 
Gn TN $55 INT GP 125 
sou) FL = TUM=1) INT GR 136 
T2 = 7,99 MINT CP 137 
CALL WRITEMOEIGNE 1G) TEM) p TIM EM) 992 RZ 9S2e0AGT2~00 T 22K) INTGR 139 
cue TO 989 INT GR 130 
ANd Tl s TEMe1)° 7 Visage tes ee aI ; INTG?140 
1? = La earn tres ary INTGR141 
CA = DA EB ADDL Ke TL T2 sOCTP FULD SFC) GCI) eA} INTG® 14? 
CALL WRITEMCIGOG2L TEM) p TTMOMD 92 R29S2 oMAGT2 900 TE deKG) INTGP 143 
TFLITEST SCC.) GA IN SEF ~OINTAR 144 
cu TO 999 INT GP 145 
FON TR = CLIK FOTIVZIBC HOUT) gam ahs ENP Gr aae 
TL = TiM-1) INT? 147 ' 
PA = OA GB AMDELXy Tle T2eM CI) CCI) ETE SECT I eA) INTGP 448 
CALL W2ITEMEUIGUR 20 TEM) e TIM CM) 992 2 952 @DAeT2 gO EF deKQ) INT GP 149 
ITCITEST.£O.1) GU TN 780 INT GR 159 
ITFST = 1 INTGP.351 
a UTA Mi 7 oela i a) Sec ann sre mee s -| Y, |h aes 
150-1? = FIT) DCiau iit tame INTGR153 7 
TA = DA GB ADMELK, 122 TOM) DEL) of (1) oF OS) 6S CUD AD zs - INTGP 154 
CALL WRITE EIGIE 20) TUM) » TIM (MD 9 Q2 R2 952 eDAT2e00 T deKQ) INT GP 155 
. tn To 999 > INTHP 1546 
AOU «T= TIM=D) a INTGR 157 
j oe TP we CLE DIRGOIDIAID EDD AP EDD : eS cee oe ee INT GRR 
TA = DA & ADDILK, ee tah GALI LAGE CD of LADKAD INTGR 159 
CALL WEITRE (IGE 30, TEM) e TTMEM) (02 9R2 982 eDAeT? VDE T YeKQ) © “= INTGRION ° 
IFCLTEST.FQ.0) GN TO 65C INTGR 161 
cn 1m 995 3 INTGR 16? 
ayo Th THN) INTGR 163 
reese biggim © °c cues iene recone + ENTOP 166 
CA = MA & ABDUL Ky TL TMD AO E1) EEL) oP LAD 9 (dD oAD INTGP 165 
CALL WR TEAL IGDE Gs TEMD y TIMEM) 902 RZ 982 9D AST? DE 1 deka) : INT OR 166 
eu To 999 INT GR 167 
‘oso oIPEST es dt Eo Re s ae “INTGR 168 
re = 1? INTGR 169 
Wm 999 NXEEN SY ee rere ace ee eee ee) emis ne 4 . ee mee em te “INTOR Tn 
LLL CONTINUE INTGRITL 
C Me ihe a pats So ey Se AE 
¢ : INTGR 173 
NCCUR( KN) = DA : . . i vas: : ' INTGP 174 
NCCURCKNEKR) = OF INT GP 175 
TREK O sO oY WEIL TSOFD NE DAS DT TINT OR TTS 
¢ : INTGP 177 
112% TRC LOPE KAEKAY 69.0) GO Th pach INT GP 17K 
ct MT : INT GP 179 
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RL = CRE LIZSIGC 1} sa? 

fi 2222 T = Qe NP 

CALL INTERPC LEC PLOTSHUED 9 2eP eRHT) 
RHE = CPHT/ZSTCOLT SY pea? 

Su" os SUM = TRE B RETOCH ED) 
ALO = HT 

COMP INUF 

NCCUR (KE) = SUM 

FOIF = KEL) -HENCP) 

WlETACRe LACGT OCs SUM 
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IFC TUIPCKAD.EG. 0) GO TH a#ran 

CALL AVPLTEX De 12 eXXe VV e275? OC) 
TFINPECOY.LT.42)) PF TURK 
IPFCTOPCKAR SKA FO,L) CALL TNFCOF(I9C aT 2 4X9) 
PFONTFCOY LT. TE RF TOP 
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RETYON 
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INT GO LA 
INTOP LAY 
TATGE 19? 
IMT GO 1A 
INTGP IN 
INTG® 1a5 
INT GP 1M4 
INT GP [AF 
INT OR TRA 
INTGR [RO 
INT GP Lan 
IMT GE AD] | 
INT GR 199 
INTGP? 194 
INTOP 194 
INTGR 195 
INT GP 196 


INT GR LOR 


SHALTE TUNE LOX 44 TEMP LX LAMLEWER CORRIDOR LX SHRY=DE INT GP 140 
VX PBHUPPFR COPPING® LX YSHENTFGRL VALUFLX LRHLEAVE CORRIDOR 


TNT GP One | 
INT GR AN} 
TNT GP On 2 
INT GP 2A 4 
INT GR 24 
INT GR ONG 
INT GP 2n4 © 
INT GE 207 
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CANXOES JNXRE 4 
IMPLICTT REAL #8 CAH, f=P 4 one ' INK 9 

IMPLICIT REAL #8 {RP JNXRE q 
IMPLICIT COMPLEX * 16 ¢ 9 JNKAR 4 
SUPAOUTINE JINEBO SEZR oe? Ly NNo Cy FN@y PNT D JNMAT 5 
DEMTMNSTON SQEPPy SE(7) 4 TEMPEZ) JNKIE 4 
C JNKGE 7 
ae “OUNKAR On 
c guia itint JNXBES 1S FALLOO AY MESSFL Te CALCULATE THE BESSELUNKAr oO 
t FUNCTIONS (© MAKIMUM ORDER.» Np ANT AF R-T, JNK RE UA 
i, one LY 
OX= OCMPLK(IE ZT) JNK ar 1? 
7=COARSIOX) JNXSF 24 
pl a weiss ane ; mee 7 INK AP LG 
AXNWECE Due te MNO) JNKAE 15 
C JNYRE 14 
. IFt z= 15.49959999999 WU%slL etl NYE LF 
Ct JNKARE 18 
a UX # C100, C200) UNKRS Jo 
C PORN A 9.7.8 HANDOTIOK OF MATH FUNCTION S<1 5. DEPT. CCMMERCE "NK AR oon 
QCCe OCMPLAIC, eID ONKAE 9Y 
OXBe COC.) & OK JNK RE 99 
OPT e{39,141502653995792 4» O00 ) » JNK OL 24 
OXANSNFLOATONN J INYBE 24 
NACHT = OX @-f ECS MOMe ORNU Et.2500,9.07) Fe CPL JNK AR 35 
QMMY 2 (46.00) 0,0) & OXNU # OXNU : JNK AE 34 
Qr oe (1,099.00) NXE 27 

QO=(QXMU <C1.9C, 0.007) 7£ OxA . WNX RE A ; 
QZA =X 8 wm 9 : » JNK RE 20 
OPT = 9p JNY AE 
¢ JNXAE Ay 
' * Mf ‘12 Te?, “acre » 2° an SESUSSARASS 88 SAE TS SUN ETE SO 
QUX LEDELOATO Pe CT 17) JNKAE WR 

QPT = =APT & ( QXNU = QT ® OXEt * (ONKU ~€ SXE 612.00 e0.D0) ) UNMSE 4% " 
eer y 2 ANNE ¢ OFA JNXAE 5 
Apt = OP & OPT UNK RT Ab 
PTPN = CHAR SNOT /OPME W) JNXBE BT 
, TEE PTON @f oP 14s L416 oo ‘ ‘ ANKE BRO 
14 CONT INYE JNXAL 39 
QP = QPNEW bal vie. JNORE 46 
12 Q¥P ew XE 6f4e00- 6.991 JNKAE 41 
13 Cont uF JNXBT 42 
JNXBFE 44 
' OXT = C709 MOP me ; po em INE he 
NOT = An INXBE 45 
395 fet,4aete, 2 ace : JNXBE 86 
QXX TeMELOATI TE CI= 10) JNKAE 67 
: ANT # SNOT ® COMMU -QNT OMT 2 © COXMN = (OKE £02.00 0 eOG) PO*QINKRE 48 
1 pJ9RK ET A978 JAXBE 49 
QONPa = O06 NOT i 3 JINX 8A 
ROTON =» COARSE QNT NONE | ) JNENE 81 
IFC PATON = ©) Lae 1?e 1? JNX AE 89 
a? CONTINUE JNXRE 53 
vo 2 QONEW Got soe ate t oper en Lt oC) 
15) AKT = OKT G 1440090600) JNY BE 55 
t6 CONTINUE | SS Tg Sars aA aa | JNXAG GH 
¢ JNKAE ST 
ANS = CHSGATE (7.00 ,7.NO0/ OPE 7X) * COPNEH * (Cn CeSs CUCHT) Jayne 5A 
Lea *® COSINCOCIT EY JNK RE 49 
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CALL DRLIME OANS » TEMP J JAYNE 6 

fhe om TEMPE LT) JNKAE AY 

INE © TEMP COP JNXAF A? 

RETIN JNK AC 64 

¢ JNKAP Ad 

Lf CONT INGE JNXRE AS 

As},na JNY RE 6A 

mi LY taeda Nn JINX SE 4T 

LT ASSFLOATC DY eA a . 7 7 JAX AE AR 

PPRaIl ZF? Ne JNX AF AO 

727/708 UNO OP 

PaS ye T (7 eee PE? Ape Dy JNK PE PY 

THOT PATAMQET 2 21a? PR) JNXAE TP? 

Suve et Oy JNKAT 73 

SUE TRO, en ee Fr 3) . SUNK SP Pa 
ANG= SFLOATONN} THE T JM K AG! 

PAN ASR AMIN JA JNYRE Th 

XONK ECA SCANG DEANNA SNXOF PT? 

MPR TESINGANGPOR NAA JNY OE TA 
CWel.ie JAY A 

cs SOL LT) = 14965 ek ot are a Faas ies ‘s 7 “Jaen an 

SRO 2) = 1.065 JNA AY 

C AMM PI 

OY § Pal) P7006 JNXAP a 

IF (MOOG Ts 2925352 JSP Pim 

3 Mal JAKE AG 

: COT I mee tat nn nim me ce a ete Se TMX A AA 

@ a? JNX AP AT 

& Ke thy JNX AP oA 

ATS=OFLOATI Kd* THE JNKAF BO 

OWet VEVeRANEL OL TEL) 18 JOE LAAT ENNET) JNK AR or 

SROM DEO TVELOSTAT) JNKBE 9] 
ASTIN ORM SINAOM EE [A Tage Te E trees mire cea sen eesti eeprom eat pasa “INK AF 

S40 2 SMe 6 SAM ANXTE 9 

SUN TeSUM TA STOM) JNXAE 4 

IF CABSESR E20 SPE LINAC 15 59S JNK AF 5 

G& TF CANSE STOP SLOLEDH$CD 75595 ANKE he, 
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WRITEC Ee, tach) JNK AP 90 
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T EMR X ONE H SMR =X ONT SUM T ANKE 144 

IN T&X 22R* SUM TE KANT ® SUMR JANN ME YN? 

KETURN INEM TA 5 
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coven + SOMITE AMBRE CUCCURT OO LION» 2PRE 2. mn www REALS... 
evita ENCE KUCKUALOOLLONs 2PaA2 } aAbaAuleob 
‘wut WUCCURLOOLZL, 257 oD. ReAJL 4a 
erie Pine \OLENIALOOLABhs ZTURNI REAVI LOO 
; {VALE € (OCCURIOOLON), 20 =). REAUI boa 
eulva GNCE INUCCUN TOK R, LATMOS) KeAQL Lb 
vis a | iA Buu! | v Ail, ENCE tauccua t DF ) LJ IKCMW }.. ee emir ON 1 a ce emt AicAw [ los 
EQuIWAL ENCE {NUCLUAI30)¢ [NALPH) READE no% 
Ajo+ SEMUIMALENGE [NOCCURIZ9)», LTAPE } won vin, cer ie Read hos . 
EQUIVALENCE (MOCCURIZ3)» ITHRST) KtAVLLo® 
woe MUM ALENCE, [NOCCURIOT)» Ler 4 Reavlio? 
EQUIVALENCE (NUCCUMC20)— MATLNG) ReaAvl log 
ma EDU VAL ENCE INOCCURU23 } yA FLA23 RESTA = — 
conehbey yt (NOCCURIOG)> MAXVAL) RKkAUIATU 
om EuUL WA ENCE inde. WaRtlo #. PA TRGRN TI cstctgars x tonctommntiricn vote naa se ReAULATL 
EQUIVALENCE (NUCCURELO}, MXTAGL? REAJLATS 
Cees « EWUIVALENCE (NUCCURILT}s MATABZ? ee ee REAJLASS 
EUUIVALENCE CWUCCUN(LBIs MAMCOD I HeAUL Le 
Hel) ami i AE QU WAL ENCE... CwGCUR 19) i] AAXACO) Stok > le eden IND Le oben mn td ere REaAul 1? 5 
EQUIVALENCE (NUCCURTAS)s NGEOM } KtaAvlizs 
iH EWULVALENCE. CNCCCUR(O3)» AUPT  ) HeAUEAT? 
FUU LV ALENCE Jno apap MUSEDP) HEAUL AT 
EUUIVALENCE ANUCCIKE 24). NPLUT J KEAVILTY 
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EQUIVALENCE (NUCCUR(L 4), NPRINT) REAVILEO 
EQUIVALENCE -INOCCURI22)—9 NIHRST). ~-——-——— | READ DE. 
EQUIVALENCE (1OCCURCOOCOL), [oP } REAULLBe 
“ c a Vat; epee zi ascend kes : = - REAULLaS | 
¢ READLIO4 
~DATA INCOLT/838a/- - - ~ KeAvllus 
C READL abo 
By Sh vdsaemtonseeten econ «|e salina Coie tah pseat ban Noone te Vosncatis ais Seooheasts aa"? Papeete wi ioc ccinnjenil aap arp eer oes cesbieas: Jia ac ALE A DLE hes 
e dcavl lob 
GU TO (let) e ILL ill —.- Sinke (sine Nae Reaulaby 
L CONTINUE Reavll yd 
CALL ZREADX . - AcAvI LL 
ca is SETUPL BHICCUR rly UCCUR, 4000) Aeaulive 
a sa eee Berghe age wows KEAULLYS... 
G “CASE LDENTIFICATICN INPUTS KEAUIIY4 
CALL SETUPIBHCASE eGeCASE J ReAuvligo 
CALL SETUPLYHUATE 96eDATE } KLAUELYO 
CALL SETUP( BHMC MU + GsEMO i. ReAulid? . 
C KEADLL94 
CG asad SP ACENEENG: BAO: STOR CONTA Sis sa oc emecictcu ese Sen tiyaies ts fa case aa KEAJTL9Y..... 
CALL SETUPEGHNPRINT e4sNPREINT) REAuL edu 
CALL SETUP OHEPR, eGelPRa }. READLEO’ 
CALL SETUPE SHZUAK 2GeZHAR ) Reaul ede 
CALL SETUP{SHZPR2  .¢GedPaz } * REAUL EOS 
CALL SETUPLAHZST 2GeZ5T } REALL £04 
- -CALL. SETUPCONTST 2. p Be TST ) = Ee A Sa cece ae ee ROAULENS _.. 
CALL SETUPIGHALST eSsALST } REAVLZOR 
C ee ie REAITZOT 
‘ ¢ INTEGRATIUN ACCURACY CONTROLS | * READLZON 
CALL SETUP(SHDELIN 9 8,0ELIN } " READLZOY 
CALL SETUPCHHCHIGH  y8rCHIGH 16) ReAvl cio 
CALL SETMPCGHCOUWN © 9 &eCOUWN elo) ----—- 2 ee ee, keAJleds... 
¢ REAVI did 
¢ INPUT TRANSITION ALTITUDE TRZTR REAVIZELS 
m CALL SETUPLAHTKZTH 2eOeTRLTR } ReadlZzl+ 
G . wee i yee e Reaves 
C ‘ INPUT LIMITS OWN THE FAIRING REGIONS REAUI elo 
som GALL SETIPOBHKLUP. 9 @eXLUP. .} mai a1 tar etenietimatninietn «cement fee a adem wine RCA SEE no 
CALL SETUP @HXLLOW elyXituw ) REAuI cas 
CALL SETUPCAHXUP + &eXUP Wf oe ks , KEAvI EL9 
er SET JP( dHXLOW rSeXLOW } READI 220 
c ” Resvleds . 
C “TRAJECTORY OPTION AND INITIAL VALUES. OF TRAJECTORY / PARAMETERSREADIZE2 
> cow GALL SETUPEBHLUPT op MLR To 2. ne ee meme ee ene nme MEAN SES. 
CALL SETUP 8HZ0 14520 ) REavl2Ze4 
oo GALL SETUPE BHGAMFO.---»8,GAMFO.).-- - 2 Lk... mae . REAULYSS 
CALL SETUPI sHVO 66,V0 F] KtAvl 226 
CALL SETUPCSHTO .. ~ ..98sT0 | Se ee REAVL ce? 
CALL SETUPU'BHXRO e6eXRO RtAvizds 
mee come GALL + SE URC BHP Oem 89 PO eee = AY nen eee nee see eerie: ater enene REAUL ZED... 
CALL SETUP SHQO 1 8,Q0 ) ReaAuiesO 
sve GHLL SETURL GHSMMO————— 9 By SMRO 0) ne ee ees oe REAJL ZS) ..... 
CALL SETUPLBHPHIO eGePHIO ) KEAVL 232 
- GALL SETUPIBHPS TO edePSIO . ) erates ReAvléds 
CALL S&TUPLBHTHEALO 4 8sTHEALO) ReAvles* 
C ae fe tO veoneens. . READLZS9 .... 
c INPUT TEST QUANTITIES FOR ANGLE OF ATTACK CYCLE. TIME. REAUI eso 
CALL SETUPCSHICRIT eUeTCRIT 3} Reasles? 
CALL SETUP(BHTECUN ¢@,TECON ) ; Reaul2ss 
f Readles¥ 
¢ 
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t NUMOER OF TABULAR INPUTS FUK EITHER THE hela pened TUNNEL REAUDL 2460 
Coon . OPTEON OR THE INPUT-TRAJECTURY.UPTIUNe. .--~- .- ReaAglebh.. — 
CALL SETUPCQHMAXVAL 94 MAXVAL) KtAvl dts 
ee Te. a. REAVI 243 ‘ 
c INPUTS FUR WIND TUNNEL UPTION KEAI1 644 
CALL SETUPCBHWIZ eBeWTZe75) KeAUl ees 
CALL SETUPLOHWTMINFE  ¢8)WTMINF, 75) REALL d40 
oweCALL SETUPCBHATRINE....9 QpHTRINF 6 75} -——-—-..-- ee ee ee READ ET... 
CALL SETUPCSHWTPTOT 9 SswTPTUT) 75) REAvL 246 
C gags ae Awad REALI E49 
ct INPUTS FOR TAdULAR INPUT TRAJECTORY REAVLZ50 
». GALL SETUPLBHTRAST «= 9 8sTRAIT, 75) Set eo nee ReAvl 25. 
CALL SETUP(BHTRAJZ 58 9TRAJSZ»75) REAVIZSZ 
sear meae CALL ~SETUPASHTRAJV —~» Bp TRAIV, 75)... REA 25 
GALL SETUPEBHTRAJW = 9 8p TRA 75) READl 254 
CALL SETUPCBHTRAJRN =» 8eTRAJAN, 75) RLAIL 255 
CALL SETUPCBHTRJALP 2 BsTRJALP275) REAUI 256 
Cc : es: -. REAM S97 
c TABULAR INPUT ANGLE OF ATTACK FOR USE WITH LOPT#1 TRAJECTORY READI<58 
coe ona CALLSETUPE BHINAL PH ..p 4eINALPH).. ee ee ee READL 2592 l 
CALL SETUPEBHALPTAB 4 8sALPTAS.75) ReAUE 260 | i 
c aa REAJIZob | 
C MASS LOSS OPTION INPUTS Réeavlebg ai 
CALL SETUPI BHMOPT e4eM0PT ) Revd 23 
CALL SETUP( BHMHEAT e4sMHEAT ) REAUl 2b4 
oa samme or GALL SETUP ( GHOSE OP..-¢ @sNOSEDP). ..- 5 Oe eee et eee es READ Lee ols, 
C REAVI 266 
‘ TWST IS A WALL TEMPERATURE USED.IN RAREFIED FLOW CALCULATIONSREAJLZo7 
CALL SETUPL SHTWST sG9TWST } Reavdl zob * 
C RkeaAvlcéo?. 
C TINLT IS INTERNAL TEMPERATURE OF VEHICLE USED TO ‘git al REAVLETO 
C . ws TEMPERATURE GRADIENT .IN. MASS LOSS CALCULATIONS USING THE . Revd 7a... 
¢ ITERATIVE METHOD. KEAUI 274 
GALL SETUPCSHTINIT . efsTINIT.) : KeAule?3 
¢ Keavl 274 a 
1g UPTIUN INOICATING THE SET OF GEOMETRIC VARIABLES BEING INPUT Rtaule7o 
1 GALL SETUSLBHNGEOM 94 NGEUM } REAL cTo 
“¢ UPTICN TO INPUT: CRY READ cTT ... 
i CALL SETUPLOHIKCMG  s4eIKCHQ ) READ 27d 
C GED ETRY JNPUTS FOR THE FIRST. CONFIGURATION ReAvLet? 
Call SETUPI SHW) OW) } ReAvl 280 
“pote GALL SETUPCSHTHETAL 9 Qe THETAL) REALL Zba i 
CALL SETUP( OHANL sGeRNL =} RkAvI eve iW 
min mieer Ae SETUP C BARB L .———g Bg RBM...) aa ea ae i ee MBA 28S if 
CALL SETUPCBHLAMDAL + 8yLAMDAL) REAUI c64 if 
reat SETUP(SHLAL eBeLAL } KeAUI dS 7 
CALL SETUP(BHMXTAGL ¢4eMXTABL) LcAVI dob 1! 
witents pit SETUP(BHCMGING + 8sCHQINL) : ReAvl cb? li 
CALL SETUP BHT ABZL rhyTABZ1 «50) REAQL 266 lf 
einen GALL SETURE HTEACGOL 5. Oe-TXCG01.50) enone ah EAVE SBF. 
CALL SETUPCSHTABIL +B TABLE +50) KeEAdl 290 
1 epi Af fb SETUPC BHTABI XL—..0BeTASIND BO ee OR Mo eae REAOL 2y1 
Chile AEADIE9z 
| oe lat oe CONFIGURATION HEATSHIELD PROPERTIES...:. —... ... KeAvl c93 
CALL SETUPCBHMATLNL +4 MATLYL) AcAUL 294 
seme coe we Ah SETURL BHTW Ys eng Bp WL. Dooce am te ee ee ee MEADE EIS. 
CALL SETUP( BHBETALL + BsBETALL) KbAUL 296 
CALL. SETUPCBHUETAZL. e8eGtTA2]) © - . ..- Fe Abaulc#? 
CALL SETUP(BHUETASL 4 8yBETA3L) 4 deaul 298 
(CALL SETUP(GHUETA4. 1B eHETAGL) ; Kuauley? 
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a caress : cavutetesdanecndbay tnt inl ALIMENNAL HAUL UY ] 
CALL SETUPI BHHREF) rysHREFL ) READL 300 
eon creel oy SETUP UBT ac ites! Mg FMM wea Psst sche araai a ane REAUL IO... 
CALL SETUP{BHRHU21 eHsRHJ2) ) REALL 902 
Bee screen oe arene tr a eee te - ME AU SOG 
CALL SETUPC SHEP 2) O,CP21 REA S04 
CALL: SETJPISHCPGL. -—- pd lPGL ‘aes hed DS 
CALL SETUPEBHNSLL styNSLI ) REAUL SOG 
woe te eee EAL Le SE TUR CHING Tb omen m ge i NS TL De tine cane oe arc AEA BLP sen 
CALL SETUPLBHNGLA rp&ysNGLL ) Rew SU 
CALE SETUPIBHNGT ~: rByNuTL 5 Rieu oo4 
CALL SETUPEBHDELHCL + 4eVELHCL) Ke mod SLO 
CALL SETUPCBHEPSILL -deEPSILa} REAL Sho 
c Real dad 
Ge -2TURN TS ALTITUOE-FUR DISCONTINUOUS CONFLSURATION C4ANGE RE AUT IL A. coo 
CALL SETUP{BHZTURN eHe2TVRN ) Ae dk om 
c Akaul oa & 
Cc GEOMETRY INPUTS FOR THE SECOND CONFIGUKATION Kekul dio 
- -CALL SETUPC8HWE sOeWe REAULOLT 
CALL SETUPCBHTHETAL 5 &sTHETAZ) KeAvloab 
wine ee GALL SETUIPCGHRNG. «9 GpHN2 2 Qe - “ “ HeAwl ole... 
CALL SETUP{ BHRB2 18skB2 ) ReAwl oe0 
- CALL SETUPIBHLAMDAZ.- 5 asl AMDAZ) REAdloda 
CALL SETUPLBHLAZ eOelAd } ReAdloed 
- CALL SETUPLOHMXTABZL 94 MATAB2)- KeAIL D2 » 
CALL SETUP(SHCMQINZ 5 &eCMQINZ) REAL se 
--- GALL -SETUPCOHTABL2 sbyTAUZ2 450). wht nena acen on +i mn REA OZ 5 
CALL SETUPCOHTXCGD2 5 GeTXCG02+50) HEAUL oe & 
CALL SETUPI BHT ADT 2 98eTABI2 +50} READE S27 
CALL SETUPCBNTABIX2 5 EpTABIX2¢50) KEAVI S24 
c KeADl sed 
¢ SECON) CUNFIGURATIUN HEATSHIELO Lidl KeAdl say 
CALL SETUP(OHMATLN2- »4eMATLNZ) ‘ ” ReAdvl doh uw 
CALL SETUPC SHTWe oeTw2 ] KeEAUl 392 
Call SETUPIBHSLTAL2 9 SebETAL2Z) REAQL Jos 
CALL SETUPLOHHETAZ2 »HetETAZ2} HeaAgt d34% 
CALL SETUP(ANHETA32 2 &eHETA32} Acadd odo 
CALL SETUPIBHUETAS2 » &e6ETAS2) WEA 356 
és » CALL SETUPLBHHREF2 oo pBeHR GF 2 bi ce ee einen mminnennnes name Head s3T.,... 
CALL SETUP(SHF2 9 BeF2 1 Heavl 3a6 
“CALL SETUPLBHRHUE2: «. p> BeRHOZ2 .)..- - ~ REAQL S39 
CALL SETUPLBHUELRH2 »¢ E,OELRHZ) KEAUL S40 
«CALL SETUPIBHCP2Z2 . ¢BeCR22 Fo. - HEAUL S44 , 
CALL SETUPI EHCP G2 eHeCPG2 = } KeAgL dee 
wenn GALL: SETUPLUHNS LS «9 WON SL dS -deniee cc  ne ne e “ KREAULISGS .... 
CALL SETUPISHNST2 eBeNST2 3 REAVI 344 
«GALL SETUPLOHAGL 2... -—» GsNbL2. J~~ =e REAUDL S43 
Call SETUPIBANGE2 eBeNGT2 3 READ Soo 
~ CALL SETUPEBHUELHC2 —— @sDEL HC 2) ee ck emi meee REAdl 347 
CALL SETUPLHHEPSIL2 + @eEPSIL2) KeAvl 34a 
‘Crees Ns ties Oo Sis ods ceoeceie ee minis Catt teeeas HEAIL S49. 
C DRAG COEFFICIENT INPUTS” KEAUL 350 
GALL SETUPLBHAAKCD «94 pMAXCD Bo cee eee a AEADI SSS 
CALL SETUPCSHHT AG ebsHTAS 4575) REAUI SS é 
CALL SETUPCONLUTAB . »8,CUTAB. » 75). ite REAUL S20 
Cc KEAUL 354 
C INCREMENT IN ORAG COEFF... INPUTS scent cil KEADL 55... 
HALL SETUPCAHMAXWCO 3 »94sMAXACO) itaudl 4>o 
CALL SETUP UHWHTAB eUsWHT Ad »75) Keaul aot 
CALL SETJPOISHWCUTAG =» BpWlUTAG, 75) z REAvI 336 
CALL SCTUPC OAWKEF ebyAWREF J ALAUI S99 
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C READL 360 
G we ec PAIRS TRING OP TCO. PROUT Son ene mecca ceeeseenn gonna - - KEADIBSL ... 
CALL SETUP(BHNTHRUST s44NTHASTI REAUE she 
CALL SETUPIBHTON a4) TON Ce ee rr REAUL SOS w 
CALL SETUPEBHTOFF eSy TUFF} ReAWL Som 
CALL SETUPT ANZON 0 9e20M ee om Réeadl dab 
CALL SETUP( BHLUFF e20FF 1 READL So6 
vaso Ee ills oh ETP THR ST a THR TD, oe ce neccenectent nmmneteamtnnn tt mmm weni epee sameeren HEAL Po « 
CALL SETUPIBHTHIHO = yaeTHTHO 425) NEmwd doe 
CALL SETUPEOHTHOELT ste TEL Te 25) Kool ao 
CALL SETUPUBHTHOELL se TMeL 2a 25) KeaAdl 370 
Cake SSTUP( Sato 0 By THO } ReAUl ala 
CALL SETUPIBHISe vbr SSP } Heal até 
» CALL SETUP BHA Oty ere eee eee Cee ni ches Kea sta 
CALL SETUPCDHPSIZET «ber STZET) Kéeaul sim 
CALL SETUPLOHTMEZET sy SeTHEZET) ReEAUL STS 
Cab SETUP AHOELY e6pDELY 3) Keaul Sita 
CALL SETUPLBHUELS pe ADELE «de o. REAULaTT 
t REAU! atu 
i Revie voson i sw PMT. FOR. AGUIAR NP UT. TM OS PERE enemies en en veneee RAUL TW 
CALL SETUPUBHIATHOS «4 )1ATHOS) Reaus add 
CALA SETUPUBHUPBNDE . ¢8,UPaN02) KeEAIL aba. 
CALL SETUPCGIONDNDZ 98) CNBNDL) , READL soe 
CALL SETUPCAHTBATML ¢6)THATMZ 50} * ReAulads 
CALG SETUFLOHTABRHO = ¢8+TABRHI 50) WeAdl Sua 
vain enim lll SET PC BHAT ABS ND 6 Ba TABSIND SOD eens teem ol AAD, BB 
¢ READ Ido 
¢ PLOTTING AND TAPE INPUTS KEAWI S07 
CALL SETUPLGHNPLUT «e4eNPLOT 45). " Heavl dub 4 
CALL SETUPCBHITAPE Ga ITAPE PO Ww KEAUL 309 
t REMWL 390 
aie ee. ig) MS REAL |CUNSFAMTE A, «61 n cennninn wamneienenmiiss uci sheet an bin mn oo eos REAM SD 
CALL SETUPL BHC o ByC } READL oDe 
CALL SETURL arta 289k ’ KEADI 393 
Sted SETUP ANGAMMA +OyGAMMA } KEAUI B24 . 
otk CaKued atte +e GRE } REAUL S95 
i vil SE art ph aHCARG +8rCAPG HEADLI¥6 
; aL TOE DIA, |b. se nnn MRE wey lf ose attests ef sme tape nen AEAVL SOT... 
CALL seturt BHLETA = ¢ BrLETA } Heal 340 
Cat SETUPCUHSIG ~~ eBeSIG » Heasl ary 
elena rUPt eH efstw Kt Avl 400 
ia (pl er yee ee ee tase A iinet wt mite ww EATER. 
G CURVE eit CONSTANTS HEADI GIDE 
(yal a He ie GALL.SE TOR dH A. Pminnfohbitetmmstt whl SLOd. Lal seadlanietiaeal eemeamadianaliahandentpmneimmmanmnnaliomammamaiaiaaaninandents ventas niateass vee AE AL nae 
CALL. SETUPY Od eBrbsdid READIED 
et ‘ KeaAviaos 
2 CALL KEAUEN( INCOL14 $99) REAULSOO 
ip oa 4 OPUTS )oa = 0, ! eT KEAQI 4 IT 
DO oo t= by 3 ReAUledy 
oti i Lied ee ere ee Hotta de osha seen eeausatinooenemsintsaeec Ne RA A AO hac w | 
10 m0 Hebdlesd | 
ee "a ge Kat. Pe) per sl aerhabcnielibseaiiroeatiearin inal oats slaeteseoriaiees men Atri meih ee ori aang vw i) def ae @ READ 4L A 1 | 
ersOR(AD East) 1OPid) = ih ReAulabe || 
98-4 OMT TINUE Le a a ee ee AcAII4LS | 
LFCTOPC SH eGo) 10P073) = 1 KeAUIOL® | 
blll tl) Gain ints TN tt + to LM? eo ha ee Ce ee ee a ne MY th 1 Nt REAVI4L5 iat ' 
rord iy = 0 ° KEAvl 4.6 | 
ih nem a n) fl M4, Ae hii ef mts ree ee oa! Acadl oat | 
ni Ba ti» Dah 27 : Reaglass 
de bee . REAILeL® id 
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th DEC POPS GRO) Babe Mb LP TEM te Be ie ee etcetera ee oe mentee ce cw amma ABAD 2 doo 
WOT, Ke Tp Ob eb Heb bed 
* ee ae ee Ce ; Riau edd 
TT CUNT TINUE Ake dt 
~ TROTTEMs Ee db -POPCTG) we be- He wee . ae Mh add oid 


PECAOPC I) Ee) eh et a ee head ene 
BAe CUNT ENUB em haw: nesiettietieee tes wie et meee ons vn eens MA Te | 
CALL HEDING(6) Hebot oct 


GO Ti 98 on Abdul oe ¥ 
94 EFC IAPLUTEQ.LE CALL AVPLT (HH eHeneH edo l} WHeaul oad 
[FIL TAPE.NESGIEND FILE 6 Ah AL ood 


VFI TAPC. Ms ‘adios a Hee eo 
te «CALL GKET. ee pictaveteeag Mansi Tone Mito, Meco stun tte beg ear areas ap Picaieys ake S OREM COE EO oe 
98 RETUAN Kel ead 
END . Real ao 
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eh I ae eer een oat sad Dearie nts onan a ow eevee WEAGY © Be on 
SUMROUTINE “READY Realy 2 
THOR LCLT REAL SE [AeHe 2) . ReAsY 3 = 
COMMONS BLKOS E40 e409 @XE401 GL 4D) S040) KPL40) GPL GO) TEGO) eGBCGOIKEA DY % 
daeSe USPe GTPeGSS, GIT. GSHeF eFPePOeFOLEs PeTOyRS Ske ZeeArEle KeAuvY » 
POELTACFAC 9C (40,10) HEALY 6 
mot GOUMMOMS BLES Me Ny loed Sa Ml eMS a NS oD Te Rel nec ce cet te eee cn meee oe ERY oe To 
COMMON /NALTFGZF SSA NALT AbauyY 4 
CUMMUN /ENO/ ITERM Aekuy oY 
t REAOY HeAUY LD 
¢ ReGUY ba . 
t aka aY 1S 
"9 SUBROUTINE READY. PERFORMS—THE. CALCULATIONS.OF. THE DAVIOON.... REAWY Lo... 
« MENDIMIZATION MET HOO WHICH ESTABLISH A DIRECTION ALONG WHIG4 [TO AcAur 44% 
th SEARCH FOR A RELATIVE MINITSUM ANDO 69% OFF AN INTERVAL IN KeauyY 15 
t THIS DIRECTION WITHIN WHICH A RELATIVE MINIMUM IS LOCATED. HEAUY Ld 
G snee s REAUY 17 
C RehaY Lb 
eee eee Da Semeeta tte : we ee ee MEY OD 
Go TO (200+ 201) 9b Reauy 2d 
200 [ Ta) READY «1. 
201 CALL MATMP (NeNeHe iy S 3 READY da 
7 - ws . ReAgY do 
#02 OO 203 Iwi nh REAWY 24 
PPE Mra Sy CCL Blom wa AM DED «ic 8esiee werecommmnrstsd Siceseeecmeeecentgceacs emcees ity toate ip dha ctonns ienrsas tolls ats ee wae w LU AIRS RR 
t Reddy 26 
204 Aa] hoe : Ready 27 
TOS CAGL MATHP (MaNeSeGeGS I Reavy ab ‘ 
200 IF (GShENEOT 227s 227 KEMUY ¢¥ 
G RKEAQY 3D 
BT PD SRE FIGS Dm ee nettenmeenennteniene: nem cnet Hott fom adamant Om me gee eae ~ READY Sa... 
WAITECGs LOCO} GS, Ey TPL, IT KEADY dz 
EL @ OMINDSeeO00— TPLE 2. HEAWY do 
209 Sia+GS REAUY 9% ' 
G s ON css cs oe ; Reavy 35 
#10 DU E12 flew KEADY Jo 
5 POUR!) MRAP O TN Fast Dh USEMRIGOIIE TT Biases: mictenninnienerenecanenecsteenetscnnineonutanene tb sane Me anata sen _WeADY 37... .. 
rl READY 96 
212 Mis? ReALY 39 
2h4 CALL FONG Ny GPe FP) XP Md) HEALY 40 
IFUITERQM.EU.O) GO TO 224 Reavy +i 
dO LOOT # Len REAUY 42 
10 MMA AN) PME! teen wes: <n raeich snmeneemnconeness co aieén nittin vin poets seis Selanne mnt to on versie « ROMY A... 
ou Ty shinai REAUY 44 
G et eee i en ice Desss vee  e HEAUY 42 
214 Mah AtMUy 4d 
+ 215 CALL MATMP [AyNeSeGPy GSP) - : ae Keauy #7 
216 UF (GSPI2IT 229.229 ALAUY 43 
‘yet vii AT wit MEPOE ) ZLO plz 9, 229. st ei na em RLAUY Y.. Ererey 
2hG WRITECGadd KeAUY 5D 
“i: (5: RUE OER ; : KeAUY 52 
bd ged eaLe 294 ao Reavy $< 
Lo Et, = 2.00981 rig ee prscre: SS re Z Reavy 39 
buy T9 210 , REAUY 54 
vai RIL. Fasre BMI Te See eet om eceiiaanse ns csateeos tap beeseoon Sore 1k hit eta hei cone wae, READY 9d 
232 00 294 Tahen ; KEALY 5a 
23% 4 ONT) eGPCT) 7 KeAuY 57 
23% TUDPeXPKLD ’ RLAUY ba 
RA VE Cob eeOV 22 2ede 223 ALAUY >Y 
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IMPLICIT REAL #@(A-Hy0-2) REUUC 4 i 
SUBRUUT INE REDUCEICUR sKREDs LAEDeLRED eWRF)--- ---.. Eo: KevuG 3 ‘ i 
CUMMON /IKCOMS XCOMI200), [COM1200) Revue 4 I 
CUMMUN OCCURSNUCCUR . ere iG Reus 5 | 
COMMUN /IDNUS/ I1D1L050)> 102450) Kivu 6 ' 
ve COMMONS LOPT/LPRUC 2 TN» NCONS » [PNT »LEX+LEMI T+ DRAND... vo atbad tots heres iw eae RE GUGE Fiza i 
CUMMGN “ENDS ITERM abvul 8 i 
DIMENSION CUP 20) ; ReEbuUG  Y | Z| 
OLMENSTON OCCUR 4000) sNOCCURL 30) RELUC 10 ' 
t : : hit se HEUUL Ld 
C : RELUL le ull 
£. . SUBROUTINE REDUCE. CAN .AULTIPLY -ANY. SPECIFIED ELEMENT .OF..THE — RELUL-435-.—.. m 
t CCLUR ARRAY BY THE REUVUCTION FACTGR WRF AND SET THE FIRST REGUS 14 | 
t ELEMENT OF THE CUP ARRAY EQUAL TO ANY SPECIFIED ELEMENT OF THE kKEbuC 15 | | 
C - OCCUR ARHAY MULTIPLIED BY THE FACTUR WRF. REQUC 46 ll 
C i RELUL AT | 
t REvUL 18 
‘cc S webataaeet Jee, +), Soreness decease ne ia aur dositae pY.aede sce SSN Saad «cate ae bi CNG: 
12 = 102013 REDUL 20 
IRED = {REO Ga 4eOue ¢). H|| 
IFCTRED.EQeLRED) GO ‘To 100 REUUL 22 iW 
IFCLTERM) 10 50, 60 ; Reuuc 23 | | 
c REVUC 24 M 
are ey eee Reet ae Ee ang tee ee a nt Oe ee en | q 
10 on 20 {| = 1, IN HELUL 2o 
LFCIOLAT) -EQ01Z) GO Tu 40 =. RELUL 27 
20 CONTINUE REDUL ch ‘ 
¢ RevUL 29 
TFC ICUM(3).EQe1) 59 TO 40 . RELUL 30 
Ges re a ce coccinea Aine amen eae means. KEV Ue obec 
40 “uccurstz) a WRF #OCCUR( IZ) RELUC 3¢ 
KRED # 0 erat tariaae fet z chae REVUL 33 | 
G0 Td 200 AcLUL 34 * 
¢ nthe. . : RELUC 35 
“#0 = CUPEL) = WOE FOCCHRET) REDUL Jo 
» BRED Oo. WK eM eSetea ce acts aelaee COMI as ES: Vay ete bake lw (oie RE OGL 
GO TO 200 ; KEDUC 3b 
C ook . . ae Rebus 89 
$0 KRED = } REDUC 40 
G 50 MAY BE USED FOA.LATER IMPLEMENTATION . : 7 i REVUL 44 
t REUUC 42 
(eg He Ds ae ae ee a EE eer i, ee ie Co RED UG a. | 
Go ‘10 200 REVUC 44 }) 
C s ning Soa | ie ee nes eee : Be otek RE ee a ae REQUE 45> i 
60 KRED | ; REVUC 46 
Pe OCR eee a Ree een REGUL 47. iq 
a REUUC 4&6 I 
8 BEES RU |S et oe ease EE oC ROI 2 Bae | 
100 KRED # Te - REOUC $0 
t. EBS fee nips ae ae ep Re ne, peor Pa fe er : REUUL 34 
200 RETURN REDUG 5¢é a 
’ END ee Ee A eres ein one t ° | AES i 
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IMPLECIT REAL #8 (A-Hy OZ) RITOD 2 
=» BUUROUT INE -RITOUT-—— n-ne me ee ee cr ee mt cere we oRITOU. 3 
c RITWU 4 
wees = DIMENSION ALMAX(200 by ALMEN( 200) ¢ TAMAX( 200) 5 TAMING 200) eZMAX(200), ALTOU . 5 
LEMENEZ00) ¢FMAN( ZOO) SHANE 2001s GceuR oven NOCC URI 30) alTuu 6 
Shy ew ae eae ss a Ne SIRE Ra eh ARN 1M Res We as ge eee ea aie ececns ert jemnrp eter RNID, MOU md iii 
EQUIVALENCE KITUU 8 
-- LEOCCURC1044) pALMAXC LD) sCOCCURLL2Z44) sALMING LID. «| - Ritus 9 . 
ZLQCCURC L444) e TAMAXE LY) pCOCCUR(LO44)  TAMIN( LI) s RIFUJ LO 
~ PEOCCURL L444) eZMAXCLD) » (OCCURC 2064) — ZMING LD)» RiTOU Li.. 
GUDECURY Ce SAT rE MAREE) Ts (OCCUR 2Z444)—¢ FMINCG1)) Ritu be 
wees eee FQOUEVALENCE : : PR ee ess mnt, eM Vice Ge fe 
LENUCCURTLILD IKMAX) 5 {NOCEURC 2)» LKMIN) e (NOCCUR (14) pNPRINT) AdTQU Ll 
RITUN 45 
CGMMCN Occur, NUCCUR RITQU Lo 
oe On cinta ates ee cee REVO 17 a. 
C RITOU Lo 
ve omsecuee see ~ SUBROUTINE..RITOUTPRINTS-.OUT. A-SUMMARYOF THE MAXIMUMS. AND... XITU 19 .. 
C MINIMUMS IN ANGLE OF ATTACK AS DETERMINED BY TRE TESTING IN RIiTouU 20 
+ C---- SUBROUTINE VIXEN e—-- - RITQU 2h 
C RITuu Zé 
LFANPAINT SE QeO)RETURN - aAlTOU Zo 
51 WRITE (651051) nlTU 24 
~~ LO5L FORMATO AHO SOX%423HMAX [MUMS—IN.. ALPHA. PRIME/40X%— @HTIME SIX elHZ oh 2k—...kiTuU 25 
L3HSMFeOXeLLHALPHA PRIME) KITOU 2a 
-¢ : snes RITUU «Tf 
OO 52 J*iey IT KMAX RITOU <b 
ALMAKOJ)=597 e295 7800SAL MAXI) - - RITOU 29 
WRETE(C Op LOSZ)TAMAX( J) p ZMAXU SI 2 FMAK(J)¢ALMAXC SD RITOU 20 
2052 FURMATILH 2UKeF Teds ZLEL4 5p 4X pF G03) . - RITUu 3h 
52 CONTINUE KiTuy se 
WRITE G,1053) RITUU 33 
x22 FURMATCLHO 50X,23H4INIMUMS I IN ALPHA PRIME/ 40K. GATIMEs OX SANZ oL2Ks RITO 34 
~ - LOHSMF e 6X_ IL LHAL PHA. PRIME)D-- vein oF tetas 9 RITOU 35. 
c RITOU 36 
pO 53. JeL eI KMIN af 2. hs sda bake eam eae atyas tons on nate ora ea OL apeiee: 
ALMIN( J 1576295 7BVO*ALMIN(J) RIT 38 
. SRETCLE, LOSEITARING JVs AMINCA) eRHEN( OI SALMENGS) RIT 39 
53 CONTINUE ¢ RiTUs 40 
- RETURN - pede arts naehecn eae eet me RITUU 41 
END RITUU 424 
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CKHOSURK, » a nme 1 tre a nee rene SINELAN Det ine ILE Me a ls gh mt IE Ew ROSHR 4... 
SUBROUTINE ROSORK RUSHK ¢ 
LMOLICEYT WEALS@B (AwHy Dod) oe ee ee ne cee ROSBR 3 < 
CUMMON OCCUR, NOCCUR ADSHR 4 
COAMON/BLKOSHHEI1600) pF E40 be D140 sELA OD eDPLAOD AIHODs — - hUSbR > 
L 0403465265, dh | ei ae hoa last TeRSe SL 9222 KUSHR 
MOLL AXX EL i ere ee cb wat ren oom act MUS BRS S0h soe 
QDELAK+ FAK sCXKL40510) HUSBK 8 
GuhIMDN YS OOPT SEO 20) KOSUR 9 
CMMON/ENDS I TERM KUSBK LO 
COMMONS LOPT/ LPROC) VMAX pNCUNS» EPNTs LEXeLIMET SL RAND KUSHK Ld 
GOMMUN /IXCUM/ XCOM{(200)5 peOnteoe RUSBR bd 
UIMENSLOM CO2Z0,20)...- Deo Shei eA oii wid ten Poa ee am Cee passe MUS BELA. cas 
QIMENSTON GECURLSOGG) » NOCCUR( 30) RUSbM Le 
EQUIVALENCE (HHiGOL)eC (add KUSBN 45 
EQUIVALENCE (XCOMN(L)» ALPHA) RUSoK lo 
EWUIVALENCE (XCUM(2 9, BETA) KUSBK 17 
EWUIVALENCE (XCUM(3 35 GAMMA} KUSBR Ld 
» EQURVAGENCE (XCOM(4°S) DEL)... te ct ee | ROD OK 19. 
EQUIVALENCE (XCOM(I5)¢ RATU) KUSBR 20 
EQUIVALENCE (XCOM(6)-TOL) HUSBR od 
€ KOSHK 22 
¢ ae ee ; : . . .* KUSBK 24 
G SUGROUTINE ROSBRK ATTEMPTS TO MINIMIZE THE VALUE OF A FUNC ROSuK £4 
Come ames THON UE Plo PZo0eee PU NMAX)) WHERE. THE P'S. ARE. ADJUSTABLE. owhUS BK 25 2 
C GESIGN PARAMETERS BY THE JOSENSROCK RUTATING eRAR ENE MINI- RUSUK 26 
G MIZATION TECHNIQUE « : ROSBR 27 
C x > RUSBR 2s “ 
¢ NMAX IS THE NUMBER OF PARAMETERS TO BE ADJUSTED. ITS MAXIMUM " RUSDR 29 
¢ VALUE IS 20. ROSBK 30 
Goi Basta Pee Si a by Sellen Wes eaeon Seen 5h Pas < ROSON 34... 
C ; RUS BK 22 
NIKTA = Q ene . ene KOSHRK 23. 
NSTAG = 0 RUSBR 34 ) 
INDEX = 0 eR eee AUSBR 35 
¢ ‘ RUSHK 36 
VSB -= GALL FE VO NMAK 5 gh) bea cee eet teste te ee tame cn nee mnt tie ae eae aaa RUS Yt eee 
TFUITERM.NE.O) GU TO 600 ROSHR 34 
* tae We egiena: Pe att , RUSHK 39 
LE(NTREA.LEsLIMIT) GO TO 610 | RuSsK 40 
, WRITEUGebLLF -- RUSBR 42 
ont FORMATE LHO 34HNTKIA GREATER THAN LIMIT IN ROSBRK) RUSHR 42 
POUIG IMME MEDAN poet cn aN rrekscc ae sb ous cathe enmmmentecsoranentiaichene pence = eh un cesi= Tass cnt dlen Aonieeiemn no AUN BK GS Ph. 
Gu TO £20 RUSUK 44 : 
G10 =TNDEX = INDEXGL RUSH 45 
¢ KUSUK 40 
be IF(NTRIASNE.O) GO TO BOS Waser RUSBR 47 
c RUSHR 4b 
be ULAS Tai .. a cc gin URE eG 
USTAG = U RUSUK 50 
ER, eer APART STBND me cis nce’ me rcoatlates mpeenbeinnericiuinn fons naeinedd of emai secesier kesh KUSHR 31 
NSUCC # 0 ROSBR 22 
ONTRIA @ LL. eS .:. | een eee eee RUSHR b3 
DU 104% t= 2yNHAX : KUSKA 54 
BAMA E (CNOND RN ONL Ps aca ee ce eee tener ee cma Plas cig ce Sereda ae aeteee SUS UN De 
¢ ; RUSUK So 
G C MATRIX LS SET TO A UNIT MATRIX KUSUK oF 
, i RUSK Sd 
NU 105 1*d_—NMAX KUSUK 5Y 
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OU 105 dale NMAX HOSUR 60 
LEC LI) 1069 ROT 9 2 Ob amen a— mem cents ce ee ROSUA 61. 
106 Cited) # 0,00 RUSH of 
-Gu TU 105 | mth 2A Bona meanest ie PS Rae de aaa pte ee HOSGH 64 
10oT Clled) * 1.00 KUSGH 64 
105 CUNTINUE Ren eee re ese RUSUA 6d 
Nel RUSBA 66 
Se (VOL VSaO0 GO oats eetecatecaacis a sane bare Seas trans wp ee oie oan aunt ea ats Se ee SER ao tiee 
GU TO 108 RUSUK 64 
Cc nen en RUS RK 69. 
C SECOND AND SUBSEQUENT TRIALS, START HERE RUSUK 70 
C 6 KUSHR Th { 
‘ RUSUK TZ | 
ROD MTR VA RANT REALE A.D) cence crteer te aioe seein nies iin ihe eaten seinen etn RUSUK 72 ._. 4 
LFIU-ULAST#(L.O0ETOL)) 109, 109110 KUSBR 74 ij 
C RUSUR 75 |! 
C LF TRIAL nas A SOEGESS» START | HERE RUSHR 7b i 
C a Setimiece = Aste KUSUK TT. | 
109 ULAST © U RUSBR To | 
Fee IN SUC Come GUC Gs Bre Dey ces aes se ee ee eo ee RUSE Th $4 
O(N) = DCN) & ECN) KUSER 80 i I 
TFINSUCC.GE+20)- GO TO. 50 ee KUSBRK do] 
E(N} = AL PHA#E(N) KUSUR 82 | 
- GU TO 108 5 se He me on RUSUK 3 yi 
c RUSHRK 64 | 
Gwe TE. TRIAL WAS A—~RATLUR Ep START-HERE-—— 2. ~~ ee ee ee RUS UR OS, i | 
110 CUNT INUE KOSUK Bo 1 
WRITEL 62 306INTR IAN NSTAG NSUCC eUeNMAX o (PUT), IT ab pNMAX) ROSaK wT 
a 3060 FORMATCLHO2Z4H NTRIA N NSTAG NSUCCLOXLHU/ LH ¢165F20. LO/iH RUSHER ot 
1 13HPCT) f= .2- TO 134th 6F20.10)) . RUSOK Bg 
OO LiL T=hyNMAX ROSdK 90 
LL PUL) = POTD = DPCT). --- ~~: RUSUK GL. — 
U = ULAST RUSUK 92 
TFOENSUCC) Li2eLizekl3 RUSBK 93 
’ 112 EIN) = -dETA#E(ND RUSUK 94 
~ G0 TO 125 ee ete ee we oe cn ee - RUSUR 95 
"$0 wWKITE(6,922) NSUCC ROSHK Go 
922 FORMATOLIHO LSHWARNING-NSUCC2. 1691X-22HGOING.TO NEXT. VARIABLE)... .RUSDR 97... 
c ; RUSUR 98 
G SEARCH IN ONE OIRECTION ---- ---- RUSUR 99. 
CG KUSBA190 
1L3 EIN) = GAMMABDIN) -.- - meee : KUSBRIOL 
C IF N = NMAXy THIS STAGE IS COMPLETE. AUSBKLI 
+ DON EQ NMA GO TOD LG ne a ee ee ee emcee KOS BRLOS WW. 
N@Nal RUSBRLI4 
ETRY me I on oa cn arene eames teers ie see tmtatatetion ee tee sensi ee ne ktm | meat - . ROSHRLOS 
NSUCC = O RUSHRLOG 
; GU TOMB 2 56 Plan ec orate a toa eS PLD . <Ao NE Ronee ote MUSBRLOT 
Cc RUS BKRLOY 
Cen meems-- STOPPING.-CRITERIALS-APPLIEO AFTER 2ND.AND.SUBSEQUENT_STAGES... RUS HKLOY .._. 
Cc KUSHKLLO 
SNL TR UMSRAG) RMR ie PUL pee a ics as ncn awa ee wl WUSBKLLI .... 
117 ULIM = UPREV = DEL*®UPREV AUSdKade 
LFC(UPREV“USTAG)6EQ.0000) GU TO.L2k. 2... . « KUSERILS 
URAT = {USTAG@ULASI)/(UPREV=USTAG) KUSoK 4b) 4 
TFCULAST —@ ULTMIL 2b L2L gL LG. me me ea ee ee ee eae nee KOSokLa5 . 
Vis LFCUSTAG = ULIM)L2Le1219119 : KUSUK ALD 
LL9 TFCURAT © KATUDL2Z00121 9121 KUSuKLE? 
C ; RUSBKLAY 
120 WKITE(6,301) KUSURILY 
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301 FORMAT( AHO 2GHROSBRK SEARCH 1S COMPLETE.) 
-WRETE (64302) ULASTs (RADI) Eel @NMAXD-— ees ind at seit. 
302 FORMAT(LHO37HEINAL PERFORMANCE FUNCTION. ULAST = E20.10/1H 
L2OHEINAL PARAMETERS ARE /(1H 6€20.10)) yo 
WRITE (09303 INTRIASNSTAG 
403 FURMAT(LHOLSHNUMBER OF TRIALS-= 15/1H 1BHNUMBER OF STAGES = 15) 
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C SEARCH 1S COMPLETE 
c 
© SEARCH WILL CONTINUE IN ANEW. COORDINATE. SYSTEM... 
121 UPREV = USTAG ; 
VL STAG BULA GT eon meee man tee ne me creer re mtn creme nme summer 
NSTAG # NSTAGEL 
¢ : sierra: . 
C DEFINE NEW COURDINATE SYSTEM. 
CALL GRAM eee Es 
WRITE( 65304) 
304. FORMAT (AHO. 20HC-MATRIX. FROM _ROSBRKA—.).———~ ---— a ri ee 
pu 200 Tal,NMAX 
200 WAITEL 69 205) (COL ed de Jab eNMAX) 
305 FORMATCLZ(1X,1PE9 02d) 
129 Nel ee. <3 + 
bird 0.00 
- NSUCCBO oo nee eee meee merece nn cennememe ts meneame cemnnnee ams ms amuetanat mammnnas ment 
108 CUNTINUE 
WRITE (Op 306)NTRIA PN NSTAGSNSUCCOUsNMAXs(P CLD e Tale NMAXD .. 
G 
C DEFINE CUURDINATES FOR NEXT TRIAL. 
125 00 122 [=LsNMAX 
OPEL PAC OL eNPFELND. © coe oe Pete ied. poeta Seamer eh a Aion, deteids Aamtte 
L222 PUL b@PCLEDPUT) 
C 
WRITE (6,500)(DP U1), Lal sNMAX) 
WKETECGsSOLINMAXS (PSTD eleLyNMAX).. ©. . = 
WRIETE(6¢503) Ne EWN) 
500 FURMATULHOL4HIN ROSORK OP #/f2H 1P6E20~7)) POE Perret a We aoe 2 Meee 
301  PUMMATCLIHOL3HPC ID Eo = 2 TO 13/tRH 1P6E20.7)) 
303 FORMATELAHO 13HIN ROSURK N # 13s2XOHE(N) =£20.10) 
| GO TO 123 
400 RETURN 
END 
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COROT ATE mn nme mnnnnennn eet hc a ci leek ats pata acts tacaneaate & late Me co 
IMPLICIT REAL #B (A-H,O=2) KUTAT 2 
oF . - SUURUUT LNE- ROTATE (DER IV yh P}—-—--— --— = 2 eee ee ee HUT AT - ee 
C AUTAT 4 
REAL 2B MX gM Vg MZ mem mame ma ne ce te cemeeceems ce eee os ROTAT 
REAL#3 NUMERL» NUMER 2 ROTAT 6 
BUC ees ose atic ee Sennen eee ene Sassanid eee en ORO AL cote 
DIMENSTON TABIX(50) pTABI(50)) TAUZE SO)» ‘OCCURL+000) +NOCCURI 30). RUTAT 8 
1 DERIViLSOd ¢XJRIZ5) eXJTI25) - — -- ... .RUCAT.. 
DIMENSION XNRG25)— XN1I25) RUTAT LO 
C - es : : RUTAD da 
EQUIVAL ENCE . KOTAT L2 
we veem DEUCCUR (O01) pAREF.-) 9 OCCUR(O02) pALPHA. Je (OCCIRICO3) sALPRIMDe.. - RUTAT 13.0... 
ZLUCCUR(OLZ) sCNALP }y€OCCURLOL3) +CMALP )s(OCCURIOLS) CD de RUTAT 14% 
30ICCUR (O20) ,CMQ be OCCUR(O21).0 PV, LOCCUR( O27) G6 de RUTAT 15 
4(UCCUR( O42) 5Pl yf UCCURIO43)—P },(UCCURL044) -PHI de RUTAT 16 
- SC OCCUR (O45) gPSI.- bef OCCURIN5S0) 49 )eCUCCURCO51) 900 de ~ RUTAT LT. 
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S{OCCUR( 2373 5XUP De * ROTAT 2B 
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C L2PT = 2 IS A SIMPLIFIED ANGLE OF ATTACK SOLUTION REQUIRING NO AUTAT 42 
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i COSLaCOSI THEALPS ROTAT 72 
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By INS ei uen Wes Jone pawns Ovan “t PNET LPH RTRNA HS WR) TENS SORNE/ 11 RUT AT ALO ll 
ve mentee OBR IVER SD RM MAUD Koco ee eee KUTATIL4 i) 
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C SUMPLIFIED “ANGLE OF ATTACK SOLUTION WHICH REQUIRES NO RUTAT Led 
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10 2.0U0#We(BETAZ*S INGO) ¥¥2) RUTAT 140 
1? CONTINUE KUT AT A4%4 
THAR==2.000#P1#2/(SETAL#V*SINGO*DSOR TL AK2*RHUINF)) RUTAT b4e 
[FCTHARLCTECONIGU TO Le RUTATL43 
IF(TBHAR.LT.1.00-10)GU TO Lé ROT AT L44 
Bi Ober ee aul eaa ae eee RUTATIGO5 .... 
¢ OSCELLATORY SOLUTION FUR ‘CYCLE TIMES, TBARy GREATER. THAN — RUTAT L400 
C TECON. RUTATA47 
‘ SMF21L.0D0/TGAR RUTAT LSS 
TEM=2 .QU0#USUKT (CAKE AK2) ®KHOINF) KUTATL49 
CALL HESSELCTEMs0.0U0,1L.00-S5p1L0eXIRe XI edd KUL AT ad0 
CALL NEUMAN(1 OU Se TEMyOeUD0 210s LeOD Ae XNRe ANE XJVCL I oKILUL) XIRC dul AT Lon 
LexJ1f 2091) RULATide 
IF(Z.NE.20)60 TO z KUT ATaba 
# UTEMUZ==0 .S5008TEM*HASETAL/Z ROVATL5S4 
VENACZZOTEMUZ#(HANRI 2D FXIRCLIEXIRIZISXNRGL DD. RUTAT L355 
NUMERL = (XIRULI *AKLFTHEALO*2 » 37690~3*BETAZ) /(Z9 JE PCAKL®RHOINE & RUlaTL5o6 
-ABETAZI) ce ee ~ RUDATIOT 
NUMER22 THEALO*DTEMDZ* (=X4R (2) ®DEXP(=AKLERHOLNF) KOT AT be 
AC2=(NUMERL-NUMEK2) /DENAC2 . RUTATLOS . 
ACL= (THEALO#DEXP( AKL ®RHUINE) -AC2*XNRULD)/XURULI RUTAT ILO 
soe ALALTSDEXPC AKL *RHOINF) /OSQRTIPL*OSQR TC AKZ*RHOINF) } KUT AT Loe 
2 CUNTINUE RUTAT Ld 
mene ALPKIMEDEXPCAKLSRHJINE SOC ACLEXIROLIGAC2#XNRGLI) | RUT AT Lb 
CN=CNALP*ALPRIM KOTAT 104 
~- - CME CHAL PAL PRIM .n ene eee ; soll tp hating den (ete: - RUTATLOS .... 
ALPHA=ALPRIM RUTATi oo 
; RETURN : SMe Was sano a3 sac oe RUT AT 467 
C RUTATL 66 
wie! sce . ENVELOPE..ANGLE OF ATTACK CALCULATION FUR TBAR LESS THAN.OR_ KUTATL09 . 
C EQUAL TO TECON. ROTATE? 
- 16 AL AL O=DEXPCAKL #RHUINF J /OSQATAP I ®DSOR TC AKZ*RHOINE) ) ROTATALAL 
LFC2.bQeZ0)ALALI*ALALO RUTAT AZZ 
- ALPKIM#0.63601L00*THEALC#ALALO ... jo rinc igs gets RUT AT ATS 
ALPRIM=ALPRIA/ALALI RUTATAT4 
- SMF =0.000 : By cy apt uae G team tees Manca shea how aa son imeasale Zacceitya chal Sacer a heed cmemen HUTATLTD .... 
C ; Fs : KUT ATaTo 
u IF ANGLE OF ATTACK ENVELOPE VALLE IS LESS THAN OR EQUAL TU RUTATLT7 
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IMPLICIT REAL*8( AH O02) 
. SUBRUUTINE SAVEUVIK pOVARS OCCUR) — - — --- ~ 
DIMENSION OVAR( 50), UCCURI4000) 


GCCUR(205~209)— AND OCCUR(222) IN LOCATIONS D¥AR(L<19) 
- RESPECTIVELY ANO LATER RESTORES THESE VALUES TO THE UCCUR 
ARRAY. 
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RETURN 
200 - OU e221 = ly 13 - 
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APPENDIX II - Overlay Information and Other Pertinent 
Information for running the Optimum Decoy 
Design Program 
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SET UP OF BINARIES FOR THE OVERLAY 


Program 2542 has been run successfully on the OS65 at Avco using 
version 14 of the operating system. 


All binaries whose names do not follow one of the OVERLAY cards 
should be first in order in the binary deck, 


These binaries should then be followed by the following six overlay 
lea ; segments, The names of the binaries included in a particular segment 
| follow the overlry card for that segment, 
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| i,0 A Preset Deck is needed to run the program. 


2,0 The subroutine LA000000 is written in Assembly Language. 


ht | 3,0 The input quantity 10P(76) should be set to 1 at AVCO and 
| to 0 at Aerospace. | 


The dummy subroutine PLT must be removed at Aerospace. 


The dummy subroutines ENDJOB, EZPLOT, FRAMEV, | 
IDFRMY, and PLTND must be removed at AVCO, ‘’ | 


